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FOREWORD 


This  dissertation  describes  digitaJ  logic  operations  and  reconfigurable 
optical  interconnects  realized  from  aluminum  gallium  arsenide  (AIGaAs) 
eiectfXK)ptic  Fredkin  gates.  The  Fredkin  gates  f^ricated  for  this  effort  were 
channel  waveguide  type  passive  and  active  AIGaAs  zero-gap  directional 
ooupleis  (ZGDCs).  The  first  AND  logic  gate  reported  in  the  open  literature 
using  AIGaAs  ZGDOs  is  demonstrated.  A  reconfigurable  optical 
imerconnect  is  demonstrated.  The  theory  behind  the  Fredkin  gate  and  the 
ZGDC  is  presented.  Several  models  predicting  the  refractive  index  of 
AIGaAs  versus  Al  oonoentration  and  wavelength  are  treated.  The  electro¬ 
optic  properties  of  AIGaAs  are  also  presented.  A  riKxJel  for  predicting  the 
electro-optic  coefficient  of  AIGaAs  versus  Al  concentration  and  wavelengih 
is  developed.  Models  for  predicting  the  interaction  length  and  behavior  of 
ZGDCs,  as  wen  as  the  switching  voltage,  are  presented  and  optimized  for 
AIGaAs.  Fabrication  techniques  for  AIGaAs  ZGDCs  are  formulated  and 
analyzed.  Device  characterization  techniques  are  presented.  The 
performance  of  the  actual  AIGaAs  ZGDCs  fabricated  for  this  effort  fit  the 


III 


models  exceptionally  well.  Fabrication  tolerances  are  also  addressed  for 
both  passive  arxl  active  devices,  issues  for  integration  of  AiQaAs  ZGDC’s 
into  practicai  circuits  are  addressed. 

Correlation  of  the  experimental  and  theoretical  results  yielded  a  mole 
fraction  of  Ai  of  30%  for  the  core  layer  and  35%  for  the  cladding  layers. 

Refractive  indices  of  3.420  ±  0.003  for  30%  Ai  and  3.368  ±  0.003  for  35% 
Ai  at  a  wavelength  X  -  827  nm  were  determined.  A  linear  electro-optic 
coefficient  of  *  -1.49  ±  0.08  x  10*’*  m/V  at  X  =  827  nm  was  inferred  for 

the  core  layer.  A  switching  voltage  of  7.6  VDC  ±  0.3  VDC  was  achieved  for 
an  AIQaAs  ZQDC  with  a  5000  pm  interaction  channel  length.  This  is  the 
lowest  switching  voRage  reported  in  the  open  literature  for  an  active  AI^Qa^. 

xAs  ZQDC  with  a  5000  pm  interaction  channel  length.  Propagation  losses 
measured  0.5  - 1  dB/cm.  Device  losses  measured  3  •  7  dB. 
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CHAPTER  I 


INTRODUCTION 

As  electronic  and  optical  data  processing  circuits  approach  multi- 
gigahertz  operation  rates,  the  need  for  mutti-gigahertz  communication 
between  these  processing  circuits  becomes  essential.  This  communication 
ndudes  cNp-to-chip  on  a  common  multichip  module,  multichip  moduie-to- 
multichip  module  on  a  common  board  and  board-to-board  through  a 
common  backplane.  The  data  rates  at  which  communication  through 
electrical  imeroonneds  can  be  transmitted  is  limited  to  the  low  gigahertz 
rates,  unless  one  is  willing  to  expend  large  amounts  of  power.  This  results 
from  the  fact  that  electrical  wires  arxl  striplines  posses  frequency  dependent 
dissipative  fosses  due  to  conductive  and  dielectric  properties,  and  frequency 
dependent  parasitic  coupling  fosses  due  to  radiative  and  surface  wave 
propagation  properties.^  The  higher  the  frequency,  the  higher  the  losses. 
This  is  where  optics  can  provide  solutions.  Opticai  interconnects  can 
transmit  data  at  opticai  frequencies  with  orders  of  magnitude  lower  losses 
than  their  electricai  counterparts.  The  power  requirements  of  opticai  and 


1 


electrical  Intercormects  versus  data  rate  is  illustrated  in  Fig.  1 .0.1 .  Fig.  1 .0.1 
was  generated  using  the  following:  Waveguide  lengths  for  both  the 
electronic  and  optical  interconnects  were  50  mm.  The  total  number  of 
interoonnecis  was  600  for  both  types.  For  electrical  interconnects  on  an 

alumina  substrate,  waveguide  widths  of  25  pm  and  waveguide  heights  of 
650  pm  were  assumed.  The  effective  and  relative  permittivities  of  the 
microstrip  used  were  6.7  and  9.9,  respectively.'  The  capacitance  used  per 

line  was  6  pF  and  an  impedance  of  50  was  assumed.'  The  power  versus 
data  rate  for  the  electrical  interconnects  was  determined  with  the  methods 
given  by  Edwards.'  For  optical  interconnects  on  a  gallium  arsenide  (GaAs) 

substrate,  waveguide  widths  of  3  pm  and  waveguide  heights  of  1  pm  were 
assumed.  Waveguide  losses  of  0.2  dB/cm  were  used  at  1  GHz  growing  as 

0.06  dB/QHz.**^'^  The  modulator  loss  used  was  5  x  10*  W/channel  at  1 
GHz  growing  as  5  x10*W/GHz.'*^  The  detector  capacHanoe  used  was  0.3 
pF  and  the  detector  amplifier  capacitanoe  used  was  1.7  pF.'*^  The  power 
used,  per  detector,  was  1  x  10*  W/channei  at  1  GHz  growing  as  1  x  lO** 
W/GHz.'*^  The  laser  power  used  was  100  mW  at  1  GHz  lessening  at  the 
same  rate  as  the  waveguide  loss  per  GHz.'*^  Finally,  20%  electrical 
efficiency  and  50%  optical  coupling  efficiency  was  assumed  for  the  laser. 
The  power  versus  data  rate  was  then  determined  by  adding  up  the  power 
usage  of  the  separate  components.'*^ 
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Data  Rate  (GHz) 


Figure  1.0.1 

Power  Requirements  for  Optical  and  Electrical 
Interconnects  versus  Data  Rate 
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In  adcfition  to  the  lower  power  consumption  and  higher  bandwidth, 
optical  interconnects  offer  isolation  from  electro-magnetic  interference  (EMI).^ 
Optical  interconnects  are  starting  to  be  utilized  in  chassis-to-chassis  and 
board-to-board  communication  using  glass  optical  fS)ers/^  On  the  smaller 
scale  of  multichip  modules  and  within-chip,  integrated  optics  is  just  beginning 
to  emerge  Integrated  optics  offers  the  addHIonal 

advantages  of  smaller  interconnect  widths  and  in-plane  crossover  capability 
with  very  low  crosstalc.  This  eliminates  the  need  for  air  bridges.  Thus,  less 
wafor  and  board  real  estate  wiH  be  required  for  optical  interconnections, 
which  provides  room  for  additional  circuits  and  devices.  Fabrication 
processes  for  integrated  optics  must  be  oompatibie  with  those  used  for  the 
manufacture  of  electronic  interconnects  and  devices.  This  is  important  from 
the  standpoints  of  practicality  and  acceptability.  Aluminum  gallium  arsenide 
(AIQaAs),  as  wen  as  other  lll-V  semiconductors,  is  an  attractive  material 
because  one  can  fabricate  not  only  electronic  integrated  circuits  in  it,  but 
also  laser  sources,  detectors,  opto-eiectronic  integrated  drcuits  and  optical 
waveguides.  AIQaAs  also  possesses  electro-optical  effects  which  can  be 
utilized  to  fabricate  optical  switches. 

Optical  switches  can  be  incorporated  with  optical  waveguides  and  other 
optical  switches  to  realize  digital  logic  operations  or,  very  importantly, 
reoonfigurable  interconnects.  Current  research  to  develop  processing 
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techniques  for  AIQaAs  optical  waveguides  and  devices  has  shown  that  the 
processing  techniques  for  these  optical  circuits  are  compatible  with  the 


techniques  used  for  the  fabrication  of  electronic  circuitry  in 
Therefore,  by  using  AIQaAs  and  GaAs,  one  can  house  the  entire  electronic/ 
opto-eiectronic  circuit  on  a  sirtgle  substrate.  With  the  incorporation  of 
AIQaAs  and  GaAs  opto-eiectronic  circuits  and  optical  waveguides,  the 
muitichip  module  or  board,  now  entirely  electronic,  could  perform  more 
efficiently  with  less  power  consumption,  move  data  around  faster  with  less 
noise,  and  use  less  real  estate. 

An  even  more  attractive  benefit  is  the  potential  to  reconfigure  the 
multichip  modules  to  either  perform  entirely  different  operations,  or  to  route 
signals  around  inoperative  or  defective  portions  of  the  module.  It  is  an  area 
that  the  United  States  Air  Force  (USAF)  and  the  Advanced  Research 
Profects  Agency  (ARPA)  is  very  much  interested  in  pursuing.  There  are 
several  ori-going  contractual  efforts  as  well  as  several  more  planned  to 
develop  integrated  optical  lnterconnects.*’^^*^’^®^-’®’^“  The  research  topic 
proposed  here  could  directly  benefit  or  reshape  such  programs. 

FredWniSile^  Conservative  Logic 

Fredkin  gate  is  a  term  given  to  a  nonlinear  dual  input/dual  output  signal 
processing  gate  that  performs  a  conditionai  crossover  of  the  two  input 
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signab  via  an  input  controi  signal*  (see  Fig.  1 .1 .1).  When  the  input  control 
signai  is  a  iogic  iow,  the  two  input  signals  cross  over  as  illustrated  in  Fig. 
1.1.2.  When  the  input  control  signal  is  a  logic  high,  the  two  input  signals 
fbHow  a  paraliei  path  as  shown  in  Fig.  1 .1 .3.  The  logic  of  the  Fredkin  gate 
can,  thus,  be  represented  by  Table  1.1.1,  given  originally  by  Fredkin  and 
ToffbU.*  Fredkin  et  al.  presented  the  unit  wire  and  the  Fredkin  gate  to  form 
the  elements  on  which  consen/ative  logic  is  based.*  Fredkin  et  al.  defined 
conservative  logic  as  any  Boolean  function  that  is  invertible  and 
conservative.*  By  invertible,  it  is  argued  that  invariarKe  under  time  reversal 
is  not  required.  By  conservative,  it  is  stated  that  the  model  for  computation 
is  only  required  to  possess  at  least  one  additive  conserved  quantity.* 
Conserved  quantities,  such  as  conservation  of  energy,  momentum  and 
angular  momentum,  are  functions  of  the  system’s  state  that  remain  cor)stant. 
With  conservative  logic,  it  is  ideally  possible  to  build  sequential  circuits  that 
dissipate  no  internal  power.*  Where  the  two  input  NAND  gate  forms  the 
universal  primitive  for  the  set  of  all  Boolean  functions,  the  two  inputAwo 
oufout  Fredkin  gate  forms  the  signai  processing  universal  primitive  for  the 
set  of  aH  conservative  logic  functions.* 


Figure  1.1.3 

Operation  of  Fredkin  Gate  with  Control  Line  High 
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Table  1.1.1 

Representation  of  Fredkin  Gate  Logic 
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Dioital  Logic  and  the  Fredkin  Gate 

The  functions  we  seek  are  often  non-reversible.  These  include 
elementary  Boolean  functions,  such  as  AND  and  OR,  as  well  as  more 
complex  functions,  such  as  adders  and  flip  flops.^^^'^  By  supplying 
constants,  along  with  the  arguments  to  one  or  more  Fredkin  gates,  one  can 
realize  many  of  these  non*invertibie  functions  using  conservative  logic.  Rgs. 
1.2.1  >  1.2.3  show  conservative  logic  realizations  of  AND,  OR  and  NAND 
functions.  Fig.  1 .2.4  illustrates  a  conservative  logic  realization  of  a  J-K  flip 
flop.  The  unit  wire,  represented  in  Figs.  1 .2.3  and  1 .2.4  with  a  non-inverting 
buffer  symbol,  is  a  conservative  storage-transmission  primitive.  It  both 
stores  data  and  moves  data  from  one  point  of  space-time  to  another,  and  it 
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AND  Operation  using  Fredkin  Gate 


u 


V 


a  (don't  cars) 


1 

b 


V 

Vi 

V 

y2 

Xg 

a  -f  b 

9  +  b  (don't  care) 


Truth  Table 


Inputs 

U 

Outputs 

V  Yi  Yi 

\mm 

1 

0 

0 

0 

1 

1 

1 

0 

1 

1 

1 

0 

1 

1 

0 

1 

1 

1 

1 

1 

Figure  122 

OR  Operation  using  Fredkin  Qate 
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NAND  Operation  using  Fredkin  Gates 
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is  sspaisited  by  one  unit  of  time.^  The  unit  wire  connects  one  conservative 
iogic  gale  to  another  from  the  conventional  sense  that  the  gate  is  an 
instantaneous  element  and  the  unit  wire  is  the  delay  element  contained  in 
an  interconnection/  The  output  conserves  the  value  present  at  the  input, 
if,  for  every  new  computation,  we  supply  a  new  set  of  constants  and  throw 
away  the  *don1  cares*  present  at  the  output,  we  will  dissipate  energy.  For 
almoet  aN  Boolean  functions,  the  number  of  gates  increases  exponentially 
with  the  number  of  input  lines  and,  thus,  the  number  of  ‘don't  cares*  at  the 
output  also  increases  exponentialiy.^'^  This  results  in  an  exponential  growth 
in  the  dissipated  energy.  Conservative  logic  reduces  the  degree  of 
dissipaled  energy  growth  since  the  number  of  gates  does  not  depend  on  the 
number  of  inputs.  This  results  in  energy  dissipation  that  only  has  linear 
growth.^^^  The  reduction  in  the  dissipation  of  energy  is  achieved  with  the 
introduction  of  inverse  conservative  logic.^ 

Inverse  conservative  logic  is  attained  by  replacing  the  conservative  logic 
gate  with  its  inverse.^  The  Fredkin  gate  is  attractive  because  K  coincides 
witti  its  inverse.  The  inverse  looks  IBce  a  mirror  image  of  the  original.  Fig. 
1.2.5  shows  a  oonsenretive  logic  realization  of  a  logic  function  and  its 
inverse  using  Fredkin  gates.  Next,  by  joining  the  originai  function  with  its 
inverse,  and  then  connecting  the  outputs  of  the  originai  Fredkin  gate  network 
to  the  inputs  of  the  inverse  Fredkin  gate  network,  one  will  get  back  the 
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Conservative  Logic  Operation  and  its  inverse  using  Fredkin  Gates 
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orls^  aigumenis  and  oonatants/  Fig.  1.2.6  illustrates  this  property  by 
interfacing  a  NAND  operation  with  its  inverse.  As  can  be  seen,  the  ‘don’t 
cares*,  namely  Vj  and  y4  in  this  case,  are  also  used  as  inputs  to  the  inverse 
logic  gate  network  in  order  to  help  generate  the  orig'mai  arguments  and 
constant.  One  would  also  Nke  to  have  the  output  of  the  furx^tion  generated 
by  the  first  Fredkin  gate  network  because  it  contains  the  results  of  the 
desired  logic  function  you  want  to  realize  in  the  first  place.  So,  if  in  addition 
to  feeding  the  output  from  the  NANO  function  gating  network  in  Fig.  1.2.6 
inio  the  input  of  Ns  inverted  gating  network  it  is  also  fed  to  the  control  input 
of  an  additional  *Spy*  Fredkin  gate,  one  will  generate  not  only  the  original 
input  arguments  and  constants  at  V4,  y^,  and  y^,  but  wHI  also  output  the 
desired  logic  function  operation  at  y,.  This  garbageless  NAND  operation  is 
demonstrated  in  Fig.  1.2.7.  The  *Spy*  Fredkin  gate  acts  simply  as  a  fan  out. 
So,  even  though  we  may  not  wish  to  make  further  use  of  the  original  input 
arguments  and  constants  found  at  V4,  y7  and  y„  and  therefore  let  the  energy 
be  diesipeted  as  heat,  the  energy  dissipated  by  these  unused  lines  will  grow 
only  as  3N,  where  N  is  the  number  of  output  iines.^  As  stated  earlier,  the 
number  of  gates  is  independent  of  the  number  of  lines,  thus,  the  energy 
dissipation  grows  Knearly,  not  exponentially.  This  can  play  a  major  role  for 
large  circuits  with  many  input  lines. 
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Figure  1^.7 

Qarbageiess  NAND  Operation  using  Conservative  and  Inverse 
Conservadve  Logic  FrecHdn  Gates  with  'Spy*  Fredkin  Gate 


An  optical  FradMn  gate  is  the  opticai  realization  of  the  Fredkin  gate  with 
the  use  of  an  opticai  switch.  The  argutnent  channeis  of  the  optical  Fredkin 
gate  are  optical.  The  control  channel  for  these  devices,  however,  can  be 
either  opticai,  eiectricai,  or  acoustical.  Therefore,  the  optical  Fredkin  gate 
need  not  be  entirely  optical,  it  also  does  not  have  to  operate  exactly  as 
Fredkin  proposed.  WhHe  some  of  the  opticai  Fredkin  gate  designs  cross  the 
argument  channels  with  a  logic  low  control  channel,  as  Fredkin  et  al. 
proposed,  other  designs  cross  the  argument  channels  with  a  logic  high 
control  channel.  We  will  designate  this  second  type  of  gate  as  a  reverse 
Fredkin  gate,  which  is  not  the  same  as  an  inverse  Fredkin  gate. 

Optica}  Fredkin  gates  have  been  proposed  and  demonstrated  by  several 
researchers  utilizing  various  optical  switching  designs  and  materials. 
Proposed  optical  implementations  of  the  Fredkin  gate  include: 

1)  Polarization  switch:  The  polarizations  of  the  two  arguments  are 
orthogonal  and  are  passed  through  an  electro-optic  modulator  controlled  by 
an  applied  electric  voltage.^'*’*''**’’®*’’*  This  is  a  reverse  Fredkin  gate 
design.  When  voltage  is  applied  to  the  modulator,  the  polarizations  of  both 

arguments  rotate  by  90**  and  switching  occurs.  A  similar  method  has  been 
demonskaled  using  a  polarization  cobr-coded  liquid  crystal  light  valve.^®*^^ 

2)  Aoousto-optb  gate:  The  arguments  are  passed  through  an  acousto- 
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Optic  deflector  arKi  the  control  line  is  an  acoustic  signal  applied  to  the 
■Aio8.i64.ts6  jg  gigQ  g  reverse  Fredkin  gate  design.  Switching 
occurs  in  the  presence  of  the  acoustical  signal. 

3)  Photorefiective  gate:  Based  on  four  wa  /e  mixing,  the  photorefractive 
gate  is  an  all  optical  switch  where  the  arguments  are  signal  beams 
counterpropagated  through  a  photorefractive  material.  The  control  line  is  the 
pump  wid  it  is  also  counterpropagated  through  the  photorefractive 
material.^’’*'^^'^®‘*®^'^®®*^®^'^“  This  s,  again,  a  reverse  Fredkin  gate  design. 
When  the  pump  is  present,  switching  occurs  via  phase  conjugation. 

4)  Waveguide  switch:  The  arguments  are  symmetric  coupled  channel 
waveguides  in  an  electroK>pticaily  active  material.  Depending  upon  the 
extent  of  the  optical  activity  of  the  material,  the  control  channel  is  either  an 
eiectricaf  or  optical  signal  applied  across  the  section  of  the  two  argument 
channels  where  coupling  oocurs.®^'^“’'“’^“ 

Some  of  the  designs  for  waveguide  switches  include:  directional  coupler 

zero-gap  directional 

coupler  switch,  which  is  a  riKxlified  version  of  the  directional  coupler 

and  crossing  channel  switch  or  X-switch.^®’^®-*-®* 
a4^6j7.36.iHi86  dlrecllonal  coupler  and  X-switch  designs  all  operate  as 
Fredkin  gates.  The  channel  waveguide  zero-gap  symmetric  directional 
coupler  switch,  which  is  the  focus  of  this  dissertation,  is  illustrated  in  Rg. 


1 .3.1 .  Researchers  have  thoroughly  investigated  the  theory  behind  channel 
waveguide  directionai  coupler  switches  and  several  of  these  devices  have 
been  demonstrated  using  lithium  niobate  (LiNbOa),  indium  phosphide  (InP), 
GaAs  and  AIGaAs  as  the  waveguiding  medium.®’’®^^'®®^’'^®®’^’®®' 

107,1 12,1 14,1 1S.1 17,1 19.142.164,166 


Directional  Coupler  Switch 

direction«t  >  •ipler  switch  is  a  device  consisting  of  parallel  channel 
waveguides  separated  by  a  finite  distance.  The  coupling  between  the 
modes  of  the  parallel  waveguides  results  in  an  exchange  of  power  between 
guided  modes  of  adjacent  waveguides.^®  This  is  referred  to  as  directional 
coupling.  The  eigenmodes  of  the  coupling  region  consist  of  or>e  symmetric 
and  one  antisymmetric  inode.^®  "'^®'®*'®®’®^*^®’’®®^  Treatment  of  waveguide 
coupling  can  be  performed  by  coupled-mode  theory.^®*"’^®’^®®  Consider  Pii 
as  the  input  power  to  waveguide  1,  Piz  as  the  input  power  to  waveguide  2, 
POf  as  the  output  power  from  waveguide  1  and  P02  as  the  output  power 
from  waveguide  2.  Let  waveguide  1  propagate  through  the  symmetric  mode 
and  waveguide  2  propagate  through  the  antisymmetric  mode.  Assume  that 
the  waveguides  are  not  too  close  so  that  the  overlap  integral  of  the  mode 
functions  will  be  small.^®’^’"'^®'^®®  With  a  single  input  to  waveguide  1 ,  Pi^  =  1 


then  Pis  »  0.  The  output  power  from  waveguides  1  and  2  are  given 


Input  Channel  1  ^  Output  Channel  1 


Figure  1.3.1 

Schematic  of  a  Zero^Bap  Directional  Coupler 
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Po,(l)  =  PI,  -  P0,(D 

=  Pl,(1  -  (x*/(x*  +  5*))(1  -  cos*[(ic*  +  5*)'®(«)]))  1 .4.1 

and 

POad)  =  Pii((K*/(K*  +  5*))(sin"[(K*  +  5")^'*(l)l)),  1 .4.2 

where  2S  s  •  p2  is  the  difference  in  the  propagation  constants  between 

two  adjacent  uncoupled  waveguides  1  and  2.  k  is  the  coupling  constant 
between  the  two  adjacent  waveguides  and  I  is  the  interaction  length  over 
which  the  two  waveguides  are  coupled.  By  definition,  k  » jc/2L,  where  I  s 

L  is  the  length  required  for  complete  coupling  of  the  power  from  one 
waveguide  to  the 

For  symmetric  waveguides,  the  phase  veiccities  in  the  two  modes  are 
equal  (i.e.,  p,  s  ^  or  5  »  0).  For  complete  coupling,  Eqs.  1.4.1  and  1.4.2 
become’*"'i®'i”-io® 

Po,(L)  s  Pi,(i  -  (1  -  cos®[«/2]))  *  0  1 .4.3 

and 


P0i(L)  =  Pi,(sin*[K/2])  =  PI,. 


1.4.4 
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From  Eqs.  1.4.3  and  1.4.4,  one  can  see  that  the  input  power  Pi,  into 
waveguide  1  will  exit  waveguide  2  as  P02  at  a  distance  L  Using  Eqs.  1 .4.3 
and  1.4.4  complete  coupling  occurs  when  L  s  ic/2k. 

Applying  an  electric  fieid  to  either  waveguide  1  or  waveguide  2,  over  the 
distanoe  where  the  waveguides  are  coupled,  will  change  that  particular 

waveguide's  propagation  constant  Total  switching  occurs  when  5  s 
This  can  be  verified  by  the  substitution  of  5  =  v^it/2L  into  Eqs. 
1.4.1  arKi  1.4.2.  By  doing  this  we  get'6.^7.i».io2.io6 

Po,(L) «  Pi,(1  -  (1/4)(1  -  COS*[K]))  «  Pi, 
and 

Po^iL)  =  Pl,((1/4)(sin*[ic]))  =  0. 

Using  the  definitions  for  6  and  x  given  above,  one  gets  total  switching  to 
occur  when^*’"’^”'^®* 

0,  -  P2)  =  VSit/L  1.4.7 

If  we  apply  an  electric  field  (£)  over  waveguide  2,  for  instance,  this  will 
HfKiuoe  a  change  in  the  refractive  index  of  waveguide  2  by  an  amount 
proportional  to  E  Mathematicaliy  this  is  given 


1.4.5 

1.4.6 


AHs « {n^i2)rE, 
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1.4.8 


where  ris  is  the  refractive  index  of  waveguide  2  and  r  is  defined  as  the 
electro-optic  coefficient  of  the  waveguiding  niaterial.  To  apply  an  electric 
field  to  one  of  the  waveguides,  we  place  electrodes  on  the  top  and  bottom 
of  the  waveguide  over  the  coupling  length  (L)  and  apply  a  voltage  (\/)  to  the 
top  electrode.  The  electric  field  Is  then  related  to  the  applied  voltage  by  V 
s  Ed,  where  dis  the  distance  between  the  top  and  bottom  electrode.  The 
total  index  (ngO  of  the  waveguide  2  region  between  the  electrodes  can  be 
defined  by 

n2'*n2  +  An2  =  n2  +  {r\^l2)rVJd.  1 .4.9 

By  definition,  p  »  Znnfk,  where  A.  is  the  wavelength  of  the  source.^^ 

If  we  let  waveguide  1  and  waveguide  2  be  identical  (n^  s  ns  »  n)  and 

subetitute  Eq.  1.4.9  into  1.4.7,  along  with  the  definition  for  p,  this 
yieWs^^'i^-io® 

(2x/Ji)(n'  -  n)(n*/2)(rWd)  =  (2tdk)(n^/2)(rWc/)  =  1 .4.10 

We  can  now  find  the  voltage  required  to  completely  switch  the  coupling.  By 
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rewriting  Eq.  1.4.10,  we  get 

»/=  (i/3XcO/(n2*/l).  1.4.11 

The  operation  of  a  directional  coupler  switch  is  illustrated  in  Fig.  1 .4.1 . 


As  with  directional  couplers,  the  interaction  region  of  a  zero-gap 
directional  coupler  wU  support  two  guided  symmetric  rrKxies  and  two  guided 
antisymmetric  modes.®*^*^^^'*'”'*  The  differerrce  between  them,  however, 
results  from  the  fact  that  we  have  zero  spacing  between  the  two  symmetric 
waveguidea  in  a  zero-gap  directional  coupler.  These  modes  win  now  be 
propagating  in  a  single  ,  waveguide.  We  will  again  have  a  voltage 
independent  phase  difference  that  accumulates  over  the  interaction  length 
Pte  two  mode  section),  due  to  the  difference  in  the  phase  velocities  of  the 
two  orthogonal  modes,  and  a  phase  difference  over  the  two  mode  section, 
brought  about  by  the  voltage  induced  change  in  refractive  index  (electro- 
optic  (EO) 

The  electro-optic  dependent  phase  difference  will  be  denoted  by  (P,  - 
s  Afi/2,  where  ^  and  are  the  propagation  constants  of  the 
symmetric  and  antisymmetric  modes,  respectively.®^  The  value  of 
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required  for  complete  coupling  of  energy  from  the  top  half  of  the  two  mode 


waveguide  to  the  bottom  is  given 


Ap  s  a/L. 


1.5.1 


aP  is  varied  by  changing  the  index  of  refraction  of  the  waveguiding  material 
in  the  twoHTXxfe  section.  We  will  denote  the  refractive  index  change,  as  we 

did  in  the  fast  section,  by  An.  The  change  in  the  difference  between  the 
propagation  constants  of  the  two  modes  can  then  be  written  as^^ 


a(aP)  »  An(aAP^.  1.5.2 

Equating  Eqs.  1 .5.1  and  1 .5.2  yieids^^ 

An  *  (ic/L)(aAP^)*\  1.5.3 

It  was  shown  in  the  last  section  that  for  linear  EO  materials,  the  refractive 
index  changes  proportionaliy  to  the  applied  voltage.  With  this  relation,  it  can 

be  seen  from  Eq.  1.5.3  that  dAp/9n  is  inversely  proportional  to  the  applied 
voltage.  dAp/dn  is  aisc  strongly  dependent  on  the  waveguide  separation  s 

for  the  dhectionai  ooupler.^^  Papuchon  and  Ro/^  showed  that  as  s 


ill! 
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ctscrMsee,  d(Ap)^  increasee,  and  as  the  value  of  increases,  the 
voltage  reqtilrsdfdr  switching  decreases.  Therefore,  the  rnaximurn  value  of 

dA^i^  and,  thus,  the  mirymum  voltage  required  to  switch  the  coupling 
guidee,  occurs  when  s  »  0.  When  s  «  0.  we  have  a  zero-gap  directional 

coupler.’^  ExperimerSs  have  shown  that  the  value  dA^^  varies  by  an  order 
of  rnagriltude  between  the  ooriventional  directional  coupler  and  the  zero-gap 
directional  oouplar.^  This  verifiee  that  zero-gap  directional  couplers  require 
leas  voltage  to  switch  the  coupling  than  conversional  directional  couplers. 

Expressions  for  the  power  ooupHng  in  zero-gap  directional  couplers  are 
sinriilar  to  thoee  for  oonveritiorial  directional  couplers.  They  are  given  by^^ 


Po,(l)  *  Pi,ooa^(Ap/2  ^  40)1}  1 .5.4 

and 

POgCD  «  Pi,sin*[(A^  -f-  ^)ll,  1 .5.5 

k>  is  the  aocumuialed,  voltage  independent  phase  difference  in  the 
Mfodties  of  the  symmetric  and  antisymmetric  modes.  With  no 

applied  ({.e..  s  O)  coupling  occurs  when  Oo  =  n/2L  When 
voiage  Is  applied  so  that  a^  »  k/L,  then  complete  switching  occurs. 
Equaling  Ap  with  Eq.  1.4.9  yields 


(2ic/X)(n’/2)(fV7d) « idL 
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1.5.6 


The  voltage  required  for  complete  switching  is  then  given  by 

(Xfl)/(n»/t),  1.5.7 

where  n  is  the  refractive  index  of  the  interaction  region.  By  comparing  Eq. 
1 .5.7  with  Eq.  1 .4.1 1 ,  one  can  see  that  the  voltage  required  to  completeiy 
switch  the  optical  input  from  one  channel  to  the  other  is  less  for  the  zero-gap 
directional  coupler.  The  zero-gap  allows  us  to  approximately  treat  the  center 
interaction  region  as  one  waveguide.  It  wW  be  shown  in  Chapter  IV  that  Eq. 
1.5.7  is  the  same  expression  derived  for  an  electro-optic  modulator. 
Operation  of  the  zero-gap  directional  coupler  switch  is  illustrated  in  Fig. 
1.5.1.  Another  advantage  of  zero-gap  directional  couplers,  in  addition  to 
operating  with  less  voltage,  is  that  the  interaction  length  is  much  shorter  than 
the  interaction  length  of  conventional  directional  couplers.*^*”  This  results 
in  the  reduction  of  real  estate  required.  Interaction  length  deperxience  on 
waveguide  separation  will  be  treated  in  Chapter  III. 
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Wavefront  Tin 


An  additional  parameter  that  needs  to  be  considered  for  both 
conventional  and  zero-gap  directional  couplers  is  wavefront  tilt  and  its  affect 

on  waveguide  coupling  in  the  branching  section.  The  larger  the  tilt  angle  a 
of  the  branchirig  section,  the  less  the  waveguide  coupling^*  (see  Fig.  1 .3.1 ). 

Tilt  introduces  phase  differences  that  decrease  waveguide  interaction  in  the 
branching  section  of  waveguide  coupling.  The  maximum  allowable  tilt  angle 
can  be  obtained  by  forcing  the  branching  region  of  waveguide  coupling  to  be 
within  a  single  lobe  of  the  standing  wave  generated  by  the  two  waves 
propagating  along  the  tilted  waveguide  portions.  This  gives  the  limit'* 

8ln(a^)  <  "kJAnw,  1 .6.1 

where  w  is  the  width  of  the  waveguides. 

Electrode  Separation 

The  separation  gap  g  between  the  electrodes  over  the  coupling  regions 
of  waveguides  1  and  2,  even  though  not  being  part  of  the  equations  above, 
playe  a  major  role  in  switch  operation  (see  Fig.  1.3.1).  It  has  been  found 
thsd  the  wider  the  gap,  the  lower  the  electric  field  strength  and  induced 
relmctive  index  change.  The  narrower  the  gap,  the  higher  the  electric  field 
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Strength  and  induced  refractive  irKiex  change.”  it  has  aiso  been  presented 
that  the  electrode  separation  gap  depends  on  the  waveguide  width.  H 
should  be  much  less  than  w,  so  that  only  the  propagation  constant  of  the 
fundamental  mode  is  stror^gly  modified.”'”  There  are  limitations,  however, 
on  how  narrow  or  wide  g  can  be.  This  is  due  to  the  physical  reason  that 
variations  in  the  electrode  gap  can  only  be  partially  compensated  by  applied 
voltage.”  Complete  compensation  is  not  possible.  Overlap  of  the  applied 
electric  field  with  the  optical  electric  field  within  the  waveguide,  is  a  factor 
that  affects  the  behavior  of  the  index  change.”  it  is  brought  about  by  the 
electrode  separation  gap  versus  the  applied  voltage.  This  leads  to  the  need 
to  determine  an  optirmjm  g  versus  waveguide  width. 

Now  that  the  Fredkin  gate  and  zero-gap  directional  coupler  have  been 
introduced,  the  theory  and  background  refererices  necessary  for  completing 
this  research  will  be  reviewed  in  Chapters  li.  III  and  IV.  Chapters  ii  and  IV 
of  the  dissertation  will  describe  the  optical  and  electro-optical  properties  of 
AIQaAs.  The  characteristics  of  conventional  and  zero-gap  directional 
couplers  will  be  presented  in  Chapter  III.  Chapter  V  will  detail  the 
experimental  procedures,  observations  and  recorded  data  of  this  work.  This 
includBs  material  growth,  device  fabrication  and  operation.  Results  are  then 
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piBsentod  in  Chapter  VI  and  conclusions,  along  with  recommendations, 
given  in  Chapter  Vli. 


m 


CHAPTER  II 

OPTICAL  PROPERTIES  OF 
ALUMINUM  GALLIUM  ARSENIDE 


By  convention,  the  nomenclature  for  aluminum  gaiiium  arsenide  is 
generally  given  as  Ai^^Qa^^,  where  x  is  the  fraction  of  aluminum  atoms 
present  in  the  material  and  l-x  is  the  fraction  of  gaiiium  atoms  present  for 
each  atom  of  arsenic  in  the  material.  Research  has  demonstrated  that  the 
optical  properties  of  Ai,,Gaf.xAs  have  a  strong  dependence  on  x  and  the 
optical  wavelength  (X).  In  this  chapter,  we  will  quantitatively  investigate  the 

refractive  index  (n)  as  a  function  of  t>oth  x  and  X,  and  then  determine  these 
parameters  for  the  specific  AlxGa^.,^  structure  and  wavelength  to  be  used 

for  the  experiment. 

Refractive  index 

A  literature  search  was  conducted  for  references  investigating  the 
refractive  index  of  AlxGa,.xAs.  Three  of  the  sources  found^*^'^  were  ones 
in  which  the  irtdex  had  actually  been  measured  as  a  function  of  both  x  and 
X  and  that  these  three  sources  are  the  most  widely  referenced  among 

38 
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At^Ga,^  integrated  optics  researchers.®**^"'”'^®®^®^  The  iiterature 
search  also  uncovered  a  few  sources  that  fit  these  experimental  data  to 
analytical  models  to  interpolate  the  refractive  index  of  various  AlxQa^.xAs 
compositions  for  which  the  index  had  not  been  measured.  Sections  2.1  -  2.3 
will  outline  all  of  the  pertinent  sources.  These  sources  will  be  used,  along 
with  information  gathered  through  conversations  with  several  researchers  in 
the  technical  area,  listed  in  the  bibliography,  to  determine  the  refractive  index 

for  the  specific  x  and  X  used  for  this  experiment. 

Before  going  onto  the  modelling  sections,  the  wavelength  and  the  Al 

concentration  of  the  waveguide  devices  to  be  fabricated  for  this  experiment 

need  to  be  chosen.  It  is  desirable  to  have  an  x  and  X  that  will  yield  as  low 
a  propagation  loss  as  possible.  Selection  of  X  was  based  on  the  type  of 
laser  diode  used  in  the  experiment.  Since  the  absorption  is  inversely 

proportionai  to  X,®®  the  higher  the  X,  the  lower  the  loss. 

Among  the  oommerdaliy  sold  laser  diodes  with  as  high  a  X  as  possible, 
Al^Qai^  laser  diodes  in  the  760  •  850  nm  range  are  readily  available  as 

well  as  inexpensive.  Ai^Ga^^  laser  diodes  with  X  >  850  nm  are  usually 

more  expensive  and  not  as  readily  available.  For  X  >  900  nm,  the  material 
used  for  the  laser  structure  included  either  indium  or  phosphorus,  depending 

on  the  desired  X. 

Relevant  and  complementary  United  States  Air  Force  interests  include: 
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(a)  Three  Terminal  Opto-eiectronic  Arrays  (TTOA)  for  which  a  laser, 
photodetector,  and  amplifier  structure  are  all  being  fabricated  into  a  single 
pixel  area  in  material;^  (b)  Integrated  Optical  Organic  Devices 

(lOOO),  for  which  polymer  based  optical  waveguides  are  being  developed 
for  operation  within  the  700  *  900  nm  wavelength  range;^  and  (c)  Optical 
Interconnect  Technology  (OIT),  for  which  integration  of  laser  source,  optical 
organic  waveguide,  photodetector  and  electronics  are  ail  being  integrated 
together  to  realize  a  complete  multichip  module  optical  interconnect 
operating  in  the  700  •  900  nm  wavelength  range.^’^^'  These  research 
programs  dictated  the  use  of  an  AiKGa,.xAs  laser  diode  with  an  operating  X 
somewhere  between  700  and  900  nm. 

After  considering  compatibility  with  the  TTOA,  lOOD  and  OIT  programs 
along  with  price,  availability,  and  monolithic  Integration  of  the  laser, 
waveguide  and  detector,  a  Sharp  Model  LT015MD0  40  mW  index  guided 

laser  diode,  operating  eX  X  -  827  nm,  was  selected.  A  spectrometric 
measurement  of  the  laser  diode  was  taken  with  an  input  current  of  81 .0  mA 

and  maintained  at  a  temperature  of  25.0^0  to  verify  X. 

With  X  chosen,  the  Ai  concentration  to  be  used  for  the  core  of  the 
waveguide,  that  would  yield  as  low  a  propagation  loss  as  possible,  needed 

to  be  determined.  At  a  given  Xthe  absorption  is  inversely  proportional  to  the 
fraction  of  Ai  atoms.^  The  higher  the  x  value,  the  lower  the  loss,  in 
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addition,  the  required  waveguide  thickness  for  a  given  X  is  inverseiy 
proportionaitox.  These  criteria  would  suggest  a  large  Al  concentration.  For 

X  >  0.350,  however,  the  bandgap  level  becomes  indirect,  resulting  in 
diminished  optical  activity,  due  to  the  fact  that  both  a  photon  and  phonon  are 
reqiNred  for  an  indirect  transition.^  This  dictates  that  x  ^  0.350.  it  has  also 
been  reported^  that  at  wavelengths  Hi  the  800  >  850  nm  range,  the  modal 
attenuations  increase  markedly  for  a  waveguide  core  layer  when  Al 
conoentrations  are  above  0.310.  Accordingly,  an  x  value  of  0.300  was 
chosen  for  the  waveguide  core. 

Based  on  the  findings  given  above,  the  Al  concentration  is  constrained 
to  the  range  0.200  ^  x  ^  0.400  when  0.7  pm  ^  A,  1 .0  |im  for  the  AI^Qa^^ 
refractive  index  modeis.  In  addition,  the  models  were  fit  to  the  measured 

indices  by  optimizing  for  an  Al  concentration  x  s  0.300  and  X  827  nm.  This 
provided  better  agreement  with  the  measured  refractive  index  in  this  limited 

range  of  x  and  X  than  when  the  entire,  unlimited,  ranges  of  x  and  X  were 
considered. 

Section  2.1 

In  1972,  Boyd^  fit  to  Sellmeir  equation  experimental  data  published  by 
llegems  and  Pearson^: 


40 

(n/  «  A  +  B/(X*  -  C)  -  DX*.  2.1 .1 

where  » the  index  of  refraction,  X  is  the  wavelength  given  in  and  A, 
B,  C,  and  D  are  the  Seiimeir  coefficients,  of  which  A,  C,  and  D  are  functions 
of  X.  Ilegems  and  Pearson  used  reflectivity  measurements  to  determine  the 
refractive  index  of  AI^Qai^  for  x  values  of  0.080, 0.210, 0.380, 0.530,  arxl 
0.820,  respectively.  The  accuracy  of  x  was  known  to  ±  0.020,^  which  is 
within  the  toieranoe  of  the  various  electronic  material  characterization 
methods.”^  Numerical  values  for  A,  B,  C,  and  D  can  be  found  for  various 
Al  concentrations  by  using  the  equations  generated  by  Boyd  as  iisted  in 
Table  2.1.1. 


Table  2.1.1* 

Seiimeir  Coefficients  for  Ai^Gai^ 
ABC  D 


10.906  -  2.92X 

0.97501 

(0.52886  •  0.735x)^ 

[0.002467 

x^O.36 

•(1.41X  +  1)] 

(0.30386  -  0.105X)* 

x^O.36 

Several  of  the  assumptions  used  by  Boyd  were:  1)  Linear  variation  of  A  with 
x;  2)  B  independent  of  x;  3)  Quadratic  variation  of  C  with  x  to  incorporate 
ttie  shift  in  bandgap  as  a  function  of  x,  which  varies  lineariy  with  a  different 
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Slope  in  the  regions  of  direct  and  indirect  band  transition;  4)  D  varies 
Unearty  with  X. 

Thus,  Eq.  2.1 .1  allows  one  to  determine  the  refractive  index  (nj  as  a 

function  of  both  x  and  X.  Fig.  2.1 .1  shows  the  plots  of  the  Sellmeir  equation 
for  two  X  values,  0.200  and  0.400,  as  a  function  of  X,  between  0.7  and  1 .0 
pm.  It  is  suggested  by  Hunsperger,*  that  Eq.  2.1 .1 ,  along  with  the  values 

for  A,  6,  C,  and  0,  can  also  be  used  to  model  the  refractive  indices 
measured  by  Casey,  Sell  and  Panish.^^  When  the  values  for  x,  given  in  that 
latter  work,  were  substituted  into  Eq.  2.1 .1 ,  however,  the  modelled  index  did 
not  fit  the  experimentally  measured  index  very  well.  For  the  indices 
measured  by  Casey  et  ai.,  a  second  model  was  found  that  fit  these 
measurements  much  better.  This  model  is  discussed  in  the  next  section. 

Section  2.2 

In  1985,  Adachl*  modelled  the  refractive  index  data  of  Ai^GSi^  obtained 
by  Casey,  Sen  and  Panish.^^  These  measured  indices  were  determined  by 
refiecfivity  measurements  for  x  values  of  0.070, 0.100, 0.150, 0.200, 0.240, 
0290  and  0.380  at  photon  energies  ranging  frcxn  1 .2  to  1 .8  electron  volts 

(eV).  The  accuracy  with  which  x  was  known  in  that  experiment  was  ± 
0.020.^  Like  the  previous  index  measurements,  this  accuracy  is  within  the 

tolerance  of  the  various  electronic  materials  characterization  methods.^^' 
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AdacN  usee  the  relation^ 


(n2)*(oa)  *  8,(0#)  +  /EjN,  2^.1 

to  model  the  refractive  index  where  (d  is  the  angular  frequency,  ni(<o)  is  the 
frequency  dependent  refractive  irKtex  and  e^Cco)  and  eaCo)  are  the  real  and 
imaginary  parts  of  the  frequency  dependent  dielectric  constant,  respectively. 
Adachfs  relation,  however,  is  valid  for  only  those  energies  below  the  direct 
band  edge.  At  or  near  the  lowest  direct  gap  (EJ.  the  imaginaiy  part  of  the 
dieleciric  constant  (^(o>))  may  be  taken  as  zero  when  Eq.  2.2.1  becomes 

(na)*(o>)  -  8i(<o)i  2.2.2 

where  the  real  part  of  the  dielectric  constant  (e^Co)))  for  zincblende  materials 
like  AixQai^,  at  energies  below  the  E^  direct  gap  can  be  expressed  using 
the  well  known  Kramers-Kronig  relation^ 

e,(«>)  =  Mf(x)  +  1/2[Ey(E.  +  *J]”/6u)}  +  2^.3 


where 


/(X)  -  (X)'*P  -  (1  +  X)“  ■  (1  -  X>'®1. 


22.4 


/(XJ  «  (XJ^  •  (1  +  XJ’"  -  (1  -  7(J'% 


2J2.5 


X  *  22.6 

and 

X»*W(Eo  +  iO.  25.7 

In  Eq.  25.3  -  25.7,  and  are  the  lowest  direct  gap  and  direct  gap 

apming  energlea,  or  critical  point  energies,  in  eV,  Vd  is  the  photon  energy 
in  eV  where  Is  Plancks  constant  and  Ao  and  Bo  are  constants  used  to  fit 

the  model  to  the  measured  indices.  Using  the  measured  values  of  the 
cddcal  point  energies  taken  from  Berolo  and  Woolley,"  Eo  and  Eo  +  can 
be  lilted  to  the  parabolic  relation" 

E  «  a  bx  Q^.  25.8 

The  constants  a,  b  and  c  were  determined  for  Eo  by  Lee>  Juravel,  Woolley 
and  SpringThorpe,*^  by  fitting  Eq.  25.8  to  the  measured  values.  Using 
these  values,  the  equation  for  Eo  can  be  given  as^^ 


Eo  «  1.425  +  1.155X  +  0.37)^. 


25.9 
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AdacN  later  detenmlned  the  constants  a.  b  and  c  for  in  the  same 
manner.  Using  Adachi*s  values,  the  equation  for  Eo  ***  is  given  as^ 

Eo  +  Ao  *  1.765  1.115X  +  0.373^.  2.2.10 

Eq  and  Eo  ^  ^  plotted  as  a  function  of  x  in  Fig.  2.2.1 .  The  fit  of  the 

curves  generated  from  Eqs.  2.2.9  and  2.2.10  with  Berolo  et  al.  is  quite  good. 
By  substituting  Eqs.  2.2.9  and  2.2.10  into  Eqs.  2.2.2  •  2.2.7  and  then  fitting 
Eq.  2.2.2  to  the  indices  measured  by  Casey  et  al.,  A,  and  B^  can  be  found. 
Dolrig  this  yields 

Ao  »  6.34  +  18.39X  +  0.78)^  2.2.1 1 

and 

Bo  =  9.38  -  10.61X  +  0.78X*.  2.2.12 

These  vaiuee  for  Ai,  arxl  Bo  are  not  the  ones  used  by  Adachi^,  however, 
they  proved  to  give  a  better  fit  to  the  measured  refractive  index  of  Casey  et 
ai.  than  Adachfs  values.  Figs.  2.2.2  •  22.5  show  plots  of  the  refractive 
index  (ns),  including  the  ±  0.020  Ai  concentration  accuracy,  for  the  x  values 

0.200, 0.240, 0.290  and  0.380  as  a  function  of  )i(D  (between  1 2  and  1 .8  eV), 
which  corresponds  to,  approximately,  1.0  pm  ^  X  ^  0.7  pm.  Figs.  2.2.2  - 


V)0 


Figure  2J 

Refractive  Index  of 
Photon  Energy  for  x  = 
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2J2.S  also  include  the  refractive  indices  measured  by  Casey  et  al.,  which 

were  substituted  Into  the  model  at  the  specific  x  and  X  for  which  the  index 
was  actually  measured.  These  substituted  values  appear  as  small  spikes 

in  an  otherwise  smooth  curve,  giving  an  indication  of  the  accuracy  of  the 
model.  The  refractive  indices,  using  the  above  nKXiel  and  curve  fitting 
constants,  are  within  the  known  accuracy  of  the  fraction  of  At  atoms.  It 
provides  a  very  good  fit  around  the  chosen  x  =  0.300  waveguide  core  and 
the  827  nm  wavelength,  =  1 .5  eV. 

Section  2.3 

in  1986,  a  new  set  of  refractive  index  measurements  were  made  by 
Aspnes,  Kelso,  Logan  and  Bhat.”  The  indices  were  measured  using  a 
rotating  analyzer  eliipsometric  technique  for  x  values  of  0.099, 0.1 98, 0.31 5, 
0.419,  0.491, 0.590,  0.700  and  0.804  at  photon  energies  ranging  from  1.5 
to  6  eV.  The  accuracy  with  which  x  was  known  was  ±  0.020.^  This 
accuracy  is  within  the  tolerance  of  the  various  electronic  materials 
characterization  methods.^^^  The  Adachi  model  presented  in  the  previous 
section  did  not  provide  a  dose  enough  fit  to  be  used  to  model  these 
measurements  taken  by  Aspnes  et  al.  However,  in  1 989,  AdachP  improved 
his  model  to  better  fit  these  experimental  data.  This  new  model  includes 
contributions  due  to  Eo  Eg  +  Aq  as  well  as  the  higher  critical  point 
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energies  E^,  +  a^,  E^,  and  Eo'  +  a^',  including  nondispersive 

contributions. 

We  will  determine  all  the  criticai  point  energies  and  nondispersive  terms 
before  presenting  the  model,  since  these  energies  will  consequently  be  used 
to  determine  the  refractive  irKiex.  The  critical  point  energies  E^  and  were 
measured  by  Aspnes  et  al.^  as  a  function  of  x.  These  measured  critical 
point  energies  and  E^,  as  well  as  E^  +  a^  and  E2,  are  given  by  Adachi.^ 
Equations  for  the  criticai  point  energies,  fit  to  the  values  found  in  Adachi,^ 
are  then  given  as 

Eo  =  1.42  +  1.195X  +  0.315)^,  2.3.1 

Eo  +  Ao  =  1 .77  +  1 .079x  +  0.420)^,  2.3.2 

E,  =  2.96  +  0.456X  +  0.136x*  +  0.213X®  +  0.067)(*  2.3.3 

and 

Eg  =  4.70.  2.3.4 

The  critical  point  energy  E^'  and  the  direct  gap  splitting  energies  a^  and  a^' 
are  not  listed  by  Adachi.  The  only  measurements  found  for  E^  +  a^.  Eq'  and 
Eo'  +  Ao'  were  taken  by  Berdo  et  al.®®  Therefore,  E^'  and  E^'  +  a^'  were 
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fitted  directiy  to  Berolo  et  al.  giving 

Eo' s  4.49  0.095X  2.3.5 

and 

Eo'  +  Ao'  =  4.65  +  0.09QX.  2.3.6 

The  E^  +  critical  point  energy  was  determined  in  a  different  way.  Only 
Berolo  et  al.  and  Aspnes  et  al.  measured  E..  The  critical  point  energies  E^ 
and  E,  +  a^,  measured  by  Berolo  et  al.,  appear  to  have  the  same  basic 
shape,  separated  only  by  a  constant  eV.  Using  the  critical  point  energy 
measurements  for  E^  and  E^  +  a^  (from  Berolo  et  al.)  and  the  critical  point 
energy  measurement  for  E^  (from  Aspnes  et  al.),  an  approximation  for  the 
criticai  point  energy  E^  -i-  a^,  to  be  used  for  this  new  model,  was  found  by 
taking  the  difference  between  E^  and  E,  +  a^  from  Berolo  et  al.  and  then 
shifting  E^  from  Aspnes  et  al.  by  the  same  eV.  E^  +  a^  was  then 
approximated  by 

El  +  A^  =  3.19  +  0.456X  +  0.136x*  +  0.210X®  +  0.048/.  2.3.7 


The  critical  point  energies  given  by  Eqs.  2.3.1  -  2.3.7  are  plotted  as  a 
function  of  x  in  Fig.  2.3.1. 


Taken  from  Section  22,  the  refractive  index  is  found  here  by  using  the 
relation: 

(HafN  -  ei((D),  2.3.8 

where  n,,  the  refractive  index,  and  the  real  part  of  the  dielectric  constant 
(et(<o))  now  include  contributions  from  the  higher  transition  energies.  The 
value  of  the  real  part  of  the  dielectric  constant  due  to  the  and 
transitions^  is  given  by  the  Kramers-Kronig  relation^ 


Uo)  -  A(EJ”{/(XJ  +  1/2(E./(E,  +  Or/(xJ}.  2-3.9 

vntiere 

f(xJ  =  (Xo)^[2  •  (1  +  Xo)^®  -  (1  -  Xo)^®H(1  -  Xo)].  2.3.10 

f(xJ  =  (XJ®[2  -  (1  +  -  (1  -  XJ^®H(1  -  XJ].  2.3.1 1 

2.3.12 


2.3.13 


1  for  y 0 
0  for  y  <  0. 
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2.3.14 


I 


In  Eqs.  2.3.9  •  2.3.14,  and  E,,  -i-  are  given  in  eV,  imd  is  the  photon 

energy  in  eV  where  Ik  is  Plancks  constant  and  A  is  a  constant  used  to  fit  the 
model  to  the  experimental  data. 

The  value  of  the  real  part  of  the  dielectric  constant  due  to  the  E^ 
transition*^  and  E^  -i-  transition  (Eq.  2.3.7)  is  given  by** 

e„(ffl)  -  -BlOcyintl  -  (x,fl  -  B2(x.,)-*ln[1  •  (x.,)*).  8.3.15 

where 

and 

X„  a  W(E,  +  Ai). 


2.3.16 

2.3.17 


In  Eqs.  2.3.15  -  2.3.17,  E^  and  E^  -i-  a^  are  given  in  eV,  and  B1  and  B2  are 
constants  used  to  fit  the  model  to  the  experimental  data. 

The  value  of  the  real  part  of  the  dielectric  constant  due  to  the  E2 
transition**  is  given  by** 


£«(«>)  =  cn  -  (X2)7ai  -  (x*)*f  +  (x2)Y). 


2.3.18 


where 

Xa^licyEs.  2.3.19 

In  Eqe.  2.3.18  -  2.3.19,  Es  is  given  in  eV,  and  C  and  y  are  constants  used 
to  fit  the  model  to  the  experimental  data. 

The  value  of  the  real  part  of  the  dielectric  constant  due  to  the  Eo'  and 
Eo'  +  Aq'  transitions"  was  found  by  fitting  the  model  to  the  index 
measurements.  This  gives 

e,o'(®)  =  D(Xo'  +  TLol,  2.320 

where 

Xo'  *  WE,'  2.3.21 

arKi 

3U'  =  W(E/  +  A,!.  2.3.22 


In  Eqs.  2.3.20  <  2.3.22,  E,'  and  +  a,'  are  given  in  eV,  Ixo  is  the  photon 
energy  in  eV  and  D  is  a  constant  used  to  fit  the  model  to  the  experimental 
data.  The  value  of  the  nondispersive  elements  of  the  dielectric  constant  due 
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to  the  indiiect  gap  transitions  is  given  by  a  constant  which  is  used  to  fit 
the  model  to  the  experimental  data.^  The  real  part  of  the  total  dielectric 
constant  due  to  all  of  the  transition  energies  is  the  sum  of  the  dielectric 
constants  for  each  of  the  transition  energies  or 

EtCca)  a  £^,(<0)  +  +  £«(«>)  +  +  £i..  2.3^3 

The  strength  parameters  A,  C,  y,  and  were  determined  as  a  function 
of  Al  concentration  by  fitting  equations  to  the  curves  presented  by  Adachi.^ 
They  are  given  as  follows: 

A  a  2.9  15.9X  +  6.5)^  +  0.1x®  +  0.1X*  +  3x® 

+  13x*  +  32/,  2.3.24 


C  a  [2.31  -  0.001X  -  0.001X*  -  0.001/  -  0.001/ 

-  0.5/  - 1 .2/  +  2.5x^l,  2.3.25 

y  a  0.137  -  0.005X  -  0.0001/  -  0.0001/  -  0.02/  -  0.16/  2.3.26 


and 


ei.a1.25-1.98x-0.4/. 


2.327 


The  strength  parameters  B1  and  B2  were  estimated  by  using  the 
expressions" 

B1  «  44[E,  +  a^]/[8o(Ei)1  2.3^8 

and 

B2  «  44(E,  +  2A^yiao(E,  +  a^)*],  2.3^ 

a^iere  ao  is  the  lattice  constant  in  A  given  by" 

a^  s  5.6533  ■¥  0.0078X.  2.3.30 

The  strength  parameter  D  was  estimated  as  a  function  of  x  by  substituting 
Eq.  2.3.20  into  Eq.  2.3.23.  varying  D  until  Eq.  2.3.23  fit  the  experimental 
measurements  from  Aspnes  et  al."  for  both  x  and  photon  energy.  This 
yields 

D  »  2.535  -  1.075X  -  0.65x*  -  0.7x®  -  0.8/  -  2.8x®.  2.3.31 


Substitution  of  Eq.  2.3.23  into  2.3.8  gives  the  refractive  index  for  the 

given  x  values  as  a  function  of  tMio.  Figs.  2.3.2  -  2.3.4  show  the  plot  of  the 
index  (ru)  as  a  function  of  Wo,  between  1.2  and  1.8  eV,  for  the  x  values 
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0.186,  0.315  and  0.419,  respectively,  including  the  ±  0.020  accuracy  of  x. 
Figs.  2.3.2  •  2.2.4  also  Include  the  refractive  indices  measured  by  Aspnes 

et  al.,  vvhich  were  substituted  into  the  model  at  the  specific  x  and  X,  for  which 
the  iTKlex  was  actually  measured.  These  substituted  values  appear  as  small 

spikes  in  an  otherwise  snrK)Oth  curve,  giving  an  indication  of  the  accuracy  of 
the  model.  As  can  be  seen,  the  refractive  indices,  using  the  above  nruxiel 
arKi  curve  fitting  constants,  are  well  within  the  known  accuracy  of  x.  The 
spikes,  representing  the  measured  values,  are  nearly  undetectable.  This 
model  provides  an  excellent  fit,  not  only  around  the  chosen  x  =  0.300 
waveguide  core  and  827  nm  wavelength,  but  also  for  the  range  of  x  from 
0.200  •  0.400  and  the  range  of  X  from  0.700  •  1.0  pm. 


Refractive  Index  of  the  Waveguide  Core 
Using  the  models  from  Sections  2.1  •  2.3,  the  refractive  index  of 
A^Qat.KAs  can  be  determined  at  fractions  of  Al  atoms  for  which  the  index 
was  not  measured.  As  stated  above,  the  target  Ai  concentration  chosen  for 
the  waveguide  core  in  this  experiment  was  x  s  0.300.  Fig.  2.4.1  is  a  plot  of 
the  Indices  versus  wavelength  for  x  =  0.300  using  the  models  outlined  in 
Sections  2.1  >  2.3.  As  can  be  seen  in  Fig.  2.4.1 ,  there  is  a  significant 
dtfferenoe  between  the  three  curves.  The  curves  generated  by  Adachi’s 
models  (ns  and  ng)  based  on  Casey  et  al.  and  Aspnes  et  al.,  respectively, 
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Rdfractive 
Photon  Energy 


Figure  2.3.4 

Refractive  IrKiex  of  AlxGai.xAs  versus 
Photon  Energy  for  x  =  0.419  ±  0.020 
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are  offset  ah  from  one  another,  but  they  possess  the  same  basic  shape.  The 
slope  of  03,  above  X  =  0.9  pm,  appears  to  be  steeper  than  that  of  r^,. 

Modelling  the  refractive  index  for  x  =  0.300  along  with  the  ±  0.020  Al 
concentration  accuracy,  incorporating  Adachi’s  models,  shows  that  the 

difference  between  the  two  index  models  is  barely  within  this  ±  0.020  Al 
concentration  accuracy.  The  minimum  index  value  (x  =  0.320)  for  Adachi^, 

based  on  Casey  et  al.,  just  overlaps  the  maximum  index  value  (x  =  0.280) 
for  Adachi,^  based  on  Aspnes  et  al.  (see  Fig.  2.4.2).  In  Rg.  2.4.1 ,  the  curve 
for  n|,  generated  by  Boyd’s  Selimeir  equation,^  based  on  llegems  et  al.,  lies 
between  the  other  two  curves  with  a  slope  that  follows  the  other  two  curves, 
as  long  as  X  >  780  nm.  Below  780  nm,  the  slope  of  n,  falls  off  considerably, 
no  longer  following  the  other  two.  The  correct  curve  is  virtualiy  impossible 
to  determine  at  this  time  and,  as  mentioned  above,  all  three  curves  are  used 
almost  equally  by  researchers.  Based  on  Figs.  2.4.1  and  2.4.2,  along  with 
the  wavelength  of  the  laser  diode  used,  and  the  fact  that  the  indices 

measured  by  Aspnes  et  al.  were  not  recorded  above  X  >  827  nm,  it  was 
decided  to  further  restrict  the  wavelength  range  to  0.8  pm  ^  X  ^  0.9  pm  for 
determination  of  the  refractive  index  of  AlxGaf.xAs.  For  this  range  of  X  and 

the  0.200  ^  X  ^  0.400  given  earlier,  the  three  individual  models  follow  their 
measured  values  qufte  well,  offset  from  each  other  only  by  an  amount  An 
which  is  within  the  accuracy  of  the  known  Al  concentration.  Taking  all  three 

models  into  account,  it  was  decided  that  the  refractive  index  used  in  this 


dissertation  would  be  the  average  of  the  three  modelled  indices,  designated 
here  as  n^.  The  upper  limit  of  accuracy  was  chosen  to  be  ris  at  x  s  0.280 
(Adachi^,  since  it  gave  the  maximum  index  of  all  the  curves,  and  the  lower 
limit  of  accuracy  was  chosen  to  be  ng  at  x  =  0.320  (Adachi^,  since  it  gave 
the  minimum  index.  Fig.  2.4.3  shows  the  plot  of  n^  of  Alo^Gao.7oAs  versus 
wavelength,  including  the  upper  and  tower  limits  of  accuracy,  ns  at  0.280  and 
ng  at  0.320.  The  n^  value  of  the  Ai0.30Gao.70As  composition  will  be  given  the 
designation  noom>  ^  refractive  index  of  the  core  layer  of  the  waveguide. 

Calculating  the  average  index  of  refraction,  at  X  =  827  nm,  for  x  ~  0.300, 
yields  the  refractive  index  of  the  core  that  will  be  used  for  this  experiment. 

This  is  given  as 

Ooo»  =  =  3.419.  2.4.1 

The  "plus”  accuracy,  as  stated  above,  will  be  ris  while  the  "minus"  accuracy 
will  be  03.  Calculating  index  ns,  at  X  =  827  nm,  for  x  =  0.280,  and  index  ns 
at  X  s  827  nm,  for  x  =  0.320,  gives 


new,  =  3.419  ±  0.023. 


2.4.2 
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Refractive  Index  of  the  Waveguide  Clad 
Based  on  work  done  by  Honeywell  for  the  Air  Force,  the  optimum 
thickness  (tg)  for  the  waveguide  core  layer  of  an  AlojoGaoTo^s  composition 
was  determined  to  be  ^  s  1  pm.^  Assuming  a  single  mode  (m  =  0)  optical 
waveguide,  with  tg  s  i  pm,  X  -  827  nm,  noo^  given  in  the  previous  section 
and  the  relation* 


“  ^\iad  ^  rn*X  /[4(tg)  (riggiig  +  ^dad)]'  2.5.1 

the  vxfex  of  the  cladding  layer  can  be  found  as 

hetad  >  [(HeoJ*  -  X*/[4(g"]'®  2.5.2 

The  Al  coTKienMion  of  the  cladding  can  then  be  determir^  by  using  Eq. 

2.5.2  along  with  for  various  x  values  until  an  index  above  the  cutoff  is 

found.  Choosing  x  =  0.350,  at  X  s  827  nm,  one  gets  a  refractive  index  for 
the  dad  slightly  less  than  the  cutoff  index  (see  Fig.  2.5.1).  The  refractive 

index  for  a  cladding  composition  of  A^jgGaaesAs.  at  X  s  827  nm,  is 
determined  by  calculating  for  x  =  0.350.  This  gives 


-  3-389. 


2.5.3 
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A  TE/TM  mode  analysis  computer  program,  acquired  from  Honeywell^  arKl 
based  on  a  paper  by  Qhatak,  Thyagarajan  and  Shenoy,”  was  run  on  a 
waveguide  structure  with  a  core  composition  of  AIo^GSotoAs  and  a  dad 
composition  of  AI^^Gao,a6As  at  X  s  827  nm  to  assess  the  losses  for  the 
various  waveguide  modes.  The  program  predicted  a  0.00007  dB/cm  TE 
mode  loss  and  a  0.00009  dB/cm  TM  mode  loss  for  the  '0th  order"  mode, 
and  a  13.22  dB/cm  TE  mode  loss  and  a  16.54  dB/cm  TM  mode  loss  for  the 
'Ist  order*  mode.  This  suggests  that  even  though  noiad  does  not  quite  meet 
the  criteria  of  Eq.  2.5.2,  the  losses  of  the  higher  order  modes,  predicted  by 
the  TE/TM  mode  analysis,  are  sufficient  to  ensure  single  mode  propagation 
at  X  s  827  nm. 

The  moiecular  beam  epitaxy  (MBE)  growth  facility,  located  in  the  Solid 
State  Electronics  Directorate  at  Wright-Patterson  Air  Force  Base,  and  the 
metaForganic  chemicai  vapor  deposition  (MCX^VD)  growth  facility,  located 
at  Honeywell,  should  be  able  to  control  the  Ai  concentration,  over  a  range 

of  0.300  -  0.350,  to  within  ±  0.005  of  the  target  composition  for  both  the  first 
layer  and  the  difference  in  x  between  successive  layers.^’^"  This  will 

provide  exceilertt  repeatability.  Both  MBE  and  MOCVD  machines  will  be 
used  for  this  dissertation.  The  claim  of  the  accuracy  of  the  control  of  x  is  a 
valid  one  since  this  level  of  control  has  also  been  derTK)nstrated  on  the 
TTOA  program^  as  well  as  on  the  lll-V  Optical  Waveguide  Technology 


program.^  With  this  accuracy,  the  Al  concentration  of  the  waveguide  core 
could  be  as  low  as  X  s  0.295  and  the  Al  concentration  of  the  waveguide  dad 
could  be  as  high  as  0.355.  Fig.  2.5.2  is  a  plot  of  ncM>  for  x  -  0.355,  and 
ncuMt  calculated  from  Eq.  2.5.2,  using  for  x  s  0.295.  As  can  be  seen, 
nctad  is  slightly  less  than  n^^  at  Xs  827  nm,  indicating  the  possible  presence 
of  higher  order  modes.  The  TE/TM  mode  analysis  program,  run  on  this 

waveguide  structure  at  X  s  827  nm,  predicted  0.00001  dB/cm  TE/TM  mode 
losses  for  the  *0th  order”  mode  and  a  2.21  dB/cm  TE  mode  loss  and  2.86 

dB/cm  TM  mode  loss  for  the  *l8t  order”  mode.  Based  on  these  results, 
which  also  indicate  that  one  might  not  be  able  to  ensure  against  the 
presanoeof  higher  order  modes,  the  target  Al  concentration  for  the  dad  will 
be  sNghlly  less  that  35%.  Trying  to  grow  a  dad  with  x  <  0.350  while 
maintaining  x  «  0.300  for  the  core  is  difficijlt.^''"  Waveguide  strudures 
using  a  target  A^GaoToAs  core  composition  along  with  a  target 
AI^,^jsQ%«gAs  clad  composition  will  be  investigated  next. 


As  stated  in  the  previous  section,  the  accuracy  of  the  Al  concentration 

can  be  controlled  to  within  ±  0.005  for  this  experiment.  Using  this  accuracy 
when  determining  n^yg  for  x  s  0.300  at  X  s  827  nm  gives  a  core  index  of 


Figure  2.52 

Single  Mode  for  A^^Gao^TosAs  Core  and 
ncM  for  Alojs6(3ao44eAs  Clad  versus  Wavelength 
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rV«,*  3.419  ±  0.003.  2.5.5 

Using  the  same  accuracy  to  determine  for  x  s  0.350  at  X  s  827  nm 
gives  a  dad  irKlex  of 

0^  =  3.389  1  0.003.  2.5.6 

One  aiso  needs  to  take  into  account  the  difference  in  the  refractive  index 
measurements  taken  by  ilegems  et  ai.,  Casey  et  ai.  and  Aspnes  et  al.  This 
will  be  done  by  designating  three  configurations  of  waveguide  structures 
using  n^  along  with  the  *pius*  and  'minus*  accuracies  (02  and  rig)  for  the 
core  and  dad  layers.  Rrst,  the  average  refradh/e  indices  {n^  for  the  core 
and  dad  will  make  up  Waveguide  Strudure  1.  Next,  the  'plus'  accuracy 
indices  (ns)  for  the  core  and  dad,  modelled  by  Adachi^,  based  on  Casey  et 
ai.,  win  make  up  Waveguide  Strudure  2.  Last,  the  'minus'  accuracy  indices 
(Oa)  for  the  core  and  dad,  modelled  by  AdachP,  based  on  Aspnes  et  ai.,  will 
make  up  Waveguide  Strudure  3.  This  gives  the  possibility  of  three 
waveguide  strudures.  The  three  waveguide  strudures,  abng  with  the 
indices  of  rsfradion,  are  presented  bebw. 

Each  of  the  three  waveguide  strudures  will  be  used  for  determining  the 
interadion  length  of  the  optbai  diredbnai  coupbr  switches  that  will  model 
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the  optical  Fredkin  gates.  Comparing  the  three  structures  should  then 
provide  a  more  thorough  analysts  than  using  only  one  of  the  models  or  one 
set  of  measurements  for  the  refractive  index. 

WsYwyKte  gtructyr?  1 

(•’wv  Xoo«i  *  0-300  ±  0.005,  =  0.350  ±  0.005) 

core^  =  AlojooGao^TQoAs  clad^  =  Alo^Gao.eso^ 

n^o^i  =  3.419  ±  0.003  n^,^,  =  3.389  ±  0.003 

(Oa.  Xoo«e  *  0.315  ±  0.005.  =  0.365  ±  0.005) 

cof^  =  Alo^f5Qao,M5As  ciad2  =  Alo.38gGao,e3gAs 

IW  *  3.403  ±  0.003  n,^  =  3.375  ±  0.003 

MVJffljWg  SKrvglMr^  9 


(^*  ~  0uS85  ±  0.005,  XotKia  =  0.335  ±  0.005) 

corB3  =  AiQ^ggGaQ,7^gAs  clad3  —  Al^^g^gGan^^^p^^As 
rioons  ^  3.438  ±  0.003  -  3.404  ±  0.003 


CHAPTER  III 


DETERMINATION  OF  THE  INTERACTION  LENGTH 
OF  THE  ZERO-GAP  DIRECTIONAL  COUPLER 


In  this  chapter,  the  interaction  length  of  an  AI^Ga^.^  zero-gap 
directional  coupler  is  determined  by  substituting  the  refractive  index  of 
AI^Ga^^  into  several  dielectric  strip  waveguide  directional  coupler  nrKxieis, 
driving  the  separation  between  the  waveguides  to  zero. 

Directional  coupling  is  the  exchange  of  power  between  the  guided 
modes  in  adjacent  parallel  waveguides  separated  by  a  finite  distance.  Rg. 
3.0.1  is  a  schematic  of  a  conventional  directional  coupler  device.  The  iridex 
of  refraction  of  the  core  region  is  n  while  that  of  the  dad  regions  is  n(1  -  a), 

which  is  less  than  the  index  of  the  core  region  by  an  amount  a.  The 
dimensions  a,  b  and  c  are,  respectively,  the  width,  thickness  and  separation 

of  the  waveguides.  The  interaction  length,  or  distance  required  for  total 
exchange  of  power  from  one  guide  to  the  other  is  The  wave  is 
considered  to  propagate  along  the  z  axis. 

Extensive  treatment  of  optical  waveguide  and  directional  coupler  theory 

has  been  addressed  by  many  researchers  for  both  continuous  two- 
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dimensional  film  waveguides  and  dielectric  strip  waveguides.  Various 
methods,  such  as  dielectric  waveguide,  coupled  mode  theory,  effective 
index,  beam  propagation  method  and  finlte^lement  analysis  have  been 
16.17.2^61-M  Analysis  of  these  methods  shows  good  agreement  between 
all  of  the  different  models  provided  that  we  are  far  from  cutoff.  This  means 
that  the  normalized  frequency,  defined  as 

B  =  (W(/it)Cn  -  n(1  •  A)r^  3.0.1 

is  greater  than  or  equal  to  1 .6  for  single  mode  operatfon.^*"  In  Eq.  3.0.1  the 
free<space  propagation  constant  of  plane  waves  is  defined  as  ks  ZidX  or 

wave  number,  where  X,  is  the  wavelength.  Since  the  waveguides  fabricated 
for  this  experiment  will  be  rectangular  ridge  type  structures,  we  will  use 

several  of  the  dielectric  strip  waveguide  methods  to  model  the  devices  for 
this  experiment.*'^^^ 

Marcaffli  Rectangular  Dielectric  Waveguide  Model 
In  1969,  Marcatilf^  presented  an  analytical  model  for  integrated  optical 
rectangular  dielectric  waveguides  and  directional  couplers.  Fig.  3.1.1 
illustrates  a  directional  coupler  with  the  following  parameters:  n^  and  n2  are 
the  refractive  indices  of  the  core  and  top  clad,  respectively;  ng  is  the 
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refractive  irKtox  of  the  outside  dad;  04  is  the  refractive  irrdex  of  the  bottom 
dad;  r^  is  the  refractive  irxfex  of  the  center  ciad;  the  width  and  thickness  of 
the  individuai  waveguides  are  a  and  b,  respectively;  the  distance  of 
separation  between  the  two  waveguides  is  c.  The  modes  of  the  dieledric 
waveguide  are  classified  as  transverse  electric  (TEp^  arxi  transverse 
magnetic  O’MpJ  modes.  The  variables,  p  =  (0.1,2....)  and  q  =  (0,1,2,...), 
indicate  the  number  of  extrema  of  the  eledric  or  magnetic  field  along  the  x 
and  y  directions,  respedively.  This  gives  the  order  of  the  mode.  For  the 
fundamental  "0th  order*  TEp,  nrK>de,  p  =  q  =  0,  and  the  transverse  eledric 
mode  is  given  the  designation  TEqq.  TEp,  (TMp^  modes  have  the  eledric 
field  component  polarized  predominantly  in  the  y  (x)  diredion  and  the 
magnetic  field  component  polarized  predominantly  in  the  x(y)  diredion.  The 
amplitudee  of  the  fields  have  only  one  maximum  in  each  diredion. 

For  weH^uided  modes,  the  field  decays  exponentially  in  the  n2,  rig,  n4 
and  regions  of  the  diredionai  coupler  waveguide  strudure.  Most  of  the 
power  travels  in  the  n^  regions  while  only  a  small  part  of  the  power  travels 
in  the  r^,  r^,  r^  arxl  n^  regions.  An  even  smaller  portion  of  the  power  travels 
in  the  shaded  areas.  Marcatili  found  that  very  iittie  error  was  introduced 
when  the  fields  along  the  edges  of  the  shaded  regions  were  not  properly 
matched  when  calculating  the  fields  in  the  n^  regions.  This  allowed  him  to 
negled  the  shaded  areas  from  consideration,  thus,  simplifying  the  boundary 
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value  problem  by  solving  only  along  the  four  sides  of  the  core  region  n^. 


TEpq  Modes 

For  the  rectangular  dielectric  waveguide,  illustrated  in  Fig.  3.1 .1 ,  the 
electric  and  magnetic  field  components  in  the  five  areas  are  given  b/' 

=  M|OOS(/cp(  •¥  a)cos(/^  +  S)exp(-ikgZ¥ktiti,  3.1 .1 

M2006(^+ a)exp[>(M^/t-il^n>0].  3.1.2 

Hga  s  M^o06(^  4-  hyex^^iik^SCgZ-koij],  3.1 .3 

M  M4COs(^  *f  ajea^iiky^y-ikjZ^katj,  3.1 .4 

Hjg  s  M^oos(^  4-  5)sin(/^  4-  ‘y)exp(-Mr^iQoO.  3.1 .5 

Hy,  =  =  Hytf  =  0,  3.1.6 

3.1.7 


Ha’ 


3.1.8 


82 


3.1.9 


E^  « -iWfsai, 

Ea  *  -[(1/(<an2*/g)3*H^>B>^, 
Ea  =  -[(1/(<i*n3*/g)3*H^>0>5. 
E^«-[(1/(<Ben/ig)3*HjaA0>1, 
^.-[(1/(<aen,*^)a*fV9x9y!. 
Ey,  =  fy(**n,*  -  k/)/iaen,^k^l 
Ey,^HJi{/^n^r-k/)/<xn^%], 


3.1.10 

3.1.11 

3.1.12 

3.1.13 

3.1.14 

3.1.15 

3.1.16 

3.1.17 

3.1.18 


3.1.19 
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£y4  -  3.1 20 

Eys  =  -  kyl)/(aEX\^k^l  3.1 .21 

a  [(itol^*)dHjrt/3>1,  3.1 .22 

Ea  »  KI^<aen2“)dH^>l,  3.1 .23 

Ea^[{9men^^HJByi,  3.1.24 

■  [(#a)er»4*)dH^>^.  3.1 .25 

»  Ktons*)dHV9>J.  3.1 .26 


wheie  are  the  amplitudes  of  the  field  in  each  respective  medium,  a 
localae  the  field  maxima  in  the  core  region  n,,  5  locates  the  field  minima  in 
the  core  region  n,,  <D  is  the  angular  frequency  arxJ  £  is  the  permittivity  of  free 

apace.  The  y  term  is  necessary  because  single  mode  operation  of  a 
cfiredionai  coupler  will  guide  two  kinds  of  TEqq  modes.  One  is  symmetric 

when  y  m  9Cf,  and  the  other  is  antisymmetric  when  y  s  Both  are  TEM 

modes  with  main  field  components  Ey  and  H^.  For  a  symmetric  mode,  the 

_ . • _  _ _ 


plane  x  «  0  is  a  magnetic  short  ctrcuit.  For  an  antisymmetric  mode,  the 
plane  x  »  0  te  an  electric  short  drcuK.  The  electric  arnl  magnetic  field 
components  for  t>oth  symmetric  and  antisymmetric  modes  are  depicted  in 
Fig.  3.1.2. 

The  propagation  constants  k^  in  each  of  the  five  regions  are 

related  Is/^ 

k^  ky^  kg  =  cQpepn/  =  k^,  3.1  JZ7 

where  vs  i ,  g,  3, 4, 5  indicates  the  region  and  p  is  the  permeability  of  free 
space.  In  order  to  match  the  fields  along  the  boundaries  between  the  core 
region  n^  and  the  top  and  bottom  dad  regions  (ns  and  nj  and  between  the 
core  region  r^  and  the  outside  and  center  dad  regions  (ng  and  n^,  Marcatili 
assumed  the  following: 

kgi»  k^^  kj^^  kg  3.1 .28 

and 

k^s  k^rt  k^ss  ky.  3.1 .29 

Sinoe  ttte  refractive  index  of  the  core  region  n^  is  larger  than  the 
refractive  indices  of  the  cladding  regions,  only  modes  of  plane  wavelets 


I  Ey  or  Hx  for  TEoo 


Symmotric 

j  Ey  or  Hx  for  TEoo 


Figure  3.1^ 

Field  Intensity  Profiles  of  the  Symmetric  and 
Antisymmetric  Modes  in  a  Rectangular  Directional  Coupler 
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impinging  at  grazing  angles  on  the  surface  of  the  core  region  n^  are  guided. 
This  implies  that  kg  <  aruJ  ky  <  k^  and,  therefore,  the  electric  field 
components  and  can  be  neglected.  Matching  the 

remaining  field  components,  and  Ey,  along  the  boundaries  of  the  core 
region  n,  gives  the  transcendental  equations*' 

M  +  (2^a)exp(-o^-C7)/(1  +  kj^^]  3.1 .30 

and 

=  (Q  +  1  )jt  •  tan*^[(n2*/nt*)(^^2)]  -  tan''[(n//h,*)(/9i4)],  3.1 .31 

where  the  transverse  propagation  constants  k^  and  ky  are  solutions  to  Eqs. 
3.1 .30  and  3.1 .31 ,  respectively,  k^  is  the  solution  to*' 

=  (p  +  1)jc  •  \an’\kg^  -  tan*'(^jjy .  3.1 .32 

A  measure  of  the  penetration  depths  of  the  field  components  in  the  various 
media  are  given  b/' 

a  -  >5.*]“  3.1.33 

and 


11m  -  -  ^1“ 


3.1.34 
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where 

-  rVa.45^’^  3.1 .35 

are  the  maximum  thicknesses  of  the  media  for  which  the  slab  supports  only 
the  fundamental  TEqo  mode. 

From  Eqs.  3.1 .27. 3.1.28  and  3.1 29  the  axial  propagation  constant  for 
each  mode  of  the  coupler  is  determined  by*' 

^  3-1  -36 

where  s  kn,  is  the  propagation  constant  of  a  plane  wave  through  a 
medium  with  a  refractive  index  n^  at  a  free-space  wavelength  X,  and  and 
ky  are  obtained  from  Eqs.  3.1.30  and  3.1.31 . 

Substituting  y  into  Eq.  3.1.30  gives  slightly  different  axial  propagation 
constants  for  the  symmetric  (y  s  90^  modes  and  antisymmetric  (y  s  (f) 

iTKxles.*' 

/^  =  +  2(/i,’O(Va)exP<-0^<l  +  3.1 .37 

and 


3.1.38 
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where 

^  3.1 .39 

is  the  axial  propagation  constant  of  the  TE  mode  of  a  single  dielectric 
waveguide  (c  -¥  «>)  (see  Fig.  3.1 .2.). 

The  interaction  length  L  for  total  power  transfer  from  one  waveguide  to 
ttie  other  and  the  coupling  coefficient  K  between  the  two  waveguides  can 
then  be  related  to  k„  and  b/' 


•iK=n/{2L)^{k„-kJ/2 

-  3.1 .40 

The  interaction  length  and  coupling  coefficient  for  the  zero-gap  directional 
coupler  are  obtained  by  driving  the  separation  between  the  two  adjacent 
waveguides  to  zero,  which  yields 

-/Af*  jU(2L)  =  2{kJlk^{^aW  +  3.1 .41 

TMp,,  Modes 

The  analysis  for  the  TMp,  modes  is  very  similar  to  that  for  the  TEp, 
fTKxJes.  The  electric  and  magnetic  field  components  for  TMp^  modes  are 
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found  by  substituting  Efor  H,H  for  E,  -e  for  \i  and  -[i  for  e  into  Eqs.  3.1 .1  • 
3.1 .26.  Foiiowing  Eqs.  3.1 .28  arKi  3.1 .29,  the  field  components  and 
can  be  neglected  for  the  analysis. 

Matching  the  field  components,  E^,  and  Hy,  along  the  boundaries  of  the 
core  region  n^  in  the  same  way  as  was  done  to  match  the  field  components 
Hji  and  Ey  for  the  TE„  analysis  gives  the  transcendental  equations®' 


V 


it 


M  *  ^xoa[1  +  (2Va)exp(.G^-2fr)/(1  +  3.1 .42 

and 

^  =  (q  +  1)ic  •  Xan'\k/\^  -  tan*'(/9i4),  3.1 .43 

where  the  transverse  propagation  constants  and  ky  are  solutions  to  Eqs. 

3.1 .42  and  3.1 .43,  respectively,  is  the  solution  to 

M  »  (P  +  1)«  -  tan*’[(n3'/ni^(/c^]  ■  3.1 .44 

and  1)2, 1)4,  ^  and  ^  are  obtained  from  Eqs.  3.1 .33,  3.1 .34  and  3.1 .35. 

The  interaction  length  L  for  total  power  transfer  from  one  waveguide  to 
the  other  and  the  coupling  coefficient  K  between  the  two  waveguides  are 
determined  by  Eq.  3.1.40  using  ky  and  ^xO  obtained  from  Eqs.  3.1 .39, 

3.1.43  and  3.1.44,  respectively.  The  interaction  length  and  coupling 


a 
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ooefficieiit  for  the  zero-gap  directional  coupler  are  obtained  by  substituting 
ky  and  k^  obtained  from  Eqs.  3.1 .39, 3.1 .43  and  3.1 .44,  into  Eq.  3.1 .41 . 

Marcuse  Rectangular  Dielectric  Waveguide  Model 
Marcuse^  has  derived  a  rectangular  dielectric  waveguide  model 
following  the  approach  developed  by  Marcatlll.  Like  the  preceding  analysis, 
the  four  shaded  regions,  (see  Rg.  3.1.1),  are  neglected.  The  analysis 
presented  by  Marcuse  starts  with  the  Maxwell  equations: 

V  X  H  =  ejfi^dE/dt  3.2.1 

and 

VxE^^dH/dt,  3.22 

where  H  is  the  magnetic  field  vector,  E  is  the  electric  field  vector,  £„  is  the 
dielectric  permittivity  of  vacuum,  is  the  magnetic  permeabiiHy  of  vacuum, 
n  is  the  refractive  index  of  the  medium  and  t  is  the  time  variable.  The  time 
and  2  dependence  equation  for  the  time  harmonic  and  mode  fields  is  given 


exp[^(off2)], 


3.2.3 
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vvhera  (d  is  the  radian  frequency  and  the  eigenvaiue  p  is  the  z  component 
of  the  propagation  vector  of  plane  waves  traveling  in  the  waveguide  core. 

The  variable  is  designated  as  the  propagation  constant  of  the  guided 
waveguide  modes  and  is  expressed  by^ 


p  s  nilccos6i, 


3.2.4 


where  k  *  2itA,  is  the  wave  number,  n^  is  the  refractive  index  of  the 
waveguide  core  and  6^  is  the  angle  at  which  the  wave  propagates  through 
the  waveguide  core  with  respect  to  the  waveguide  dad  ns  (see  Fig.  3.2.1). 
For  the  wave  to  propagate  through  the  core  medium  with  a  refractive  index 
of  n^,  e,  must  not  exceed  the  critical  angle  for  total  internal  reflection  given 


cose,c  =  Vn,.  3-2-5 

Using  the  equations  above,  the  transverse  electric  and  magnetic  field 
components  for  the  TE„  and  'm„  modes  can  be  given  in  terms  of  the 
longitudinai  field  components  by^ 

E,  -  -<#K,*)((p3E;ax)  +  3.2.6 
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E,  .  35.7 

H,  « -(fl^*)((p3H^j4  -  3.2.8 

and 

H,  =  -(«<,*)l(p3H^))  +  (Bin|\3E^x)],  35.9 

where 

Kj  =  3^.10 

and  j  s  1 ,  2,  3,  4,  5  corresponds  to  the  values  for  the  fields  in  the  five 
regions  of  the  waveguide.  The  longitudinal  field  components  and  must 
then  satisfy  the  wave  equation" 

(aV9J^  +  (3V9/)  +  KjV  =  0.  3.2.1 1 

TE^  Modes 

For  TEp,  modes,  the  electric  field  component  is  polarized  predominantly 
in  the  y  direction.  The  amplitude  coefficients  of  E,  and  are  adjusted  so 
that  one  of  the  transverse  field  components  vanishes.  The  following  set  of 
longitudinal  field  components  satisfies  the  wave  equation  (Eq.  3.2.1 1)  and 
the  trmverse  field  component  equations  (Eqs.  3.2.6  •  3.2.9).  They  describe 
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a  TE„  mode  in  the  core  region 

Eg  s  Acoe[^^(/fT\)]cos[/^x-44)]  3.2.12 

and 

H,  =  -A(^Mj^*n,*(Wfj/p/g8in[/^Ti))sin[/f^  3.2.1 3 

with 

=  n^*A*  -  P*  =  /^*  +  k,^;  3.2.14 

and  that  to  the  top  dad  region  r^  yields^ 

Eg  »  Ac08{lyf»lC06[^x  +  ^lexp<*72y)  3.2.15 

and 

Hg  =  3.2.16 

with 

-  P*  =  /^*  -  Y2*;  3.2.17 

and  that  to  the  outside  dad  region  yields*^ 

Eg  «  A(n,*M3*)coa[<g(»f^]co8C/^(>^  32.1 8 
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and 

H,  .  4A(e^^nA<fy3/p/gco8[iaa+^)l 

•  8ln[l^(y«i)lexp[-Tf3(jf-a)l], 


with 


and  that  to  the  bottom  dad  region  04  yields*^ 


^  «  Acoe[/^n-^]a)sIMJfH)]exrtY4(>«-b)] 
and 

•  8in(MJ«^]exp[Y4(y*'i>)D. 

with 

K**.n4*<?-p*«/^*-T4*; 
and  that  to  the  center  dad  region  r%  yields^ 


£,  «  A(n,X*)oo8(/^^5)<X)e{^Tt)Joxp<y^ 


3^.19 


3.2^ 


3^^1 


3^J22 


3.2^3 


3^^4 


and 
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with 

K6**ne*A*-p*«/^*-Y5^  3^.26 

In  Eqs.  3.2.12. 3.2.13,  3.2.15,  3J2.16,  3.2.18,  3.2.19,  3.2.21,  3.2.22,  3.2.24 
and  3.2.25,  variable  A  is  the  amplitude  of  the  field. 

The  tranevefse  electric  field  component  Ey  can  satisfy  one  of  the 
boundary  conditions  by  using  the  proper  field  amplitudes.  The  other 
boundary  oorxitlons  can  be  matched  by  forcing  the  longitudinal  magnetic 
field  component  H,  in  the  top  and  bottom  dad  regions  (n^  and  nj  to  pass 
continuously  through  the  core  boundary  at  ys  0  arfo  -b,  dictating  that 
be  continuouB.  Note  that  Exit’d  neglected  sinoe  it  is  small  compared  to 
the  other  field  components.  These  conditions  lead  to  the  equations*^ 

- 1^]  -  (n/ZiKjoosI/^it  -  b)]  *  0  35.27 

and 

(n//lp8in(^)  +  (ri2‘/i&)006<l9r))  =  0.  3558 

Expansion  of  the  sine  and  cosine  functions  in  the  above  equations  yields 
the  e^ienvalue  equabon*^ 


tttn{k^) « •  n,  VJ- 


3^^ 


This  is  the  eigenvalue  equation  of  TMp,  modes  in  an  infinite  slab  waveguide 
for  which  the  transverse  electric  field  is  normal  to  the  core  boundaries  at 
ys  0  and  y«  -b.  Using  Eqs.  3^.17  and  3^^3  to  express  Ys  and  Yi  in  terms 
of 


%  *  [(nt*  -  n**)^  -  Vl' 


3^.30 


y4-C(ni*-n4V-^?J^®. 


3.2.31 


one  can  detennine  k,  using  Eq.  3.2.29. 

At  the  waveguide  core  boundaries  in  the  outside  and  center  dad  regions 
(rt,  and  n^,  Hg  must  be  equal  and  continuous  on  both  sides  of  the  dielectric 
imeifaoe  at  x>  a  and  xs  0.  The  transverse  eiectric  field  component  Ey  is 
already  continuous  by  the  choice  of  field  amplitudes.  This  leads  to  the 
foiowing  equations^ 


-  Ys«»l^«0]  =  0 


3.3.32 


kMkJO  -  YsOOS(/^  =  0. 


3.2.33 
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Expanskm  of  the  sine  and  cosine  functions  in  the  above  equations  yields  the 
eigenvalue  equation'^ 

tan(M)  •  -  YaYs).  3^.34 

for  TEp,  modes  in  an  infinite  slab  waveguide.  Since  Ey  is  the  dominant 
electric  field  component  the  field  appears  as  a  TEp,  mode  with  respect  to  the 
oulaida  and  canter  dad  regions  n,  and  ng. 

Using  Eqs.  3.2.20  and  3.2.26  to  express  y,  and  Ys  in  terms  of  ky, 

Ys  «  [(Hi*  -  n,*)A*  -  kyT  3.2.35 

and 

Ye  *  [(n*  -  3.2.36 


one  can  dstermine  using  Eq.  3J2.34. 

The  propagation  constant  §  can  be  obtained  for  TEp,  mode  propagation 
using  Eq.  3.2.14,  with  the  values  for  ky  and  k^  determined  from  Eqs.  3.2.29 
arxf  3.2.34,  respectively.*' 


p  -  (n,V  -  {k*  +  V)l” 


3^37 
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The  interaction  length  L  for  total  power  transfer  from  one  waveguide  to 
the  other  and  the  coupling  coefficient  K  between  the  two  waveguides  are 
then  found  by  substituting  p  for  k^tor  and  l/yg  for  ^  in  Eq.  3.1.40. 
This  yields® 

k/{2L)  «  2(/i*/p)[(1/ays))exp(-cy6)/(1  +  3^.38 

The  interaction  langth  and  coupling  coefficient  for  the  zero-gap  directional 
coupler  are  obtained  by  driving  the  separation  between  the  two  adjacent 
waveguides  to  zero  which  yields, 

-«C«  it/{2L)  «  2(/<,*/p)(1/ay5y(1  +  i<*/ys*).  3.2.39 

TMp,  Modes 

The  TMp,  modee  are  derived  in  a  similar  manner  to  the  TEp,  modes,  with 
the  electric  field  component  polarized  predominantly  in  the  X  direction.  The 
folmMng  set  of  longitudinal  field  components  satisfies  the  wave  equation 
(Eq.  3.2.11)  and  the  transverse  field  component  equations  (Eqs.  3.2.6  • 
3.2<.9).  They  descitie  a  TMp,  mode  in  the  core  region  n^® 

B006[/9(>*H|)]C0S{^XfQ]  3.2.40 
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and 


and  that  to  the  top  dad  region  n2  yielde” 


and 

and  that  to  the  outside  dad  region  ns  yiekfe" 

Bco^kJiat^YiosikJiy^^^ 

and 

H,  » -{B(e^^  V(*V^3)co8C/aa^)l 
•  8inI«^(>>Hi)]exp{-Tr3(x-a)ll; 


and  that  to  the  bottom  dad  region  n4  yields” 


E^^B{ni^fn2^)oos[kJi^-biyxis[^^^ 

and 


3^.41 

3^.42 

3^.43 

3^.44 

3^.45 

3.2.46 
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•  sinCI^JH^lexrtTiCyfb)]];  3^.47 

and  that  to  the  center  dad  region  yields^ 

£g  s  Bcos(/t;,4)oo®[M>^’l)l®^Tfe^  3.2.48 

and 

Hr  *  B(ayMj^%*(W^6)ooe(^^8ln[<^(>^  3.2.49 

In  Eqa.  3.2.40  -  3.2.49,  variable  B  i$  the  ampHtude  of  the  field.  The  values 
for  K1^  K,^  K3^K4*  and  Ks^  are  obtained  by  using  Eqs.  3.2.14, 3.2.17, 3.2.20, 
3.2.23  arxJ  3.2.26,  respectively. 

The  transverse  electric  field  component  can  satisfy  cne  of  the 
boundary  conditions  by  the  proper  choice  of  the  field  amplitudes.  The  other 
boundary  conditions  can  be  matched  by  forcing  the  longitudinal  magnetic 
field  component  in  the  top  and  bottom  dad  regions  (ng  and  nj  to  pass 
continuously  through  the  core  boundary  at  o  and  ys  -p,  dictating  that  Hy 
be  continuous.  Note  that  ^  can  be  negleded  since  it  is  small  compared  to 
the  dher  field  components.  These  conditions  lead  to  the  equations” 

k^kjib  +  n)l  -  y^cosikjib  +  q)]  =  0 


3.2.50 
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arid 

k^slin{k^)  ■¥  y^D08{k^) »  0.  3.2.51 

Expansion  of  the  sine  and  cosine  functions  in  the  equations  above  yieids  the 
eigenvalue  equation*^ 

\an{k/>)  =  +  Y4)/(/^*  -  Y2Y1).  3.2.52 

for  TEp,  modes  in  an  infinite  slab  waveguide.  Using  Eqs.  3.2.30  and  32.31 , 
which  express  and  Y4  in  terms  of  ky, 
one  can  determine  ky  using  Eq.  32.52. 

At  the  waveguide  core  boundaries  in  the  outside  and  center  clad  regions 
(n,  arxi  n^,  must  be  equal  and  continuous  on  both  sides  of  the  dielectric 
interface  at  xa  a  and  xa  0.  The  transverse  electric  field  component  E,  is 
abeady  oontirtuous  by  the  choice  of  field  amplitudes.  This  leads  to  the 
fotknvirig  equatioris:^ 

(n,V«a*«4(5  -  -  (th^frdooslkja  -  41  -  0 

and 

(n  */fg8in(f^  +  (n,*/fc)cos(f^  =  0. 


3.2.53 

32.54 


Expansion  of  the  sine  and  cosine  functions  in  the  equations  above  yields  the 
eigenvalue  equation^ 

tan(M)  *  n,*/f^n3^s  +  -  n,  Vs),  3^.65 


for  the  TMf^  nKXles  in  an  infinite  slab  waveguide.  Thus,  the  field  appears  as 
a  TM„  mode  with  respect  to  the  outside  and  center  dad  regions  (ng  and  ng). 
Using  Eqs.  3.2.35  and  3.2.36,  which  express  73  and  Yg  in  terms  of  one 
can  determine  ^  using  Eq.  3.2.55.  The  propagation  constant  p  can  be 
obtained  for  TM„  mode  propagation  by  substituting  and  k^  obtained  by 
Eqs.  3.2.52  and  3.2.55,  into  Eq.  3.2.37. 

The  interaction  length  L  for  total  power  transfer  from  one  waveguide  to 
the  other  and  the  coupling  coefficient  /C  between  the  two  waveguides  are 
then  found  by  substitution  of  the  above  TM^  mode  variables  into  Eq.  3.2.38. 
The  interaction  length  and  coupling  coefficient  for  the  zero-gap  directional 
coupler  are  obtained  by  using  Eq.  3.2.39. 


In  a  1975  paper.  Bums  and  Milton  discussed  mode  conversion  in  planar 
(Select  separating  waveguides."  This  has  since  become  known  as  the 
effective  index  method,  and  is  treated  in  detail  by  Nishihara,  Haruna  and 


SuhaiB  for  both  two  and  three  dimensional  dielectric  waveguide  structures.^ 
The  effective  irKlex  method  is  based  on  the  concept  of  the  effective  index  of 
the  guided  mode.  The  guided  mode  propagating  along  the  z  axis  sees  the 
effective  index  {N)  defined  by^*^ 

P  =  foV,  3.3.1 

v^re  p  is  the  axial  propagation  constant,  /r  is  the  free  space  propagation 
constant  (defined  in  the  first  section  of  this  chapter)  and  N  =  n^cose^  (see 
Fig.  3J2.1  and  Eq.  3.2.4.). 

The  effective  index  method  for  the  dielectric  waveguide  structure  to  be 
fabrfoaled  for  this  experiment  (illustrated  in  Fig.  3.1.1)  is  performed  by 
dividing  the  three  dimensional  waveguide  into  two  two-dimensional 
waveguides.  One  is  2-0  Waveguide  1,  with  light  confinement  in  the  y 
direGtion  and  the  second  is  2-D  Waveguide  2,  with  light  confinement  In  the 
X  direction  (see  Fig.  3.3.1).  The  axial  propagation  constant  for  the  TEp, 
modes  of  a  three  dimensional  waveguide  structure  is  detemiined  by  the 
foNowing  procedure: 

The  effective  index  of  2-D  Waveguide  1  is  found  first  by  applying 
TMp,  mode  analysis  to  the  top  and  bottom  dad  regions  (ris  and  nj  with 
waveguide  thickness  b.  Effective  index  Al,  is  then  used  to  find  the  effective 
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index  of  2>D  Waveguide  2  by  applying  TEp^  mode  analysis  to  the  outside 
and  center  dad  regions  {n^  and  n^)  with  waveguide  width  a.  The  axial 
propagation  constant  is  finally  obtained  by  substituting  A4  for  /V  in  Eq.  3.3.1 . 

The  axial  propagation  constant  for  the  TMp,  modes  is  determined  in  a 
manner  similar  to  the  one  used  for  the  TEp,  modes.  First,  effective  index 
of  2-D  Waveguide  1  is  found  by  applying  TEp,  mode  analysis  to  the  top  and 
bottom  clad  regions  (n2  and  nj  with  waveguide  thickness  b.  Effective  index 
W,  is  then  used  to  find  effective  index  of  2-D  Waveguide  2  by  applying 
TMpq  mode  analysis  to  the  outside  and  center  clad  regions  (Oa  and  r^)  with 
waveguide  width  a. 


TEpg  Modes 

As  outlined  in  the  preceding  section,  2-D  Waveguide  1  represents  the 
TMp,  modes  in  an  infinite  slab  waveguide  of  thickness  band  2-D  Waveguide 
2  represents  the  TEp,  modes  in  an  infinite  slab  waveguide  of  width  a.  Since 
Ep  is  the  dominant  electric  field  component,  the  field  appears  as  a  TEp,, 
mode  with  respect  to  the  outside  and  center  dad  regions  (ng  and  ng).  The 
analysis  can  be  represented  by  the  following:^ 


Pa  =  /</V/TM*(nain,in4);TE*(n3l/V,ln5)}], 


3.3.2 
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where  is  the  TEp,  axial  propagation  constant  for  the  three  dimensional 
cHeiectric  waveguide. 

Rrst.  we  start  with  the  wave  equation  for  the  TMp,  modes  of  2-D 
Waveguide  1 ,  given  b/' 


3.3.3 


where 

Ey  s  3.3.4 

and 

E,  s  iVlmjr?)idHJBy,  3.3.5 

where  ^  is  the  axial  propagation  constant  for  2-D  waveguide  1. 
The  field  solutions  can  then  be  written  as*' 

Hx,  a  H,006(^  +  ^2)>  3.3.6 

at-b<y<0fbrthe  core  region  n^, 


3.3.7 


108 

at  y  >  0  for  the  top  dad  region  n^,  and 


Hj^*H4exp[Y4(n6)]. 


3.3.8 


at  y  <  •/>  for  the  bottom  dad  region  n4.  The  varieties  Yj  and  y*  are  the 
transverse  propagation  constants  of  a  piane  wave  in  the  top  and  bottom 
dadding  ragiona,  respedivaly,  and  ^2  >a  the  phase  shift  suffered  by  a  wave 
polarized  witti  its  magnetic  fieid  vedor  parallel  to  the  interface  between  the 
core  region  n^  and  the  top  dad  region  n^. 

In  addition  to  the  phase  shift  due  to  refledion  of  the  wave  from  the 
interface  between  n^  and  r\2>  there  is  also  a  phase  shift  ^4  suffered  by  the 
reflection  of  the  wave  from  the  Interface  between  the  core  region  n^  and  the 
bottom  dad  region  n4.  Using  Eqs.  3.3.3  •  3.3.8 ,  these  phase  shifts  can  be 
defined  maihematlcatty  by" 

^2  « -2arctari[(ni*/h2*)(P,*  -  •  p,*)'®]  3.3.9 

arKi 

04  a  -2arclan[(n|*/h4“)(P^*  -  -  Pi*)'®].  3.3.10 

Using  the  differences  between  the  phase  fronts,  which  are  due  to  the 
reflection  of  the  wave  at  the  two  interfaces  and  given  mathematically  by" 

_ - _ _  -  .  _ _  .  _  . . 
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8,  « (teo66,)[(1Aanet)  -  tanej  s  (£i^8in6i)(co8^i  -  sin^i)  3.3.1 1 
and 

$2  X  A^slnei,  3.3.12 

the  condition  that  t>oth  reflected  waves  contribute  to  the  same  plane  waves 
can  be  expressed  by  the  relation*^ 

n,(82  -  8i)k  +  02  +  ^4  *  2(a  +  1)jt,  3.3.1 3 

where  q  x  (0,1.2....)  indicates  the  number  of  extrema  of  the  magnetic  field 
in  the /direction. 

The  eigenvaiue  equation,  obtained  from  Eqs.  3.3.9  •  3.3.13,  can  then  be 
written  as^ 

^  X  (g  +  I)ic  -  tan'’((n2*/h^*)(^^2)J  -  tan*H(n//hi*)(Vy4)J,  3.3.14 

where 

3.3.15 

Y2  =  3.3.16 

and 
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3.3.17 


are  the  tranaveree  propagation  constants. 

Using  Eq.  3.3.1,  the  transverse  propagation  constants,  given  by  Eqs. 
3.3.15  -  3.3.17,  can  be  expressed  in  terms  of  the  effective  index  of  2-D 
Waveguide  1: 

3.3.18 


Ti  «  .  na*)"®  3.3.19 

and 

3-3-20 

Next,  the  wave  equation  for  the  TEp,  modes  of  2>D  Waveguide  2  is 
givetYb/' 


a*^j^  +  (/W-pa*)fy  =  0, 


3.3.21 


where 


3.3.22 


and 
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3.3^3 


The  field  eolutione  can  be  written  as*' 

Eyt  •  £,coe(A^  *  4  3.3^4 

at  0  <  X  <  a  for  the  core  region  now  defined  as 

^  3.3^5 

at  X  <  0  for  the  omeide  dad  region  n„  and 

^  3.3^ 

at  X  >  a  for  foe  center  dad  region  n,.  The  variabies  7^  and  75  are  foe 
transverae  propagation  constants  of  a  plane  wave  in  foe  outside  and  center 
cladding  regione,  reepectively,  and  ^3  is  foe  phase  shift  suffered  by  a  wave 
polarizad  with  la  electric  field  vedor  parallel  to  the  interface  between  foe 
core  region  At,  and  foe  outside  dad  region  ns. 

In  addWon  to  foe  phaee  shift  due  to  reflection  of  the  wave  from  foe 
intaifooe  between  At,  and  ns,  there  is  also  a  phase  shift  suffered  by  the 
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nHh&c/Son  of  the  wave  from  the  iriterface  between  the  core  region  and  the 
oemer  dad  region  r^.  Ueing  Eqe.  3.3^1  -  3.3^  .  theee  phase  shifts  can 
be  defined  rnathernatically 

« -2arctan[(fc*  -  r\,***)'"/(n,*/^  -  3.327 

and 

4s « -2afclan{(ft,*  -  ns*/t*)"*/(n,*i((*  -  3.328 

Using  the  dfiferences  between  the  phase  fronts,  which  are  due  to  the 
refiectton  of  the  wave  at  the  two  interfaces  and  given  mathematically  b/^ 

s, « (acoae^HO^anfi,)  •  tane^J « (a%lne,)<co8^i  •  dn^Ot)  3.3.29 
and 

8^  «  aMnei,  3.3.30 

ttie  oondMon  that  both  refiected  waves  contribute  to  the  same  plane  waves 
can  be  ergaeeeed  with  the  relation*^ 

nt(Ss  •  9i)k+  ^$*^g^2(p+  1)a,  3.3.31 

where  p  « (0,1.2,...)  Indicates  the  number  of  extrema  of  the  electric  field  in 


the  xcRredion. 


The  eigenvalue  equation,  obtained  from  Eqs.  3.3^  -  3.3.31 ,  can  then 
be  written  as*' 

M  *  (P  'I  )k  *  tan*^(/^3)  -  tan*^(/fyy5),  3.3.32 

where 

3.3.33 

*&  *  (ft**  -  h**^“  3.3.34 

and 

Y.«(ft**-ni*A?)'"  3.3.35 

are  the  transveree  propagation  constants. 

Using  Eq.  3.3.1,  the  transverse  propagation  constants,  given  by  Eqs. 
3.3.33  -  3.3.35,  can  be  expressed  in  terms  of  the  effective  index  of  2-0 
Waveguide  2: 

kiN,^  ‘  3.3.36 


3.357 
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and 

75  -  {Ni*  •  ns*)"".  3.3.38 

Having  aolved  for  and  the  symmetric  propagation  constant  and 
antisymmetric  propagation  constant  corresponding  to  the  symmetric  and 
antisymmelric  TEqq  modes  of  the  directional  coupler,  respectively,  are  found 

by- 


P.  -  W1  +  a(/i*/p,*)(1/«lfc)exp(-crJ/(1  +  3^.39 

and 

P.  *  wi  •  2«^*/p,»)(1/ays)exp(-cysy(1  +  3.1 .40 

where  c  is  9ie  dislanoe  of  separation  between  the  two  waveguides. 

The  interaction  length  for  total  power  transfer  from  one  waveguide  to 
ttie  other  and  the  coupling  coefficient  K  between  the  two  waveguides  can 
then  be  rslaied  to  and  by*' 


-IK*  it/{2L)  *  (p.  -  PJ/2.  3.3.41 

The  interaction  length  and  coupling  coefficient  for  the  zero-gap 


dirsctional  coupler  are  obtained  by  driving  the  separation  between  the  two 


aclMfit  waveguides  to  zero.  [Thie  yields, 

I 

P. -  W1  +  +  k,‘frs))]  3.3.42 

and 

P.  «  W1  -  2(/^*/p,*)(1/ay8)(1/(1  +  /c>5*))],  3.3.43 

which  are  then  substituted  back  into  Eq.  3.3.41 . 

TMp,,  Modes 

For  mode  analysis,  2<0  Waveguide  1  represents  the  TE„  modes 
in  an  inflnile  stab  waveguide  of  thickness  b  and  2-D  Waveguide  2  represents 
the  modes  in  an  infinite  slab  waveguide  of  width  a.  With  as  the 
dombiant  magnetic  field  component,  the  field  appears  as  a  'TM„  mode  with 
laepoct  to  the  outside  and  center  dad  regions  (ri,  and  rig).  The  analysis  can 
be  represented  by  the  fbilowing:^ 

P,  .  <(AMTE*(n,lt»,ln4):TMnn,IN,lnJ)J.  3.3.44 

where  is  tfie  TM„  axial  propagation  constant  for  the  three  dimensional 
dMectrle  waveguide.  Following  the  same  procedure  as  was  presented  in 
tie  previous  section,  the  eigenvalue  equation  for  the  TEp,  modes  of  2-D 
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Waveguide  1  can  be  written  as*' 


M  *  (P  +  1  )J'  ■  tan*'[(n3*/ni*)(/0'3)]  -  tan'^[(ns*/n,*)(Mys)]»  3.3.45 

> 


where 

3.3.46 

%  =  Oi*  -  3.3.47 

and 

Tk  -  3.3.48 


Ueing  Eq.  3.3.1,  the  transverse  propagation  constants,  given  by  Eqs. 
3.3.46  •  3.3^,  can  be  expressed  in  terms  of  the  effective  index  of  2-D 
Waveguide  1: 


3.3.49 


and 


3.3.50 


%  -  W*  -  n,*)’*. 


3.3.51 
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The  eigenvalue  equation  for  the  TMp^  modes  of  2’D  Waveguide  2  can 
be  written  as*' 


*  (q  +  1)ii  -  tan’^CVra)  -  tan'i/f/yj, 


3.3.52 


where 


3.3.53 


and 

T.  *  (P2*  - 


3.3.54 

3.3.55 


Using  Eq.  3.3.1,  the  transverse  propagation  constants,  given  by  Eqs. 
3.3.53  -  3.3.55,  can  be  expressed  in  terms  of  the  effective  index  of  2-D 
Waveguide  2: 


3.3.56 


Ya  =  3.3.57 

and 

Y4  *  (^4**  ■  3.3.58 


Having  solved  for  ^  and  the  symmetric  propagation  constant  and 
antisymmetric  propagation  constant  for  the  TMgo  modes  of  the  directional 
coupler  are  found  by  substituting  ^  and  determined  here,  into  Eqs.  3.3.39 
and  3.3.40. 

The  interaction  length  L  for  total  power  transfer  from  one  waveguide  to 
the  other  and  the  coupling  coefficient  K  between  the  two  waveguides,  can 
then  be  related  to  and  by  Eq.  3.3.41 .  The  interaction  length  and 
coupling  coeff^^^ent  for  the  zero-gap  directional  coupler  are  obtained  using 
'  3.3.42and  43.  which  are  then  substituted  back  into  the  equation  for 

the  zero-gap  directionai  coupler,  Eq.  3.3.41. 


Fleck,  Morris  and  Felt  proposed,  in  1976,  a  propagating  beam  method  for 
computing  fields  of  an  optical  beam  passing  through  a  medium.^  The 
method  calculates  the  characteristics  of  an  input  beam  as  it  propagates 
through  a  medium  over  a  small  distance  z,  and  then  corrects  for  the 
variations  of  the  refractive  index  as  seen  by  the  beam  over  this  distance  z. 
it  has  already  been  applied  to  several  optical  fiber  and  integrated  optic 
applications.®^^’"’* 

The  beam  propagation  method  is  derived  for  a  scalar  field.  The  theory 
is  restricted  to  small  changes  in  the  refractive  index.  The  first  part  of  the 


119 

derivation  assumes  the  propagation  of  a  high  frequency  beam  through  an 
inhomogeneous  medium.  It  begins  with  the  wave  equation^^ 

a  0,  3.4.1 

where  ^  represents  the  scalar  field,  k  is  the  free  space  wave  number  defined 
at  the  beginning  of  this  chapter  and  n(r)  a  noCr)^  +  An^(r)  is  the  refractive 
index  of  the  medium  consisting  of  an  unperturbed  part  (no(r)^  and  a 
perturbation  (An^(r)).  Splitting  n^(r)  into  its  respective  parts,  Eq.  3.4.1  takes 
the  form^** 

+  A*no*(r)0  =  -k*An*(rX> «  s(r),  3.4.2 

where  the  right  hand  side  of  the  equation,  designated  as  s(r),  is  considered 
a  source  function.  Next,  the  unperturbed  part  of  the  index  is  chosen  so  that 
tfie  new  wave  equation^^ 

VSy  +  /?r\,*(r)v  a  o,  3.4.3 

along  with  the  radiation  conditions  at  za  o»,  has  a  solution.  If  ly  is  known  for 
z  a  ^  where  ^  is  the  position  z  a  o,  then  y  and  its  derivative,  with  respect 


to  z.  can  be  solved  for  ail  values  of  zby  use  of  an  operator  (d).  d  acts  with 
respect  to  only  the  transverse  coordinates  This  leads  to^" 

dv/3z=  3.4.4 

For  i\  s  constant,  the  scalar  field  can  be  represented  in  the  angular 
spectrum  domain 

4>oe 

=  //  v(x,y,2)expl-^fyf  +  f^dxdy,  3.4.5 

•oo 

where  is  the  Fourier  transform  of  the  scalar  field  (y(x,y,2)),  and 

and  are  the  transverse  propagation  constants.  Taking  the  derivative  of 
Eq.  3.4.5,  with  respect  to  z,  yields^^’” 

3.4.6 

where  kg  s  •  k^^^  is  the  axial  propagation  constant.  Using  Eq. 

3.4.6,  the  integral  operator  6  can  be  derived. 

For  a  given  coordinate  (2),  the  field  (^)  is  split  into  two  parts.  The  first 
part  (^f)  is  generated  by  the  sources  in  the  region  where  z'  <  z,  and  the 
second  part  (^9)  generated  by  the  sources  where  z'  >  z.  Using  the 
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function:^^ 


»  0  for  z  <  z' 
o^frlO  =  1/2fbrzaz' 

e^(rlr')  =  1  for  z  >  zl,  3.4.7 


an  explicit  expression  for  and  ^2  can  be  obtained.  Defining  G  as  Greens 
function  of  Eq.  3.3.3,  the  wave  equation  can  be  expressed  b/'* 

♦,(r)  =  // /  Q(rir')e,(rlr')s(r')dV'.  3.4.8 

Taking  the  derivative  ot  f  „  with  respect  to  z,  yields” 


=  3Q02(rlr')e,(rlr')s(r')dV' 

-oo 

+  ///e(rlr-)8(z-y)s(ndV'. 

•OP 


3.4.9 


where  5(z  >  zO  is  the  step  size.  The  first  integral  of  equation  3.4.9 
represents  the  unperturbed  medlizn.  It  can  be  expressed  in  terms  of  fusing 
Eq.  3.4.4.” 


3.4.10 


122 

The  second  integtai  of  Eq.  3.4.9  includes  both  the  unperturbed  medium  and 
the  perturbation.  With  the  help  of  Eq.  3.4.2,  the  second  integral  in  Eq.  3.4.9 
can  be  written  as^^ 

/]^(x,y,2Ur,/,2)(-A*)An*(y^,y\2)<Kx',/,2)dx'd/.  3.4.1 1 

•ee 

Equcto)  3.4.1 1  can  then  be  expressed  in  terms  of  an  operator  (0)  acting  on 
♦  ty" 

34.12 

where,  like  A,  operator  0  acts  with  respect  to  only  the  transverse  coordinates 
(x;y).  Substitution  of  Eqs.  3.4.10  and  3.4.12  into  Eq.  3.4.9  yields^'* 

=  A^|(x',/,4  +  O0(x',/,2).  3.4  1 3 

Neglecting  the  influence  of  the  reflected  flekJ  (^2)  ^  ^1*  ^ 

substituted  for  0  in  Eq.  3.4.13.  The  first-order  differential  equation  for 
then  becomes^'* 


d0i^2(x',/,2)  =  fi0,(x',/,2)  +  Q<|>,(x',/,;^ 


3.4.14 
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where  0^1  is  the  correction  term. 

From  Eq.  3.4.14,  the  propagation  of  field  (^0  is  Qiven  by  two  terms. 
One  describes  propagation  in  an  unperturbed  medium  (di>0  while  the  other 
describes  a  correction  term  representing  the  influence  of  An.  Using 
Fourier  transforms  and  convolution  along  with  the  assumption  that  the 
unperturbed  media  is  homogeneous,  the  correction  term  in  Eq.  3.4.14  can 
be  simpHfied.^^ 


=  -(iW2no)An*<|.,(X'./.i).  3.4.15 

By  combining  Eqs.  3.4.14  and  3.4.15,  the  first-order  differential  equation  for 
becomes^^ 

=  S<>,(x',/.4  -  (»iaio)An*4,(x',/.4.  3.4.16 

In  order  to  solve  Eq.  3.4.16,  it  will  be  assumed  that  is  known  at  position 
Zb  2^.  If  An^  a  0,  Eq.  3.4.16  reduces  to" 

delBz(KyfZf^  -  5e(x,y,^,  3.4.17 


where  z{x,y,z^  =  is  given.  Therefore,  e  represents  the  field  as  it 
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profMiQales  thnxigh  a  madium  with  rafracth/e  indax  (no(r)).  Introducing  a 

oorradion  factor  (axpCT))  to  an  approximata  solution  to  tha  diffarantial 
Eq.  3.4.16.  This  solution  is  darivad  b/'* 

♦,(r) «  axp[r(r))e(r),  3.4.18 

with  T{x,y^  «  0.  Combining  Eqs.  3.4.16  and  3.4.18  yialds^^ 

dr/Bz  *  +  [6(axp(r)e)  >  axp(r)6(e)yoxp(r)e.  3.4.1 9 

K  has  boon  damonstratad  that  the  second  tarni  in  Eq.  3.4.19  can  be 
neglactad,  provided  Eq.  3.4.19  is  integrated  over  a  small  cMstanca  (z  •  z^7* 
TNsyMds 

r  « ’iU2n^^{z  •  Zi)f  3.4.20 

for  (z  -  ^  amali.  Combining  Eqs.  3.4.18  and  3.4.20  yieids^^ 

A2)  a  e(x,y^  -i-  A2)axpC>(ii^2no)An^A2),  3.4.21 

which  is  tha  aquation  for  tha  correction  term. 
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For  r^(r)  indopondent  of  2,  the  eigenmode  of  Eq.  3.4.3  can  be 
determined.  This  allows  one  to  calculate  the  propagation  of 

e(x;y^  +  3.4.22 

n 

where 

A„»//  E(i;y^t(f„(x;)^d*dy.  3.4.23 

Using  Eq.  3.4.22  and  multiplication  with  the  correction  term,  defined  by  Eq. 
3.4.21,  alQWB  the  propagation  of  Hght  in  small  steps. 

Appicatian  of  the  beam  propagation  method  to  media  with  smaH 
refractive  iTKlex  variations  allows  one  to  analyze  beam  propagation  through 
integrated  waveguide  structuree.  For  this  case,  the  operator  ft  can  be 
erqweeeed  in  terms  of  Fourier  transforms  (see  Eqs.  3.4.5  and  3.4.6).  The 
propagation  step  defined  by  Eq.  3.4.22  can,  thus,  be  calculated  using  fast 
Fburier  transforms  (FFTs).  Performing  a  two  dirnensionalFFT  extends  the 
beam  propagation  rnethod  to  three  dirnensions,  allowing  one  to  analyze  the 
dMedrlc  waveguide  structure  Hiustrated  in  Fig.  3.1.1. 

The  beam  propagation  method  for  integrated  optical  waveguide 
appMcations  has  been  prograrnrned  into  a  SUN  Sparc  station  by  A.  Ticknor, 
of  Loefeheed  Missiles  and  Space  Company.^  Although  it  does  not  generate 
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Iwil  mmiMm  for  Ih6  initraction  length  of  directional  o(M^ 
dMaclrlewaregiiida  and  effective  index  methods  do,  the  beam  propagation 
method  doee  provide  an  exoeiient  prediction  of  device  performance.  It 
caiculiftee  the  ampMiidee  of  the  fields  as  they  propagate  through  the 
medium.  As  enexampie,  the  beam  propagation  method  computer  program 
acquired  from  Lockheed  was  used  to  generate  Fig.  3.4.1  which  gives  a 
visual  prediction  of  the  behavior  of  a  Hght  wave  propagating  through  an 
Aj^^Ga^iAe  waveguide  wkh  an  S4)end.  The  fbiiowing  parameters  were 
aubstltieed  into  the  beam  propagation  method  using  the  indices  from 
Wavegukis  Structure  1  (Chapter  II):  the  waveguide  width  a  =  3  pm,  the 
waveguide  core  thicknees  b«  1  pm,  the  operating  wavelength  X  »  827  nm, 
the  teffacdve  index  of  the  core  n, »  3.419,  the  refractive  index  of  the  top  and 
bottom  dad  regions  « rt4  s  3.389,  the  refractive  index  of  the  outside  dad 
regk)neri,srv»  1.000.  zero  propagation  loss  was  assumed  for  the  straight 

guidee.  A  step  sIzb  of  0.5  pm  was  used  to  calculate  the  propagating  field 
ampKudss.  The  8-bend  was  designed  using  the  formula”^ 

-  {hfl^z  -  (/i^)sin[(2]i/i)i3,  3.4.24 

where  /  is  the  longitudinai  separation  in  the  z  direction  aiid  h  is  the  lateral 
oAmI  in  the  xdkection.  Eq.  3.4.24  aHows  one  to  minimize  the  radiative  loss 
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and  maximize  the  mode  confinement.  For  the  S-bend  structure  illustrated 
in  Fig.  3.4.1,  /i  >  20  )im  and  /«  500  tun. 


CHAPTER  IV 


ELECTRO-OPTIC  PROPERTIES  OF 
ALUMINUM  GALLIUM  ARSENIDE 


A4Qai.xAs  is  an  optically  isotropic  crystal  belonging  to  the  cubic  43m 
symmetry  class  known  as  zincblerxle.^  Certain  Isotropic  crystals,  including 
Ai^Qa^^,  possess  electro-optic  (EO)  properties  in  which  a  polarization  arni, 
hence,  a  birefringence  can  be  induced  when  an  electric  field  (£)  is  applied 
to  the  crystal  (Pockets  eflect)."'^^  The  Pockels  effect  is  a  linear  EO  effect; 
the  induced  birefringence  is  directly  proportional  to  the  first  power  of  the 
applied  electric  field.  Thus,  Alj^Ga^^  exhibits  anisotropic  behavior. 

In  an  anisotropic  crystal  the  polarization  induced  by  E  is  dependent  on 
the  diiectlon  of  Therefore,  the  speed  at  which  light  travels  through 
an  anisotropic  crystal  is  dependent  on  the  direction  of  propagation.  A  light 
beam  propagating  through  an  anisotropic  crystal  along  the  optic  axis  (0)  is 
defined  as  an  ordinary  ray.^*^^  The  refractive  index  seen  by  this 
propagating  beam  is  defined  as  the  ordinary  refractive  irKtex  (n,|  or  rtj.  n„ 
is  independent  of  the  direction  of  the  polarization.”'’^^  A  light  beam 
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propagating  through  an  anisotropic  crystai  at  any  other  direction  to  O  is 
defined  as  an  extraordinary  ray/°*  The  refractive  index  seen  by  this 
propagating  beam  is  defined  as  the  extraordinary  refractive  index  (nj^  or  r>J. 
nj^  is  dependent  on  the  direction  of  poiarization Hence,  n^^  is  dependent 
on  the  direction  E.  Since  Ai^Qa^.^  has  one  optic  axis  it  behaves  iike  a 
specific  group  of  anisotropic  crystais  known  as  uniaxial  crystais/°°'^°^  If  we 
define  the  z  axis  as  the  optic  axis  then  n^  =  n,,  and  n^  =  n^  =  n^^.  The  index 
surfaces  can,  therefore,  be  represented  by  a  sphere  with  radius  arKi  an 
eiHpsoid  with  radius  n„  parailei  to  O  and  radius  nj^  normal  to  the  O  (see  Fig 
4.0.1).^**  For  nj,  >  n|,  the  uniaxial  crystai  is  known  as  a  positive  type  and 
for  nj.  <  nil  the  uniaxial  crystal  is  known  as  a  negative  type/^  As  can  be 
seen,  the  sphere  and  eliipse  are  tangent  to  each  other  only  at  the  points 
along  the  optic  axis.  The  amount  of  induced  birefringence  is  defined  as  the 
difference  between  these  two  refractive  indices  and,  as  stated  above,  is 
linearly  proportional  to  Pockels  effect  is  given  mathematicaily 

n,|  -  nj^  as  mE,  4.0.1 

where  m  is  a  proportionaiity  constant.  Modulating  the  electric  field  aibws 
one  to  modulate  the  input  signal  and  provide  the  basis  for  building  an 
etectro-optic  switch. 
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In  this  chapter  we  will  determine  the  induced  birefringence  and  derive 
an  expression  for  electro-optic  modulation  using  the  EO  properties  of  the 
specific  AtyGa^xAs  structure  to  be  fabricated  for  this  experiment. 

Linear  Electro-Ootic  Effect 

We  will  begin  our  analysis  with  the  constitutive  equation”^’^ 

0=e^E+P.  4.1.1 

In  Eq.  4.1 .1 ,  D  is  the  electric  displacement  or  electric  flux  density  vector,  E 
is  the  electric  field  vector,  P  is  the  electric  polarization  vector  and  Sq  is  the 
permittivity  of  vacuum.  In  many  isotropic  materials  the  induced  polarization 
is  directly  proportional  to  the  applied  field  strength,  except  for  the  case  of 
very  high  fields.  We  can  write:^^ 


4.1.2 


where  %  is  the  dielectric  susceptibility  of  the  medium.  Combining  Eqs.  4.1 .1 
and  4.1.2  yields^®°’^°® 

D=eE,  4.1.3 
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where  e  »  £o{1  x)  the  permittivity  of  the  medium. 

The  analysis  is  similar  in  anisotropic  materials.  Eq.  4.1.2  can  be 
expressed  in  its  tensor  form  as^^''°° 

Xl2^ 

Py  *  fioOCai^x  +  XaPy  +  Xas^i)  4.1 .5 

and 

Pz  ”  ^(Xai^  ^  XaaPy  ^  Xsa^j^*  4.1 .6 

where  X|  ere  the  coefficients  of  a  3x3  array  called  the  dielectric  susceptibility 
tensor  of  the  medium  for  which  /s  1, 2  or  3  and  /s  i,  2  or  3.  The  suffixes 
/  and  /  designate  which  axes  the  coefficients  correspond  to  since  1 
represents  x;  2  represents  y  and  3  represents  z.  Choosing  x,  y  and  z  so 
that  the  diagonal  terms  vanish  gives  the  principal  dielectric  axes  of  the 
crystal:^® 

4.1.7 

4' 

Py  ~  ^rXaz^y  4.1 .8 


and 
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4.1.9 


Combining  Eq.  4.1 .1  with  Eqs.  4.1 .7  -  4.1 .9  yieids’*’’®°’’“ 

4.1.10 

Dy  ®  4.1  .1  1 

and 

D,  =  C33£,  4.1.12 


where 

®11  ®  ^0  +  Xli)» 

£22  ~  X22)  4.1.14 

and 

^83  “  CoO  +  X33) 


are  the  diagonal  coefficients  of  a  3x3  array  which  is  referred  to  as  the 
permittivity  tensor  of  the  medium. 

If  we  introduce  the  phase  velocity  (u  =  where  p  is  the 

permeabiiity  of  the  medium,  it  becomes  evident  that  for  anisotropic  crystals, 


x>  depends  on  the  direction  of  the  appiied  fieid.  Thus  we  have 


4.1.16 


^y=(M«22)‘ 


4.1.17 


4.1.18 


For  an  anisotropic  crystal,  at  least  one  of  the  permittivity  coefficients  is  not 
equal  to  the  other  two  (i.e.,  ^  £22  or  £33).  Only  the  polarization  parallel  to 

the  appHed  electric  field  and,  hence,  the  permittivity  parallel  to  E,  will 
contribute  to  the  phase  velocity.  For  example,  if  an  input  beam  is 
propagating  along  the  2  axis  and  E  is  appiied  parallel  to  the  y  axis,  E22  will 
be  the  permittivity,  a  polarization  (P^  will  be  induced  and  the  phase  velocity 
wiH  be  Vy  (see  Eq.  4.1.17). 

The  index  of  refraction  is  defined  as  the  reciprocal  of  the  relative  phase 
velocily.*“’*®'^“ 


n  =  c/o. 


4.1.19 


where  c  is  the  speed  of  light  in  vacuum.  Going  to  tensor  form,  we  can  write 
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the  principal  refractive  indices 

n,=:cA)^=(MAii)^^  4.1^0 

n^  =  c/Oy  =  4.1 51 

and 

n,  =  C/Oy  a  (PAas)^®  4.1 52 

In  Eqs.  4.1 51  -  4.1 .22,  p,  a  ^/^  is  the  relative  permeability  of  the  medium, 
Pe  isthe  permeabitity  of  vacuum  and  e^a  c/Eo  are  the  coefficients  of  a  3x3 
anay  called  the  relative  permittivity  tensor  of  the  medium.  For 
nonconductive  media  the  relative  permeability  (p^)  can  be  assumed  unit/* 
and  the  principai  refractive  indices  then  become^*^ 


n,a 


4.153 


ny«(c22/e«)'®  4.154 

and 

4.155 


Using  the  definition  for  the  refractive  index  along  with  the  analysis  for 
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anisotropic  crystals  given  above,  it  becomes  evident  that  one  can  induce  a 
change  in  the  refractive  index  with  an  applied  electric  field.  This  change  will, 
in  turn,  induce  birefringence,  a  phenomenon  known  as  the  Pockeis  or  linear 
electro-optic  effect. 

Index  Ellipsoid 

Because  AlxGai.xAs  has  been  shown  to  exhibit  Pockeis  effect,^'^  we 
can  use  the  anisotropic  analysis  to  determine  the  induced  birefringence.  We 
wHi  start  with  the  general  index  ellipsoid  for  a  crystal  in  its  simplest  form. 
This  wHi  allow  us  to  determine  the  direction  of  polarization  as  well  as  the 
corresponding  refractive  indices  of  the  crystal.  The  index  ellipsoid  is  given 
in  the  principal  coordinate  system 

+  /hy  +  =  1 .  4.2.1 

In  Eq.  4.2.1,  x,  y  and  z  are  the  principal  axes  (the  direction  in  the  crystal 
along  which  the  electric  field  vector  O  and  the  electric  displacement  vector 
(13)  are  paralleO  arxl  suxi  1/n/  are  the  principal  refractive  irKlices 

along  their  respective  axes.  The  lengths  of  the  major  axes  of  the  index 
ellipsoid,  represented  by  Eq.  4.2.1 ,  are  2n^  2nyand  20^,  parallel  to  directions 
x;  y  and  z,  respectively.^*^  The  index  ellipsoid  is  used  to  find  the  two 
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rafractive  indices  and  two  oorrasponding  directions  of  Dassodated  with  the 
two  independent  plane  waves  propagating  along  an  arbitrary  direction  s  in 
the  crystal.  These  indices  are  determined  using  Eq.  4.2.1 .  Visually  depicted 
in  Fig.  4.2.1,  the  major  and  minor  axes  of  the  ellipsoid  are  in  the  directions 
Df  and  D2,  respectively,  with  the  length  of  equal  to  20^  and  the  length  of 
equal  to  2n2. 

Next  we  wili  introduce  the  optical  dielectric  impermeability  tensor  of  a 
crystal.  The  coefficients  of  this  tensor  deperKi  on  the  distiibution  of 
bond  charges  in  the  materiai.'*'^*’”  qi  are  found  by  taking  the  redprocai  of 
the  relaliveperrnittivity  or  dielectric  constant.^  has  been  defined  in  terms 
of  the  refractive  index  of  the  crystal  as^ 

q|se^f>1/h/  4.2.2 

(see  Eqs.  4.1 .23  •  4.1 .25).  Application  of  an  electric  field  £  to  the  crystal 
resulls  in  a  redistrfoution  of  the  bond  charges  of  the  crystal  and  a  possible 
defomiation  of  the  ion  lattice.  This  causes  a  change  in  the  optical 
impermeability  tensor  and.  hence,  a  change  the  refractive  index  of  the  ^ 
crystal.  This  allows  us  to  write  the  index  ellipsoid  in  the  presence  of  an 
electric  field:^ 
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{\/n^  +  +  (1M,V  +  A^Ji^  -fA*/  +  +  A^yz 

•»>  A«2r4- A^-i>  A^-f  Var-i>  Aaiors  1.  4.2.3 

Fbr  E  s  0,  A,  thfXHjgh  A,  s  0,  and  Eq.  4.2.3  reduces  to  Eq.  4.2.1. 
Producing  the  notations 

j^»xr«(11)-4l 

y->y-(22)-»2 

y>2Z>(33)-»3 
>zx  (23)  -4  4 

Z)r«(31)-4  5  4.2.4 

jy  *  (12)  -♦  6 
jya  (32) -4  7 
X2x(13)-4  8 
ja  a  (21) -4  9, 
along  with 

(lAnP),  X  1/h/  X  A, 

(1/ii%  X  1/h/  X  Aa 
(lAi*),  X  1/n/  X  Ag 
(1/in?)4  X  1/n^,  X  A* 

(1/n*)5xl/hA*As 


4.2.5 
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(1/in*)e*  IMj/lyS  Ae 

(1/h*)7  =  1/n/iy  =  A, 

(1/n*)*  =  1/n/>,  s  Ae 
(l/n*),  s  l/n/t,  =  Ag, 


Eq.  4^.3  becomes 

[1/n,*  +  (1/n*),]A*  +  [1/n/  +  (l/h*)*!/  +  [1/n/  +  (l/n^y/ 

+  C(1/n*)4lyz+  [(1/n^sl«+  {(1/n*)«]xK+ 

+  [(1/n*)8lxz  +  [(1/n*)9lyx  =  1 .  4^.6 

It  has  been  shown  that  is  symmetric,  provided  that  the  medium  is 
loseiees  and  optically  active/*  This  allows  the  permutation  of  /  and  y,  and 
the  notation  presented  in  Eq.  4.2.5  and  Eq.  4.2.6  reduces  to 


4.2.7 

I 

I 

i 


;^«xxs  (11)  1 

/sXy=(22)-^2 
i?a=  2?=  (33)  ->3 
yz*  zy^  (23)  =  (32)  -♦  4 
2>f«xz*(31)  =  (13)->5 
xys  yx=  (12)  =  (21)  6, 


=  1/n/ 

(iyh*)2  ss  1/n/ 

(1/h*)3  =  1/n/ 

42.B 

s  l/n^n^  =  1/h/iy 
(1/r^s  =  =  1/n/i, 

(1/n*)3  a  l/Oj/iy  s  1/n/i*. 

Substituting  Eq.  4^.7  and  Eq.  4^.8  into  Eq.  4.2.6  yields'*^ 

[lAV*  +  (l/n^ilV*  +  [1/n/  +  (l/n*)^)^  +  [1/n/  + 

+  2{(1/n*)4]yz  +  2[(1/n*)8lxz  +  2[(1/n*)el;y  =  1 .  4.2.9 


The  electro-optic  (EO)  coefficients  are  traditionaily  defined  as^^^^ 

■  il#(0) »  All,  =  4.3.1 

in  Eq.  4.3.1,  E  is  the  applied  electric  field,  are  the  coefficients  of  a  3x9 
array  csdied  the  linear  electro-optic  (Pockets)  tensor  of  the  medium  and 
are  the  coefficients  of  a  6x6  array  called  the  quadratic  electro-optic  (Kerr) 
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tensor  of  the  medium  for  which  /r  s  1 ,  2.  or  3  and  /  s  1 ,  2,  or  3.  In  the 
equation  above,  we  have  omitted  terms  higher  than  the  quadratic.  These 
higher  order  effects  are  so  weak  compared  with  the  linear  and  quadratic 
effects,  that  they  can  be  neglected  for  most  apptications.^^’^'^  For 
centrosymmetric  crystals  (crystals  possessing  inversion  symmetry  for  all 
physicai  properties'"),  the  linear  EO  effect  vanishes  (i.e.,  -  0).  This  is 

caused  by  the  spatial  inversion  symmetry  of  such  crystals.  Therefore,  the 
quadratic  EO  effect  will  dominate.'^  For  noncentrosymmetric  crystals 
(crystals  that  do  not  possess  inversion  symmetry  for  all  physical 
properties'"),  both  the  linear  and  quadratic  effects  will  be  present.'^  From 
the  section  above,  the  EO  effects  will  depend  on  the  ratio  of  the  applied 
electric  field  (£)  to  the  intra-atomic  electric  field  (E^  binding  the  charged 
particles.  This  ratio  can  be  estimated  by  finding  the  magnitude  of  at  a 
position  (7)  arising  from  a  point  charge  (q)  located  at  a  distance  (u)  from  T. 
E„  at  7  is  given  by*^ 

E„=(1/4jteo)c^r/.  4.3.2 

The  minimum  field  inside  the  atom  can  be  determined  by  substituting  the 
lattice  constant  ao  of  the  crystal  for  u  and  the  elementary  electron  charge  for 
q.  In  this  experiment,  AioaoGaoToAs  was  used  for  the  core.  The  lattice 
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constant  for  is  given  by  Eq.  2.3.30.  For  x  =  0.300,  - 

8.8541 8782  xIO'^F/m  and  e*  1.6021892x10*’®  C,’®®  foe  minimum  electric 
field  inside  foe  atom  is  found  to  be  »  4.5  x  10^  V/cm.  This  is  expected 
to  be  far  greater  than  foe  applied  field.  As  a  result,  the  quadratic  EO  effect 
should  be  small  erxHjgh,  when  compared  to  foe  linear  EO  effect,  to  allow  us 
to  neglect  foe  terms.  Eq.  4.3.1  becomes 

-  ^|(0)  -  Aq,  =  /pE*.  4.3.3 

The  permutation  property  along  with  Eqs.  4.2.2,  4.2.7,  4.2.8  and  4.3.3 
yields®^’®* 

3 

A(1/h*),«JEr*E*,  4.3.4 

where,  for  example,  A(1/n®)4  =  +  r^E^  +  r^E^.  Substitution  of  Eq.  4.3.4 

into  Eq.  4.2.9  gives  the  general  index  ellipsoid  in  foe  presence  of  an  electric 
field  in  terms  of 


[1V  +  r,^  +  [1M;  +  +  [1M/  + 

*  2[rtt£J>z  +  + AryS/y  =  i  • 


4.3.5 
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in  oitier  to  determine  the  EO  tensor  for  crystals  belonging  to  the  cubic 
43m  class  we  will  start  with  the  EO  tensor  for  the  tridinic  1  crystal  group  in 
matrix  form:^“ 


'’ll 

^12 

^3 

^21 

^22 

^23 

^31 

^32 

^33 

r4i 

^43 

^S1 

^S2 

^63 

^61 

^62 

^63 

Neglecting  the  permutation  property  for  the  moment,  Eq.  4.3.6  expands  to 


^11 

^112 

^13 

^999 

^223 

^331 

^332 

^333 

^232 

^321 

^322 

^323 

^131 

^132 

^133 

^311 

^312 

^313 

^121 

^122 

^123 

^11 

^212 

^213 
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Since  zincblende  crystals  have  axes  of  fourfold  symmetry  along  the  cube 
edges,  we  can  let  z  be  the  axis  of  two  fold  rotation  which  can  be 


represented  by  the  matrix'^ 

-10  0 

a,  =  0  -1  0 

0  0  1 


Eq.  4.3.8  dictates  that  the  suffixes  of  the  EO  coefficients  change  as  follows: 
1  -1 , 2  ->  -2  and  3  3.’“  Substituting  this  property  into  Eq.  4.3.7  yields 


“^111 

”^112 

fli3 

“^221 

“^222 

^223 

“^331 

"^332 

^333 

^231 

^232 

~^Z33 

^321 

^322 

“^323 

^131 

^132 

"^133 

^311 

^312 

“^313 

“^121 

“^122 

^123 

*■^211 

"^212 

^213 

where 

^111  "  “^111* 
=  ^123  ^(-1X-2X3)  ”  ^123> 
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^  ^(-2K-2X3)  * 

etc.  Setting  s  ±r^  leads  to  the  following  condition:  only  if  = 

0.^®  Equating  Eqs.  4.3.7  and  4.3.9  with  this  condition  yields 

0  0 

0  0 

0  0 

^231  ^232 

^321  ^322 

^131  ^132 

^311  ^312 

0  0 

0  0 

Shce  cubic  43m  class  crystals  abo  have  axes  of  fourfold  symmetry,  which 
includes  three  mutually  perpendicular  twofold  axes  x,  y  and  z,  this  allows  us 
to  perform  a  triple  product: 

iiixl^  =  =  Jikxtj.  4.3.1 1 

Performing  this  triple  product  with  the  coefficients  in  Eq.  4.3.10  yields 
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= 


4.3.12 


where 

Ofc  ®  '113  =  0, 

fifc  =  Taa^  *  -fgai  =  0, 

V=  ^231  “>>'(^^=^231. 

etc.  Finally,  we  need  to  take  into  account  the  axes  of  threefold  symmetry. 
The  simplest  threefold  symmetry  operation  is  x  -» y,  y->  z  and  z  x.  This 
is  represented  by  the  matrix^^ 


=  0 


4.3.13 
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Eq.  4.3.13  dictates  that  the  suffixes  of  the  EO  coefficients  change  as  foilows: 
1  2, 2  3  and  3  ->  1  Using  the  threefold  symmetry  operation  gives 


~  ^123  ^(2X3X1)  "  ^231 


4.3.14 


and 


®  ^(3X1X2)  ®  ^312* 


4.3.15 


This  proves  that 


^123  ®  ^231  ~  ^312*  4.2.16 

Going  back  to  contracted  notation,  the  electro-optic  tensor  (rj^  for 
crystals  belonging  to  the  cubic  43m  group  can  be  represented  by  the  well 
known  matrix^* 

0  0  0 

0  0  0 

0  0  0 

r,,  0  0 

0  0 
0  0  r,, 


Substituting  the  EO  tensor  for  zincblende  crystals  back  into  the  index 


ellipsoid  (Eq.  4.3.5)  gives^^ 

[1/n/];^  +  [1/n/]/  +  [1/n/]^ 

+  +  2xzr^^Ey  +  2xyr^^Eg  =  1 .  4.2.1 8 

Taking  into  account  that  cubic  43m  crystals  are  opticaily  isotropic  allows  us 
to  further  simplify  the  index  ellipsoid.  This  yields^ 


(1/n*)[>^  +  /  +  ^  +  2yzr^^E^  +  2xzr^yEy  +  2xyr^^Eg  *  1 .  4.2.19 


Induced  Refractive  Indices  and  Birefringence 
We  need  to  determine  the  direction  and  magnitude  of  the  principal  axes 
of  the  new  ellipsoid  represented  by  Eq.  4.3.19  since  they  may  not  coincide 
with  X,  y,  z.  For  this  experiment,  we  are  fabricating  an  EO  directional 
coupler  switch.  This  dictates  that  the  electrodes  be  placed  parallel  to  the 
waveguides  along  the  optic  axis.  Therefore,  propagation  of  the  input  beam 
will  be  along  the  <01 1>  (i.e.,  z)  direction  through  the  directional  coupler  of 
length  L,  and  the  electrodes  will  be  placed  on  the  top  and  bottom  of  the 
coupler  running  along  z  and  separated  by  a  distance  of.  The  electric  field  E 
will  be  applied  along  the  <100>  (i.e.,  y)  direction.  With  E  applied  normal  to 
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the  direction  of  the  input  beam,  we  have  what  is  defined  as  transverse 
modulation.^*  A  transverse  modulator  using  a  zincbiende  crystal  in  the 
orientation  used  for  this  experiment  is  illustrated  in  Fig.  4.4.1 .  Since  the 
electric  field  is  in  the  y  direction,  we  have 

and 

Ey=E.  4.4.1 

This  is  similar  to  one  of  the  three  cases  presented  by  Namba**  with  E  along 
the  <001  >  direction  and  the  input  beam  propagating  along  the  <100> 
direction.  These  were  defined  by  Namba  as  the  x  arxl  y  directions, 
respectively.  Scribe  arxi  break  experiments  performed  on  GaAs  wafers  have 
derrKK^strated,  however,  that  the  preferential  cleave  planes  for  GaAs  are 
(011),  (oTi),  (oiT),  (oTT),  (101),  (loT),  (Toi),  (ToT),  (iio),  (iTo),  (Tio)  and 
(TTO)  as  shown  in  Fig.  4.4.2.^**  This  suggests  that  even  though  Namba’s 
case  can  be  solved  theoretically,  it  cannot  practically  be  demonstrated.  This 
is  the  argument  used  for  the  crystal  orientation  chosen  for  this  experiment. 
Substituting  Eq.  4.4.1  into  Eq.  4.3.19  yields  the  index  ellipsoid 


(1/n^[x*  +  /  +  2^  +  2xzr„E^  1. 


4.4.2 


(010) 


(Oil) 


(001) 


Figure  4.42 

Representation  of  a  GaAs  Crystal 
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In  order  to  detonnine  the  direction  and  magnitude  of  the  principal  axes  of  the 
new  eMpeoid,  we  first  need  to  convert  Eq.  4.4.2  into  its  diagonal  form  by 
transforming  the  coordinate  system  (x.y.2)  into  a  new  coordinate  system 
(x^y y).  This  is  done  be  rotating  the  axes  perpendicular  to  the  applied 
electric  field  by  45^^'***  For  this  experiment  we  will  rotate  the  xand  zaxes, 
which  gives 

xs  x'C08(45®)  -  Z'sin(45*),  4.4.3 

y*  /  4.4.4 

and 

z  =  x'8ln(45®)  +  z'cos(45®).  4.4.5 

Substituting  Eqs.  4.4.3, 4.4.4  and  4.4.5  into  Eq.  4.4.2  yields 

[(1/n*)  +  +  [1/nV  +  K1/n*)  -  r^iEJz**  =  1 .  4.4.6 

Using  Eq.  4.3.4,  we  have 


[(1/n*)  +  A(1/h*)lx'®  +  [1/n*]y'*  +  [(1/n^  -  A{Vn^)]z'^  =  1 . 


4.4.7 


From  Eq.  4.4.7  we  can  write 


=  1/n*  +  A(1/n*),  4.4.8 

1/(n,^)*  *  1/n*  4.4.9 

and 

M{T\rf  =  1/n*  -  A(1/n^.  4.4.10 

Using  Eq.  4.0.1  for  Pockels  effect  along  with  Eqs.  4.2.2,  4.3.3  and  4.3.4 
yields  the  induced  change  in  the  refractive  index 

An  =  mV4,E’  =  m  A(1/n*),  4.4.1 1 

where  m'  is  a  proportionality  factor  we  will  use  to  equate  An  and  A(1/n^.  If 
we  assume  r4i£ <  1/n^,  (i.e.,  A(1/n^  sufficiently  small),  this  allows  us  to  use 
the  reiation^* 

A(1/n*) »  3(1/n*)  =  -(2/n®)3n.  4.4.12 


Multiplying  both  sides  by  -(n^/2)  yields 


an  -  -(n*/2)a<1/n^. 


4.4.13 


Equating  Eqs.  4.2.11  and  4.2.13  leads 

An  «  m'  -  -(n®/2)A(1/n").  4.4.14 

Using  Eq.  4.3.4  once  again,  we  can  write 

An  -  -<n*/2)r4^£  4.4.1 5 

l\w  leads  to  the  new  indices^^ 

ny  =  n  +  An  =  n-  {r\^l2)r^^E,  4.4.16 

n^r  =  n  4.4.17 

and 

n^  =  n-  An  =  n  +  {T\^f2)r^^E.  4.4.18 


The  induced  birefringence  for  the  new  axes  with  the  input  beam  propagating 
perpendicular  to  the  applied  electric  field  (i.e.,  for  transverse  modulation)  can 
now  be  given  by^®'®^ 


=  (n^l2)r^^E. 


4.4.19 


Eq.  4.4.19  is  the  definition  for  the  imJuced  birefringence  that  we  wiil  use  to 
describe  the  phase  retardation  experienced  by  the  input  beam  as  it 
propagates  through  the  crystal. 

Electro<3otic  Modulation 

For  our  experiment  we  will  consider  an  optical  field  that  is  incident  and 
normal  to  the  x'y'  plane  propagating  along  the  zaxis  with  E parallel  to  the 
y  axis.  At  the  input  plane  z  s  o  the  optical  field  can  be  resolved  into  two 
mutually  orthogonal  components  polarized  along  x'  and  y'.  The  x' 
component  propagates  as^^^" 


Ex'  s  AexpCKcot  -  (a]/c)nx'^],  4.5.1 

where  o)  is  the  angular  frequency,  t  is  time  and  A  is  a  constant.  The  / 
component  propagates  as^°^^^ 

Ey  s  Aexp[/(ost  -  ((iyc)ny'Z)].  4.5.2 

The  difference  in  phase  between  these  two  components  at  the  output  plane 
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Zb  L  i8  caUed  the  phase  retardation.^°^^**  It  is  given  by  the  difference  of  the 
e)0onential  terms  in  Eqs.  4.5.1  and  4.5.2.  Performing  the  subtraction  yields 

♦ 

=  («j/c)[nyr  -  ny]L,  4.5.3 

where  [n^r  -  n^]  is  the  induced  birefringence  found  in  the  previous  section. 
Substituting  for  the  induced  birefringence  using  Eq.  4.4.19  gives^°^^°^ 

=  ((a/c)[(n®/2)f4,^L  4.5.4 

Using  the  relation  tdo  s  2ftfK  with  X  defiried  as  the  wavelength,  the  phase 
retardation  becomes^^”* 

B  VL)l{Xd^,  4,5.5 

where  V^is  the  voltage  applied  between  electrodes  separated  by  a  distarrce 
d.  For  this  experiment,  we  want  a  phase  retardation  of  180**  or  =  ic. 

The  voltage  necessary  to  realize  a  n  phase  change  is  defined  as  half  wave  » 
voitage.^°^’^°^  Setting  the  left  side  of  Eq.  4.5.5  equal  to  n  gives  us 


V,  = 


4.5.6 
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14  is  the  voltage  necessary  to  realize  half  wave  potarization  in  zincblerKle 
crystals.  Eq.  4.5.6  is  also  the  same  relation  that  was  given  for  complete 
switching  using  a  zero-gap  directional  coupler  (see  Eq.  1.5.7).  For 
conventional  directional  coupler  switches  the  voltage  required  for  switching 
the  input  beam  completely  from  one  channel  to  the  other,  requires  a  phase 
change  of  VSic.  SubstHutir^  into  Eq.  4.5.5  yields 

=  (V3Xd)/(nV4,L).  4.5.7 

Eq.  4.5.7  is  the  same  relation  that  was  given  for  complete  switching  using 
a  conventional  directional  coupler  (see  Eq.  1 .4.1 1 ).  Since  we  are  fabricating 
zero-gap  directional  couplers,  Eq.  4.5.6  is  the  relation  we  will  use  for  the 
switching  voltage  in  this  experiment. 

Electro-optic  modulators  and  directional  couplers  using  QaAs  have  been 
modelled,  fabricated  and/or  analyzed  extensively  over  the  years  by  many 
researchers. The  EO  coefficient  for  GaAs  has 
been  measured  or  modelled  by  several  researchers  as  weil.^'^'^^^^ 
Conducting  an  exhaustive  search  of  the  open  literature  for  AI^Gai.,^ 
waveguide  modulators  and  directional  couplers,  however,  has  turned  up  only 
a  few  papers.®’^*"^^  For  the  EO  coefficient  for  AlxGai.xAs,  only  one  paper 
found  ttiat  gives  a  measured  value  of  This  was  for  a  specific  Ai 


160 

concentration  (x  s  o.17)  and  at  a  specific  wavelength  (X  =  1 .1523  pm).  In 
addition,  no  models  for  the  determining  the  EO  coefficient  for 

were  discovered.  Therefore,  for  the  Al^Ga^xAs  composition  used  in  this 
experiment  will  have  to  be  estimated.  This  was  also  the  case  for  those 
previous  works  on  Ai^Gaf.^  modulators.®’^®®"’'*^  What  is  different  from  those 
earlier  papers,  however,  is  the  fact  that  we  have  a  measured  value  of  for 
Aio.^7Gao.g3As.^^  This  will  allow  us  to  model  for  other  concentrations  of 
AixGa,.xAs  and  should  provide  a  closer  estimate. 


CHAPTER  V 


EXPERIMENTAL  PROCEDURES  AND  DEVICE  CHARACTERIZATION 

In  the  preceding  chapters,  the  theory  behind  optical  Fredkin  gates  arni 
zero-gap  directional  couplers  was  presented.  In  this  chapter,  the  actual 
experimental  procedures  used  to  fabricate  and  evaluate  these  devices  will 
be  detailed.  AlxGa^.xAs  material  growth  and  characterization  as  well  as 
waveguide  and  device  fabrication  and  performance  will  be  addressed. 
Optical  digital  logic  will  also  be  demonstrated.  In  addition,  a  method  to 
determine  the  linear  electro-optic  coefficient  for  AI^Ga^.^  will  be  developed. 


Material  Growth 

The  AI^Ga^^  structure  used  to  fabricate  Fredkin  gates  for  this 
experiment  was  grown  using  both  molecular  beam  epitaxy  (MBE)  and  metal- 
organic  chemical  vapor  deposition  (MOCVD). 

MBE  is  an  epitaxial  growth  technique  that  utilizes  the  reaction  of  one  or 
more  thermal  beams  of  atoms  or  molecules  that  have  a  crystalline  surface 
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under  uRra-high  vacuum  (i.e..  <  1 0*'*’  Torr).  Molecular  beams  are  generated 
from  thermal  effusion  cells  known  as  Knudsen  sources/^  The  molecules 
interact  on  a  heated  crystalline  substrate  producing  a  single-crystal  layer. 
Each  of  the  constituent  elements  and  dopants  required  for  making  up  the 
composition  of  the  grown  film  is  emitted  from  its  own  Knudsen  source.  Fig. 
5.1 .1  is  a  schematic  diagram  of  the  MBE  process  for  growth. 

Each  source  is  independently  temperature  controlled  and  arranged  around 
the  substrate  in  such  a  way  as  to  ensure  that  the  uniformity  of  both  film 
composition  arKi  thickness  is  optimized.^^  The  addition  of  mechanical 
shutters  for  each  source  as  well  as  rotating  the  substrate,  further  optimizes 
the  growth.^*^^  Growth  using  MBE  enhances  the  electrical  and  optical 
properties  of  the  grown  fiims.^^  The  growth  rate  for  MBE  is  typically  1 
pm/hour.^'^^  Growth  temperatures  using  MBE  range  from  550  -  650'’C  for 
Ai^Ga^^  deposition.®^’^^ 

MOCVD  growth  of  lll-V  materials  utilizes  a  liquid  source  of  the  group  III 
components  (Al  and  Ga)  that  are  vaporized  in  a  temperature  controlled 
stainless  steel  bubbler.^*^'^®®  The  group  III  vapor  is  transported  to  a  process 
tube  by  a  hydrogen  (Hy  carrier  gas  and  pyrolyzed  with  the  group  V 
component  (As).^  This  forms  an  epitaxial  layer  on  the  substrate.  Rg.  5.1 .2 
is  a  schematic  diagram  of  an  MOCVD  reactor  for  AlxGai.xAs  growth.  Growth 
rates  using  MOCVD  are  controlled  by  adjusting  the  Hs  flow  through  the 
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Rgure5.1^'^ 

MOCVD  Growth  Process  for  AI^Qai^ 
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bubbler  and  typically  range  from  1-10  pm/hour.'^*^"  Film  layers  as  thin  as 
25  A  can  be  deposited.^^'^”  Growth  temperatures  using  MOCVD  range  from 
550  -  750°C  for  Al.Gai.xAs  deposition/^''^  P-type  dopants  are  introduced 
using  zinc  (Zn)  and  n-type  dopants  are  introduced  using  tellurium  (Te).^^ 

The  AI.Gai..As  waveguide  structure  showing  the  target  Al  concentrations 
and  layer  thicknesses  is  itiustrated  in  Fig.  5.1 .3.  The  target  thickness  of  the 
Alo^GaoToAs  core  region  grown  was  1  pm  in  order  to  provide  single  mode 
operation  at  827  nm  (see  Chaffer  II).  The  target  thickness  for  the  top 
Aio^GSo.66As  dad  layer  grown  was  1  pm  in  order  to  provide  low  modal 
attenuation.^  The  target  thickness  of  the  bottom  A^^GSoesAs  clad  layer 

grown  ranged  from  2  •  4  pm  in  order  to  provide  isolation  between  the  device 
arxi  the  substrate.  It  has  been  shown  that  a  thicker  bottom  clad  lowers  the 

propagation  loss.^  The  altemating  Alo^Gao.72As  and  GaAs  super  lattice 
buffer  layers  grown  between  the  GaAs  substrate  and  the  bottom  dad  layer 
have  been  shown  to  provide  improved  surface  morphology  of  the 
subsequently  grown  layers.®’®*’^® 

The  MBE  chambers  used  for  the  Al.Gai..As  growth  were  a  Varian  360, 
a  Vartan  GEN  1.5  and  a  Vartan  GEN  2.  All  three  MBE’s  had  40  cc  oven 
sizes.  They  were  capable  of  processing  a  single  2*  wafer.  An  Al 
concentration  uniformity  of  6.5%  had  been  measured  over  a  1  *  radius  for  ail 
three  machines.®  Samples  562  and  563  were  grown  on  Nippon  Mining  2* 
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Figure  5.1 .3 

Target  AixGai.xAs  Structure 
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GaAs  substrates  using  the  GEN  1 .5.  Samples  571  and  572  were  grown  on 
Outokumpu  Sernttronic  2*  GaAs  substrates  using  the  GEN  1 .5.  Sample  179 
was  grown  on  a  Spectrum  Technology  2*  GaAs  substrate  using  the  Varian 
360.  Samples  1108  and  1109  were  grown  on  Sumitomo  Electric  2*  GaAs 
substrates  using  the  GEN  2.  The  substrates  for  samples  562, 563, 571 , 572 
and  179  were  liquid  encapsulated  (LEG)  type  semi-insulating  with  no 
dopants.  The  substrates  for  samples  1108  and  1111  were  LEG  type  semi- 
insuiating  doped  with  indium  (In).  The  substrates  for  samples  562, 563, 571 

and  572  had  (100)  orientation  tilted  6°  off  toward  the  nearest  <110> 
direction.  The  substrates  for  samples  1108,  1111  and  179  had  (100) 

orientation  tilted  2^  off  toward  the  nearest  <110>  direction.  The  etch  pit 
density  of  the  substrates  for  samples  1108  and  1111  was  ^  1  x  lO’ 

defects/cnf .  The  etch  pit  density  of  the  substrates  for  samples  562,  563, 
571,  572  and  179  was  ^  1  x  1(f  defects/cm^  The  substrates  were  all 

polished  on  both  sides.  The  thickness  of  ail  substrates  was  450  pm  ±  10 
pm.  The  target  thickness  for  the  bottom  dad  grown  was  2  pm  for  samples 
562,563, 571, 572, 1106  and  1111,  and  3  pm  for  179.  Table  5.1.1  lists  the 
MBE  growth  parameters.  During  growth,  the  substrates  were  rotated  at  15 

rpm.  The  Arsenic  source  used  was  cracked  As  at  a  pressure  of  6.5  x  10^ 
Torr.  The  samples  were  grown  on  top  of  the  (100)  plane  for  ail 

substralee. 
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Number 


Table  5.1.1 

Parametere  for  AlxGa^^  Structure  using 
Varian  360,  GEN  1 .5  and  GEN  2  MBE 


Substrate 

Temperature 


Growth  Rate 
(monolayers/sec) 


Al 


0.30 


1108 

617 

0.70 

0.55 

0.30 

1111 

593 

0.70 

0.55 

0.30 

179 

630 

0.49 

0.39 

0.21 

The  MOCVD  chamber  used  for  the  growth  was  an  AixtronAiX  200.  The 
AIX  200  was  capable  of  processing  a  single  3*  wafer.  An  Ai  concentration 
uniformly  of  0.5%  has  been  measured  over  a  3*  diameter  using  the  AIX 
200.^*”  Samplee  121  and  128  were  grown  on  Mitsubishi  Monsanto  Kasei  3* 
GaAs  aubetralee.  Both  substrates  were  LEC  type  semi>insulating  with  no 
dopants.  The  substrates  for  samples  121  and  128  had  (100)  orientation. 
The  etch  pit  density  of  both  the  substrates  was  ^1x10^  defects/cm^.  The 
substrates  were  polished  on  the  (1 00)  plane  (front  side)  only.  Thebacksides 
were  lapped  and  briefly  etched.  The  thickness  of  both  substrates  was  600 
|im  ±20  pm.  The  ALCaa^^  samples  were  grown  on  top  of  the  (100)  plane 
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for  both  substrates.  In  addition  to  the  structure  shown  in  Fig.  5.1 .3.  a  0.1 
layer  of  QaAs  was  grown  between  the  core  and  each  clad  layer  of  samples 

121  and  128  to  provide  better  lattice  matching  between  the  various  Ai^Ga^. 
xAs  layers.^"  Table  5.1.2  lists  the  MOCVD  growth  parameters. 


Table  5.1.2 

Growth  Parameters  for  AixGai.xAs  Structure  using 
the  AIX  200  MOCVD 

Substrate  Flow  Rate 

Sample  Temperature  (sccm/sec) 


Number 

(“C) 

^*030 

Ga 

AsHa 

121 

760 

17.80 

13.90 

8.00 

95/75* 

128 

740 

17.80 

13.90 

8.00 

*(ciad/core) 


Photoluminescenoe  (PL)  measurements  were  taken  to  determine  the 
aluminum  concentration  present  in  the  core  and  dad  layers  of  the  samples. 
For  samples  562,  563,  571,  572, 1108, 1111  and  179,  PL  measurements 

weretd(anat3*K.  For  samples  121  and  128,  PL  measurements  were  taken 
at  both  300*1C  and  Selected  plots  of  these  PL  measurements  can  be 
seen  in  Appervlix  A.  At  300^,  the  Al  concentration  (x)  can  be  calculated 
by  equating  the  bencHiap  energy  at  300°K  and  the  measured  photon 

^^97  ^  where  the  peak  intensity  levels  occur  on  the  PL  curve: 
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Epj_  s  ^gpooiQ  ®  ^  .424  +  1  ^47x.  5.1 .1 

At  3*K,  X  can  be  calcuiated  by  equating  at  3°K  and 

EpL  *  »  1.518  +  1.228X.  5.1 .2 

Table  5.1.3  lists  the  Al  concentrations  determined  from  the  PL 
measurements  taken  on  the  various  samples  using  Eqs.  5.1 .1  and  5.1 .2. 
From  Table  5.1.3,  one  can  see  good  agreement  between  the  Al 
concentrations  calcuiated  from  the  PL  measurements  taken  at  the  two 
temperatures.  The  difference  between  the  two  measurements  is  within  the 
imown  accuracy  of  other  electronic  materials  characterization  techniques.^** 
The  Al  concentrations  for  samples  121  and  128  appear  very  dose  to  the 
target  core  and  dad  values.  Hence,  these  samples  were  chosen  as  the 
ones  used  to  fabricate  the  actual  waveguides  and  devices.  The  Ai 
coTKentrations  of  samples  562, 563, 571 , 572, 1 1 08, 1 1 1 1  and  179  were  far 
frx>m  the  target  and  not  acceptable  for  device  fabrication.  These  samples 
were  used  to  formulate  the  processing  techniques.  The  multiple  peaks 
measured  on  samples  562,  563,  571.  572  and  179  (see  Table  5.1.3  and 
Appendix  A)  were  most  likely  due  to  antimony  (Sb)  contamination.*^ 
Contamination  occurred  because  the  MBE  growth  machines  were  used  by 
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Number 

(“K) 

X 

X 

121 

ESI 

0.298 

0.351 

121 

3 

0.277 

0.332 

128 

Bsa 

0.296 

0.348 

128 

3 

0.274 

0.327 

562 

3 

0.361 

0.381 

563 

3 

0.324 

0.354 

571 

3 

0.356 

0.376 

572 

3 

0.360 

1108 

3 

0.225 

0.337 

1111 

3 

0.370 

179 

3 

0.242 

0.260 

0.413 


0.404 


0.437 


0.427 


0.438 


0.300 


several  researchers  growing  many  different  materials. 

Using  the  values  for  x  derived  from  the  PL  measurements,  the  refractive 
index  of  the  core  and  cladding  layers  of  samples  121  and  128  were 

calculated  with  the  models  from  Chapter  II.  Taking  into  account  the  ±  2% 
accuracy  of  the  known  Al  concentration  along  with  the  two  different  PL 

measurements  produced,  the  possible  waveguide  structures  are  given 


L* 


172 


Waveguide  Structure  121-30Q«>K-1 
(Xeowi  »  0.298  ±  0.005,  XcM,  s  0.351  ±  0.005) 
core^  =  Alo^QaQ,7Q2As  clad^  =  Aig^^Gao^^oAs 
riooni  s  3.420  ±  0.003  rv^i  -  3*388  ±  0.003 


Waveguide  Structure  121-300*»K»2 

(XooNa  -  0.313  ±  0.005,  s  0.366  ±  0.005) 

00102  *  AlojfsQagggjrAs  ciad2  =  Aioje6^^.634As 

=  3.404  ±  0.003  s  3.374  ±  0.003 


‘^veouide  Structure  121-300°K-3 


(x^  «  0.283  ±  0.005,  Xotads  =  0.336  ±  0.005) 


*  AiQ^223Qaa7i7A8  dadg  =  Alo^Gao4Q4As 

s  3.439  ±  0.004  =  3.403  ±  0.003 


^Vavegukte  Structwe  121-3°K-1 


(3^1  =  0.277  ±  0.005,  Xctodi  =  0.332  ±  0.005) 

cxxBi  »  Alo^Gao.723A8  dad^  =  Aig  332Gao.e68As 

rVoMi  s  3.433  ±  0.003  s  3.399  ±  0.003 
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Waveguide  Structure  l21-3«*K-2 

(>W  *  0202  ±  0.005,  XBtad2  =  0*3^7  ±  0.005) 

00^2  ®  ^^0JM8®®0.708^®  Ciad2  —  ^io.347^®0.8S3^® 

r\^  s  3.414  ±  0.003  =  3.385  ±  0.003 


aveaukte  Structufe 


121-3^K>3 


(Xoo«3  =  0.262  ±  0.005,  =  0.317  ±  0.005) 


00^3  ®  Aln^GawTiMiAS  CiaCl3  —  Ain  a,-yQ8<^ 

=  3.454  ±  0.004  =  3.416  ±  0.003 


*Vaveauide  Structure  128-300”K»1 


(Xo«,i  =  0.296  ±  0.005,  Xbm,  =  0.348  ±  0.005) 


corei  =  AlQ2960i^.7O4A8  dad^  »  Alna^Ga^^^^yAs 
=  3.421  ±  0.003  i\tadi  =  3.390  ±  0.003 


•Waveguide  Structure  128-300°K-2 
(JW  =  0.31 1  ±  0.005,  Xetodz  =  0.363  ±  0.005) 
core2  =  AlQ_3ffGao,3aBAs  ciad2  = 


^001.8  =  3.405  ±  0.003  r\^  =  3.376  ±  0.003 
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Waveguide  ^cture  128-300*K-3 

(}<oon3  ^  0^81  ±  0.005,  s  0.333  ±  0.005) 

core3  *  Alo^^Gaoj^gAs  dadg  —  ^io.333^^.667^^ 
s  3.441  ±  0.004  =  3.405  ±  0.003 

Waveflyjd^  i.g9-9°K-i 

(Xoo»i  =  0.274  ±  0.005,  Xctadi  =  0.327  ±  0.005) 

oore^  as  Alo^4Gao,72eAs  ciad^  =  Alo.327Gao673^ 
I^ooni  »  3.434  ±  0.003  as  3.402  ±  0.003 

Wa>fflWK!a-SlrMdvr^J.^g:a!K-g! 

(^nne  0.289  ±  0.005,  x^  3  0.342  ±  0.005) 

coje2  =  Alo^GdojiiAs  dad2  =  Aig 342680.668^ 

3  3.416  ±  0.003  ru^  =  3.388  ±  0.003 

(Xoo»3  =  0.259  ±  0.005,  Xa^  =  0.312  ±  0.005) 
ooTBa  =  Aio^Gao744As  ciadg  =  ^>412680.688^8 
Heo^a  3  3.456  ±  0.004  3  3.41 9  ±  0.003 
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The  measured  thickness  of  the  core  and  cladding  layers  of  samples  121 
and  128  are  given  in  Table  5.1 .4.  They  were  inferred  during  MOCVD  growth 
based  on  the  growth  rates.^" 


Table  5.1.4 

Core  and  Clad  Thickness  for  AlxGa^^xAs 
Samples  121  and  128 


Top 

Bottom 

Sample 

core 

Clad 

Clad 

Number 

(urn) 

(urn) 

(urn) 

121 

0.90  ±0.04 

0.95  ±0.04 

3.95  ±0.04 

128 

0.90  ±  0.04 

0.95  ±0.04 

3.95  ±0.04 

PgtsnmingfiQn  gf 

The  EO  coefficient  is  generally  defined  as  the  undamped  or  constant 
stress  EO  coeffidant  (r^4i).^  r\  is  a  value  that  can  and  has  been 
measured.  It  is  referred  to  as  the  undamped  coefficient  because  it  is  often 
measured  at  low  frequencies  in  comparison  to  acoustic  resonator 
frequencies.  r\^  is  broken  down  into  two  parts:  a  constant  strain  or 
damped  value  (r®4,)  and  a  piezoeiedric  induced  value  (rp^j:^22.i24.i34.i37 


af®  +rP 
'  41  ®  '  41  +  '  41* 


5.2.1 
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r^4i>  ^41  find  r  4,  are  given  in  10'^^  nW.  In  particular,  r\^  for  GaAs  has  been 
measured  at  several  different  wavelengths.^'*’'“’'“'^*®’^“’^®  ’®^ 

The  linear  photoelastic  constant  (o^),  given  in  10'^^  cm^/dyne,  is  a 
measure  of  the  internal  stress  due  to  epitaxial  growth.  It  has  been  recorded 
for  GaAa  at  several  wavelengths.^^'^^'^^  This  provides  the  data  necessary 
to  model  r\^,  arfo  for  GaAs  as  a  function  of 

wavelength.^^*^*^’^^®'^®'^*''®''^^  0^44  has  also  been  measured  for 

A^4aQfioj2A8-^^^  This  allows  r  44  to  be  modelled  for  AlxGai.xAs  as  a  function 
of  both  wavelength  and  Al  concentratfon.^’^’^^  Since  r\^  has  been 
measured  for  A^47Qao.a3As/^  this  allows  us  to  model  for  AI^Ga^^  as 
a  function  of  both  X  and  x  and,  hence,  determine  r\^  for  AlxGai.xAs  as  a 
function  of  both  X  and  x. 

We  will  start  with  a  model  for  for  GaAs  presented  by  Higginbotham, 
Cardona  and  PoNak^^  and  Adachi.^^’^^'’^  The  photoelastic  effect  can  be 
expressed  as^-^*' 


0^  *  Atf/X  a 
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where  X  is  the  stress  in  10^^  dyne/cm?,  e, »  Ej  Is  the  component  of  the 
dielectric  constant  in  the  absence  of  stress  (see  Chapter  II),  P||,  is  the 
component  of  the  fourth-rank  photoelastic  tensor  (constants)  and  is  the 
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component  of  the  elastic  compliance  tensor  in  10'^^  cm^/dyne.  aE|  is  the 
change  In  the  real  part  of  the  dielectric  constants  parallel  and  perpendicular 
to  the  stress  given 


Ae|  =  e^ii  - 


5^.3 


The  linear  photoeiastic  tensor  in  a  cubic  43m  class  crystal  can  be  written 
as^“ 


P11 

Pl2 

Pl2 

0 

0 

0 

Pl2 

Pl1 

Pl2 

0 

0 

0 

Pl2 

Pl2 

Pl1 

0 

0 

0 

0 

0 

0 

p^ 

0 

0 

0 

0 

0 

0 

P44 
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0 

0 

0 

0 

0 

P44 

5^.4 


Similarly,  the  elastic  compliance  tensor  in  a  cubic  43m  class  crystal  can 
100 


be  written  as 
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S., 

Sb 

S« 

0 

0 

0 

S„ 

S.1 

0 

0 

0 

S„ 

S„ 

S.. 

0 

0 

0 

0 

0 

0 

S44 

0 

0 

0 

0 

0 

0 

S44 

0 

0 

0 

0 

0 

0 

S44 

5J2.5 


From  Eqs.  4.3.18, 5.2.4  and  5.2.5,  one  can  see  that  and  S44  are  the  only 
components  that  will  contribute  to  the  EO  coeffictents  and  r\).  We  can 
write^** 


5.2.6 


where  e,4  s  S44X14  is  the  strain  component  connected  with  the  stress.^^  Pi^ 
can  also  be  defined  using  the  inverse  dielectric  oonstant^^ 


a(1/B^(C0))|  s 'AEi/fEiSj)  sr  £p||^0|l^  5.2.7 

where  8,(0))  is  the  real  part  of  the  dielectric  constant  given  by  Eq.  2.2.3. 
Substitution  of  Eqs.  5.2.4  -  5.2.7  into  Eq.  5.2.2  leads  to^^ 
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o^44(a^  ■  -e,(oi))V44S44. 


The  first-order  change  in  due  to  applied  stress  is  given  by' 


5^.8 


s  £[(d£i/dM^AM/ +  (de^ydEg^EgJ, 


5^.9 


where  is  the  interband  transition  or  gap  energy  of  the  /th  transition  (i.e., 
Eo,  Eo  Ao.  E),  etc.)  and  M/  is  a  strength  parameter  used  to  fit  the  measured 
data.  Substituting  Eq.  2.2.3  into  Eq.  5.2.9,  Adachi  obtained  the  expression 
for  the  photoelastic  coefficient:^ 


*  Co{-g(x)  +  (4EyAo)Ef  (X)  •  (E</(Eo  +  Ao))^ 
•/(XJD  +  D,, 


5.2.10 


where 


g(x)»(x)^2-(1+xr'®-(1-X)n 


5.2.11 


Eb»  Ee  ***  /(x)  ^  /(Xw)  ^  determined  using  Eqs.  2.2.4  •  2.2.1 0.  Co  and 

Do  are  sbengih  parameters  found  by  fitting  Eq.  5.2.10  to  the  elasto-optic 
measurements  recorded  by  Higginbotham  et  ai.  for  GaAs.'^^  With  the 
measurements  taken  by  van  derZiel  et  al.'°'  for  Alo  ^sGao^As,  we  have  been 
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able  to  determine  and  as  a  function  of  Al  concentratbn  Doing 
that  yields 

Co  =  -0J21  +0.19X  5.2.12 

and 

Do  =  2.03-1.95x.  5.2.13 

where  C,  arxl  are  in  lOr^^  cmVdyne.  With  $44  given  by" 

S44  s  (1 .68  +  0.02x)  X 10'“  cmVdyne  5.2.14 

and  0,4  given  by" 

0,4  =  -(0.160  -  0.065X)  C/m*,  5.2.15 

we  can  calcuiale  r^4,  as  a  function  of  both  wavelength  and  Al  concentration. 
Figs.  5.2.1  and  5.2.2  show  plots  of  1^4,  versus  x  (at  A.  s  827  nm)  and  X  (for 
X  s  0.300),  respectively. 

Using  Eq.  4.3.5  and  the  definition  for  the  refractive  index,  the  linear  EO 
effect  can  be  defined  by 


0^0 


Figure  52.2 

Piezoelectric  C^fficient  for  AlxGat.xAs 
versus  Wavelength  for  x  »  0.300 
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A  (1/e)| «  f%^ir 

Substituting  Eq.  5.2.7  into  Eq.  5.2.16  yieids^^ 


5.2.16 


r®,  s:A(1/ei)/E=  -Ae,(1/e^^.  5.2.17 

it  should  be  noted  that  as^  in  Eq.  5.2.17,  used  to  calculate  r®4f,  is  not  the 
same  one  used  to  calculate  ae^,  used  to  calculate  r%^,  defined 
hereafter  as  e^®.  can  be  obtained  by  taking  into  account  the  changes  in  the 
lowest  direct  gap  energtes:^^'®' 

AE,®  »  0e^E>Eo  + 

+  +  {^E,/(aE^AJ](AE^  +  •  •  •,  5.2.18 

where  E^  is  given  by  Eq.  2.2.9,  is  given  by  Eq.  2.2.11  and  aEq  and  aAq 
are  the  electric  field  induced  changes  in  E^  and  A^,  respectively. 
Considering  only  the  linear  EO  effects  aliows  us  to  neglect  the  higher  order 
terms.  Eq.  5.2.18  can  be  reduced  to'^^®' 


E,®  =  0eV3Eo)  Eo  +  OsVa/g  A,. 


5.2.19 
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Adachi  conveyed  that  the  changes  in  the  band  parameters  can  be  written  in 
terms  of  the  first-order  Stark-iike  effect  as'^'^ 


5^.21 

where  a^  and  b^  are  strength  parameters.  Substituting  Eqs.  5.2.20  and 
into  Eq.  55.19  yields’” 

Ae,*(<i>)e.  >  (3s^E>E.  =  (a^XA/EJEsOl)  5552 

and 

Ae,®(©)^  «  Os^AJiA,  =  bi^/(x).  5.223 

where  g(x)  and  fix)  are  defined  by  Eqs.  5.2.11  and  2.2.4,  respectively. 
Substituting  Eqs.  52.22  and  5223  into  Eq.  5.2.17  gives^^'^^ 

r*4i(<»)  =  -[1/(Ci(«>)rt[(a,/2)[A/Eo]g(X)  +  b,/(x)  +  FJ.  5224 


aEo  =  a,E 
and 

aAo  =  b,£, 


where  Fo  is  a  strength  parameter  that  accounts  for  the  nondispersive,  higher 
gap  energies.  Adachi  pointed  out^®'  that  g(x)  shows  a  very  sharp  dispersion 
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near  compared  with  /(x).  Therefore,  Aei^(o>)Ao  can  be  included  in  the 
contribution  of  the  higher  gap  energies.  Doing  this  yieids^^^ 

r’,(a>}  =  ^1/(e,(B)))“JE,•g(x)  +  5555 

with 

E;  =  (a^/2)[A/EJ.  5.2.26 

The  parameters  Eo*  and  were  adjusted  to  fit  the  experimental  data.  r\^ 
for  GaAs,  measured  by  Faist  and  Reinhart,'^*  and  r\^  for  Aio.^7Gao.a3As, 
measured  by  Glick,  Reinhart  and  Martin, were  the  coefficients  chosen  to 
fit  Eq.  5.2.25.  r\  measured  by  Faist  et  al.  was  used  (instead  of  those 
measured  by  Sugie  and  Tada^^  and  Suzuki  and  Tada^^  primarily  because 
the  Sugie  et  al.  and  Suzuki  et  al.  measurements  were  taken  on  heavily  Cr- 
doped  QaAs  crystals.  The  Faist  et  al.  measurements  were  taken  on 
undoped  QaAs  grown  on  a  Si-doped  GaAs  substrate.  The  AIq  .,7Gao.e3As 
layer  investigated  by  Qlick  et  al.  was  also  grown  on  a  Si-doped  GaAs 
substrate.  In  addition,  one  of  the  researchers  who  measured  for  GaAs 
also  measured  for  Aio,^7Gao.s3As.  We  are  looking  for  comparisons 
between  the  EO  coefficient  for  QaAs  and  AlxGai.xAs.  it  is  believed  that  more 
consistent  results  could  be  achieved  using  measurements  taken  on  both 
GaAs  and  A4Gai.xAs  by  the  same  researcher,  using  substrates  doped  with 


186 


the  same  material.  Substituting  the  measured  data,  =  -1 .68  x  10  '^  m/V 
at  X  s  1.15  ^m  for  X  s  0.00,  r\^  ®  -1.72  x  10’’^  m/V  at  X  =  1.09  urn  for  x  * 
0.00  and  =  -1.43  x  10'^*  rrW  at  X  =  1.1523  urn  for  x  =  0.17,  into  Eqs. 
5.2.1  and  5.2.25,  the  strength  parameters  were  found  to  be 

Eo‘  =  -1.375[Ao/EJ  5.2.27 

and 

F;a18.15-38.5x,  5.2.28 


with  Eo*  and  f;  in  10  '*  m/V. 

The  model  presented  above  is  valid  only  for  those  energies  below  the 
direct  band  edge.  Therefore,  we  will  limit  the  Ai  concentration  for  this  model 
to  include  only  those  values  that  ensure  that  we  are  below  the  band  edge. 

For  this  experiment,  where  X  »  827  nm,  x  must  be  greater  than  0.10  in  order 
to  use  this  model.  Fig.  5.2.3  plots  and  versus  x  at  X  s  827  nm.  Fig. 

5.2.4  plots  and  versus  X  for  x  =  0.300.  The  PL  measured  x  values 
for  the  core  layers  of  Samples  121  and  128  were  substituted  into  the  model 
to  determine  the  linear  EO  coefficient  to  be  used  for  this  experiment.  The 
r  41  values  at  X  «  827  nm  using  these  four  Ai  concentrations  are  given  in 

Table  5.2.1.  Also  included  in  the  table  is  the  estimated  ±  5%  accuracy  for 
the  measured  EO  coefficients."*''^  Table  5.2.1  provides  the  values  of  the 
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Table  5J2A 

Linear  ElectroOptic  Coefficients  for  AlxQai.xAs 
Samples  121  and  128  at  X  s  827  nm 


PL 

Sample  Temperature 
Number  (°K) 


'■»« 

(X  10“  nW) 


121 

300 

298 

-1.49  ±0.08 

121 

3 

277 

-1.56  ±0.08 

128 

300 

296 

-1.50  ±0.08 

128 

3 

274 

-1.57  ±0.08 

linear  EO  ooefficfent  that  can  be  used  in  Eqs.  4.5.4  and  4.5.5  to  calculate 
the  phase  change. 


This  section  details  the  techniques  used  to  fabricate  the 
waveguides  and  directional  couplers  for  this  experiment  as  well  as  the 
testing  methods  used  to  characterize  the  devices. 


Straight  and  Crossthrough  Waveguides 
Integrated  AlxGa^^K^s  optical  waveguides  and  crossthroughs  were 
fabricated  in  sample  121  using  processes  compatible  to  those  for 
electronics.  The  mask  used  to  pattern  the  waveguide  devices  was  a  clear 
field  type  designed  using  an  Omni-Dex  Design  Workshop  2000  graphics 


programming  environment  CAD  system.  A  chromeKm-glass  mask  was 
generated  from  this  iayout  using  a  Jeol  JBX>50il  eiectron  beam  fine 
lithography  system.  This  type  mask  was  used  to  pattern  the  waveguides 
and  devices  in  positive  photoresist  deposited  on  top  of  the  wafer  using  a 
Karl  Suss  MJB-3  mask  aligner  and  505  optical  energy  controlk  The 
waveguides  were  then  formed  by  etching  the  material  not  protected  with 
photoresist  using  a  Plasma-Therm  Series  700  reactive  ion  etch  (RIE) 
machine.  RIE  is  a  dry  etch  technique  that  uses  directional  ion  bombardmen: 
of  the  material  to  be  etched.^"  The  reactive  species  generated  from  a 
plasma  are  used  to  etch  the  material.^"  RIE  is  anisotropic  (i.e.,  the  vertical 
etch  rate  is  much  faster  than  the  lateral  rate).^^  This  leads  to  very  straight 
vertical  sidawails.  The  step-by-step  procedures  for  the  waveguide 
febifoation  are  outlined  in  Appendix  B. 

Thetarget  width  of  the  straight  and  crossthrough  waveguides  was  3  pm. 
and  the  target  etch  depth  was  2  pm  (see  Fig.  5.1 .3).  The  actual  parameters 
realized  after  fabrication  were  measured  using  a  Cambridge  Stereoscan  250 

Mk2  scanning  electron  microscope  (SEM),  a  Jeol  JSM-IC845A  SEM  and  a 
Tencor  Aipha-Slep  250  surfece  profiler.  The  actual  width  of  the  waveguides 

measured  2.7  pm  -  2.8  pm  (see  Fig  5.3.1).  The  actual  etch  depth  measured 
1.9  pm  -  2.1  pm.  Fig.  5.3.2  illustrates  the  sidewall  morphology  that  was 
achieved  using  RIE.  The  waveguides  were  characterized  using  a  Cohu 


Figure  5.3.1 

SEM  Aerial  View  of  an  AlxQai.xAs  Straight  Waveguide 
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6315  CCD  video  camera  and  a  Data  Translation  DT3852  frame  processor. 
The  827  nm  Sharp  LT015M0O  laser  diode  was  collimated  using  a  14.5  mm 
focal  length  (FL)  lens  and  passed  through  a  3X  anamorphic  prism  pair  to 
transform  the  elliptical  beam  to  circular.  It  was  focused  into  the  waveguide 

with  an  8.0  mm  FL  lens.  This  produced  a  1  pm  spot  size  at  the  waveguide 
input.  The  numerical  aperture  (NA)  of  the  8.0  mm  FL  lens  was  0.50,  giving 

a  possible  coupling  efficiency  of  84%.^**  The  beam  was  viewed  with  the 
CCD  camera  from  above  as  H  propagated  through  the  waveguide.  The 
image  was  captured  using  a  frame  processor.  Fig.  5.3.3  is  a  photograph  of 
the  waveguide  characterization  system.  A  captured  image  of  the  laser 
propagating  through  a  straight  waveguide  using  this  system  can  be  seen  in 
Fig.  5.3.4.  Images  of  the  propagating  signal  were  captured  at  the  input, 
center  and  output  of  the  waveguide.  The  average  pixel  intensities  of  these 
images  were  then  analyzed  with  Data  Translation  GLOBAL  LAB  image 
processing  software.  Since  the  coordinates  of  each  image  were  known,  the 
pixel  intensities  could  be  translated  into  propagation  loss.'”  This  frame 
processor  system  was  capable  of  detecting  256  shades  of  gray.  The  field 

of  view  of  the  CCD  camera  was  700  pm.  Therefore,  losses  as  small  as  0.2 
dB/cm  could  be  detected.  Using  this  method,  the  typical  propagation  loss 

of  the  straight  waveguides  measured  between  0.5  dB/cm  and  1.0  dB/cm. 
The  coupling  efficiency  for  end  fire  coupling  of  the  laser  into  the  waveguides 


Figure  5.3.4 

CCD  Image  of  Laser  Beam  Propagating  through 
an  AlyGa^^  Straight  Waveguide 
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was  also  measuiad.  It  was  found  that  a  75%  coupling  efficiency  coukJ  be 
attained  with  careful  alignment. 

Waveguide  crossthroughs  were  tested  using  the  methods  described 
above.  Fig.  5.3.5  is  an  aerial  photomicrograph  of  a  crossthrough  taken 
using  a  Carl  Zeiss  Axioplan  MC80  microscope.  Fig.  5.3.6  is  a  captured  CCD 
image  of  a  crossthrough  with  the  laser  beam  propagating  in  one  of  the 
waveguides.  The  second  waveguide  was  normal  to  the  first.  Comparing 
light  intensities  in  both  waveguides,  no  crossover  from  the  first  waveguide 
to  the  second  was  observed.  Wtth  the  limitation  of  gray  scale,  an  exact 
measure  of  the  crossover  could  not  be  performed  using  the  DT3852  frame 
processor.  The  minimum  signal  that  could  be  detected  was -24  dB.  Since 
no  intensify  in  the  perpendicular  waveguide  was  detected,  the  maximum 
crossover  for  these  crossthrough  structures  was  ^  -24  dB. 

Passive  Zero-Gap  Directional  Couplers 
Passive  (minus  switching  electrodes)  AI^Ga^xAs  zero-gap  directional 
couplers  CZGDC’s)  were  fabricated  in  sample  121  using  the  same  techniques 
presented  in  the  preceding  section  and  Appendix  B.  The  target  width  of  the 

input  and  output  guides  was  3  pm,  while  the  target  width  of  the  interaction 
section  was  6  pm.  To  experimerTtally  determine  the  interaction  length  for 
complete  coupling  of  the  ZGCCs,  several  structures  were  fabricated  with 
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various  interaction  channel  lengths.  The  interaction  chanriel  lengths  ranged 

from  100  ^m  •  1000  in  increments  of  100  ^m.  Fig.  5.4.1  is  a  schematic 
of  the  section  of  the  mask  for  a  passive  ZGDC  with  a  200  pm  interaction 
channel  length.  As  can  be  seen  by  the  dimensions,  a  wavefront  tilt  angle 

of  a  s  0.91**  was  chosen  (see  Eq.  1.6.1).  Fig.  5.4.2  is  a  photomicrograph 
of  f^ricated  passive  ZQDC’s  with  200  pm  interaction  channel  lengths  taken 
with  the  Zeiss  Axioplan.  The  actual  width  of  the  input  and  output  channels 

of  the  devices  measured  2.7  pm  •  2.8  pm  (see  Fig.  5.4.3).  The  actual  width 
of  the  interaction  channel  measured  5.5  pm  -5.6  pm  (see  Fig.  5.4.4).  The 
etch  depth  measured  1 .9  pm  -  2.1  pm.  The  processing  equipment  used  for 
this  experiment  could  not  resolve  separations  less  than  0.5  pm.  The  two 
input  channels  joined  together  below  0.5  pm  of  separation.  The  same  was 
true  for  the  two  output  charmeis.  The  channels  did  not  converge  to  a  point. 

This  added  60  pm  to  the  length  of  the  interaction  channel.  To  account  for 
this,  the  length  of  the  interaction  channels  of  the  ZGDC’s  on  the  mask  were 

reduced  by  60  pm.  The  actual  lengths  were  verified  by  measuring  the 
separation  between  Input  Channel  1  and  Input  Channel  2  where  they  Joined 

the  interaction  channel  (see  Fig.  5.4.5). 

To  characterize  the  structures,  the  laser  diode  was  launched  into  one  of 
the  input  channels  of  the  couplers  and  the  intensity  of  the  beam  was 
measured  at  both  output  channels.  Figs.  5.4.6  and  5.4.7  are  captured 
images  of  the  laser  beam  propagating  through  passive  ZGDC's  with  200  pm 


Figure  5.42. 

Photomicrograph  Aerial  View  of  Passive  Ai^Qa^xAs  Zero-Gap 
DirecUonal  Couplers  with  200  pm  interaction  Lengths 


Rgure  5.4.5 

SEM  Aerial  View  of  the  Intereectlon  of  the  Input  Channels 
and  the  Inleraetion  Channel  of  a  Passive  AI^Gai^ 
ZerchQap  Directional  Coupler 
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and  400  ^  interaction  channel  lengths,  respectively.  For  the  ZQDC  with 
a  200  pm  interaction  channel  length,  the  beam  enters  Input  Channel  2  and 
exits  Output  Channel  1 .  For  the  ZGDC  with  a  400  pm  interaction  channel 
lenglh.  the  beam  enters  Input  Channel  1  and  exits  Output  Channel  1. 

Output  intensity  measurements  were  taken  of  aH  the  various  interaction 
channel  length  passive  ZQDC’s.  These  measurements  were  then 
normalized  and  fit  to  the  sin?(Oo^)  cunre  corresponding  to  Eq.  1 .5.5  by 
varying  L  By  doing  this,  it  was  observed  that  an  optimum  fit  occurred  at  an 
irfieraction  length  of  L  s  196  pm  ±  3  pm  (see  Fig.  5.4.8).  This 
eiqperimentaiiy  inferred  interaction  length  for  complete  coupling  wiH  be 
compared  to  the  model  in  the  next  chapter.  In  addition  to  determining  the 
interaction  length,  the  total  device  loss  was  measured  by  comparing  the 
Input  and  output  Intensities  of  the  various  passive  ZGDC’s.  Including  the 
waveguide  propagation  loss,  the  total  loss  of  these  devices  was  found  to 
range  between  -3  dB  and  -7  dB. 

Active  Zero-Gap  Directional  Couplers 
Active  ZGDCe  (switching  electrodes  present)  were  fabricated  in  sample 
121  using  the  process  presented  in  Appendix  B.  The  target  width  of  the 

input  and  output  channels  was  3  pm,  and  the  target  width  of  the  center 
Interaction  channel  was  6  pm.  An  initial  calculation  using  the  data  from  the 


Figure  5.4.6 

CCO  image  of  Laser  Beam  Propagating  through  a  Passive  AI^Ga^^ 
2efO-<^  Direotionai  Coupler  a  200  pm  Interaction  Length 


Figure  5.4.7 

CCD  Image  of  Laser  Beam  Propagating  through  a  Passive  AI^Qai^ 
ZSro^^  Directional  Coupler  a  400  pm  Interaction  Length 
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Interaction  Channel  Length  Oun) 
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Figure  5.4.8 

d  Intensity  Plotted  Agednst  Sin^(<>oL)  with  Z.  =  196  versus 
Interaction  Channel  Length  for  Passive  AlxGai.xAs 
Zero>Gap  Directionai  Coupler 
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pfwiout  sections  suggested  that  a  very  high  switching  voltage  would  be 

required  ifa  196  ^m  interaction  length  was  used.  To  keep  voltages  low  and 
prevent  device  breakdown,  the  length  of  the  interaction  channel  was 

extended  to  5000  pm.  in  addition,  the  placement  of  the  voltage  and  ground 
electrodes  had  to  be  different  from  that  for  the  EO  rrKxtulator  presented  in 

Chapter  iV,  or  for  the  active  ZGDC*s  presented  in  Chapter  I.  This  was  due 
to  the  fact  that  the  electrodes  could  not  be  placed  on  the  top  arKt  the  bottom 
of  the  core  layer  because  of  the  cladding  layers.  Also,  if  the  electrodes  were 
plaoed  on  top  of  the  top  dad  and  the  bottom  of  the  wafer,  the  separation  d 
between  the  electrodes  would  be  too  large,  requiring  an  excessive  voltage 

for  modulation  even  for  a  device  5000  pm  in  length  (see  Eq.  4.5.6).  This 
would  reeull  in  device  breakdown.  The  electrodes  would  need  to  be  placed 

on  top  of  the  top  dad  over  the  interaction  channel  similar  to  that  shown  in 
Fig.  1 .3.1 .  This  scheme  would  have  to  be  modified,  however,  since  there 
was  an  air  inteiface  on  the  outside  walls  of  the  ZQDC’s  fabricated  for  this 
experiment.  It  was  decided  to  place  the  voltage  eledrode  over  the 
imeraction  channel,  as  one  would  for  an  EO  modulator,  and  the  ground 
electrode  beside  the  interaction  channel  on  top  of  a  mesa.  Fig.  5.5.1  is  a 
schematic  of  the  mask  layout  for  an  active  ZGDC  with  an  interadion  length 

of  5000  pm  showing  the  electrode  placement.  The  target  dimensions  for  the 
voltoge  eledrode  were  2  pm  wide  X  4996  pm  long  centered  above  the  6  pm 


10um 


24  um 


4999  urn 


Figure  5.5.1 

Schematic  of  Mask  Layout  for  (a)  the  Waveguide  Pattern  arfo  (b)  the 
ilechode  Palfom  for  an  Active  AI^Gai^  Zero-Gap  Directional 
Coupler  with  a  5000  pm  Interaction  Len^ 


X  5000  interaction  channel.  This  was  done  to  give  a  2  ^  space  from 
the  edge  of  the  metal  to  the  edge  of  the  waveguide,  and  to  ensure  that  the 

metal  would  not  touch  the  sides  of  the  waveguides.  A  contact  pad  for  the 

voltage  electrode  was  fabricated  with  target  dimensions  of  24  pm  x  24  pm 
centered  above  a  25  pm  x  25  pm  mesa.  In  the  last  section,  it  was  shown 
that  there  is  little  if  any  crossover  between  waveguides  normal  to  each  other, 

so  tt  was  decided  to  attach  the  voltage  electrode  to  the  pad  via  an  electrode 
strip  centered  above  a  mesa  intersectirig  the  interaction  channel  normal  to 
the  channel  (see  Fig.  5.5.1).  The  target  width  of  the  intersecting  electrode 
was  the  same  as  that  for  the  voltage  electrode,  and  the  target  width  of  the 
Intersecting  mesa  was  the  same  as  that  for  the  ZQDC  interaction  channel. 

The  target  dimensions  for  the  ground  electrode  were  24  pm  wide  x  4999  pm 
long  centered  above  a  25  pm  x  5000  pm  meea.  The  contact  pad  for  the 
voltage  electrode  and  the  ground  electrode  were  both  designed  to  be  1  pm 
narrower  than  the  mesas  that  they  were  above.  Metal  on  the  sides  of  the 

voltage  pad  and  ground  electrode  mesas  was  not  as  critical.  This  would 

give  a  lateral  electrode  separation  of  9  pm  and  prevent  any  possible  shorts 
from  metal  whiskers  that  might  be  left  after  fabrication.  The  contact  pad  for 

the  voltage  electrode  and  the  groimd  electrode  were  made  large  to  ease  the 
task  of  placement  of  the  Alessi  micropositioner  probe  tips.  The  mask 
generafod  for  the  electrodes  was  a  dark  field  type.  The  mask  generated  for 
the  ZQDC  was  a  dear  field  type.  The  metal  used  for  the  electrodes 
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oonsMod  of  a  trilayer  of  300  A  of  Titanium  (Ti)/400  A  platinum  (R)/3500  A 
of  gold  (Au)  depoaited  in  that  order  using  a  Temescal  FC-1800  electron 

beam  evaporator.  Figure  5.5.2  is  an  aerial  view  of  the  center  section  of  the 
fabricated  device  showing  the  voltage  and  ground  electrodes,  including  the 


pad  for  the  voltage  electrode  and  the  intersecting  electrode. 

The  actual  width  of  the  input  aixi  output  channels  measured  22  pm  - 
2.3  pm  (see  Fig.  5.5.3).  The  actual  width  of  the  interaction  channel 
measured  4.6  pm  •  4.7  pm  (see  Fig  5.5.4).  The  actual  width  of  the  voltage 
electrode  above  the  interaction  channel  measured  1.5  pm  centered  above 
the  interaction  channel  (see  Fig.  5.5.4).  The  etch  depth  measured  1 .9  pm  - 
2.1  pm. 

These  devices  were  tested  by  launching  the  LT01 5MDO  laser  diode  into 
one  of  the  input  channels,  applying  a  DC  voltage,  and  observing  the  two 
output  channels.  Fig.  5.5.5  is  an  image  showing  the  output  of  the  device 
with  the  laser  launched  into  Input  Channel  2  with  no  voltage  applied.  As  in 
the  last  section,  the  intensities  at  Outputs  1  and  2  were  normalized  and  fit 
to  the  8ln^(OoL).  in  doing  this,  it  was  found  that  an  optimum  fit  occurred  at 
t  X  123  pm  ±  0.1  pm.  The  fit  was  much  more  sensitive  to  small  changes 
in  L  widi  an  interaction  channel  length  of  5000  pm  than  it  was  when  the 
length  of  die  interaction  channel  was  closer  to  L. 

A  small  voltage  was  then  applied  to  the  voltage  electrode  and  increased 
gradually  to  19  VDC,  making  observations  of  any  changes  at  the  outputs. 


Figure  5.5.4 

SEM  Aerial  View  of  the  Interaction  Channel  Section  of 
Active  A^Qa^^  Zero-Gap  Directional  Coupler 
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The  intensity  at  Output  Channel  1  did  not  appear  to  change  until  the  voltage 
level  reached  3  VDC.  Above  3  VDC,  the  intensity  at  Output  Channel  1 
began  to  drop  off,  reaching  a  minimum  at  »  10  VDC.  The  intensity  at 
Output  Channel  1  stayed  at  a  minimum  up  to  V  s  16  VDC  and  then 
irwreased  slightly  up  to  Vs  19  VDC.  The  intensity  at  Output  Channel  2  did 
not  appear  to  change  until  the  voltage  level  reached  10  VDC.  Above  10 
VDC.  the  intensity  at  Output  Channel  2  began  to  increase,  reaching  a 
maximum  at  V  s  18  VDC.  Fig.  5.5.6  is  a  CCD  image  showing  the  output 
with  18  VDC  applied.  Since  the  waveguide  characterization  system  could 
rxit  compensate  for  the  polarization  rotation  induced  by  the  applied  voltage, 
the  switching  voltage  could  not  be  measured  directly.  Therefore,  the 
imeraities  observed  at  Output  Channels  1  and  2  were  fit  to  the 
corresponding  sin^  and  cos^  curves  (see  Eqs.  1 .5.5  and  1 .5.4),  taking  into 
account  the  polarization  rotation  induced  by  the  applied  voltage.  The 
switching  voltage  was  determined  for  the  TE^o  mode.  This  is  illustrated  in 

Figs.  5.5.7  and  5.5.8.  A  switching  voltage  of  10.5  VDC  ±  0.3  VDC  was 
inferred  using  this  method.  The  loss  of  the  active  ZQDC’s  measured  >3  dB 

to  -7  dB,  which  is  comparable  to  that  for  the  passive  devices.  In  Chapter  VI, 
the  measured  switching  voltage  will  be  compared  to  the  predicted  voltage 
derived  by  using  the  models. 
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0  Applied  Voltage  (VDC)  20 


Figure  5.5.7 

Sifi^(Oo^)  of  TEqo  Mode  Plotted  Against  Sini^(OoL)  and 
Polarization  Rotation  versus  Applied  Voltage 
for  Output  Channel  1 


I 


Applied  Voltage  (VDC) 
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Figure  5.5.8 

CoaFi^JL)  of  Mode  Plotted  Against  CoaF{<P^L)  and 
Polarization  Rotation  versus  Applied  Voltage 
for  Output  Channel  2 


An  AND  gate  was  realized  using  an  active  zero-gap  directional  coupler 
by  letting  Input  Channel  2  represent  AND  input  ”a*  and  the  control  electrode 
represent  AND  input  'b*  (see  Fig.  1.2.1).  Input  Channel  1  was  held  low  (no 
input  beam).  Light  into  the  waveguide  channel  was  defined  as  a  logic  1. 
Absence  of  light  was  defined  as  a  logic  0.  The  voltage  that  allowed 
complele  switching  of  the  active  ZGDC  represented  a  logic  1.  Absence  of 
voltage  represented  a  logic  0.  Fig.  5.6.1  shows  the  CCD  images  captured 
at  Output  Channels  1  and  2  with  varfcHJS  combinations  of  input.  Ifwedefine 
Output  Chanrwi  2  as  the  output  of  the  logic  function  and  let  the  presence  of 
Ight  represent  a  logic  1  and  the  absence  of  light  represent  a  logic  0  then,  as 
can  be  seen  in  Fig.  5.6.1,  a  digital  AND  function  is  realized.  In  Chapter  VI, 
more  oompKcated  digital  functions  and  reconfigurable  interconnects  using 
active  zero-gap  directional  couplers  will  be  addressed. 


CHAPTER  VI 
RESULTS 

In  this  chapter,  the  characteristics  of  the  zero-gap  directional  couplers 
measured  in  Chapter  V  are  compared  with  those  predicted  usirtg  the  models 
presented  in  Chapters  II  •  IV.  The  first  section  equates  the  interaction  length 


versus  refractive  irxiexarKl  waveguide  width.  The  second  section  compares 
the  measured  and  modelled  switching  voltage.  The  last  section  deals  with 
digital  logic  and  reconfigurable  interconriects. 

Interaction  Length 

It  has  been  shown  that  if  we  ensure  we  are  far  from  cutoff,  then  ail  of 
the  various  models  predicting  interaction  length  presented  in  Chapter  III 
show  good  agreement.^'”**  To  determine  that  we  are  far  from  cutoff,  we  will 
subedlute  the  various  core  and  dad  indices  for  sample  121 ,  the  measured 
width  for  the  input  stkI  output  channels  of  both  the  passive  and  active  zero- 
gap  directional  coupler  and  X  s  827  nm  into  Eq.  3.0.1  to  confirm  that  the 

normalized  frequency  B  ^  1 .6.  In  doing  this,  we  find  B  =  1 .6  -  2.6.  Hence, 
the  Marcatm  model*'  will  be  used  to  predict  the  interaction  length  of  the 
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ZQDCs,  and  the  beam  propagation  method  (BPM)"  will  be  used  to  model 
their  behavior. 

Passive  Zero-Gap  Directional  Coupler 
As  the  separation  between  the  two  waveguides  of  a  conventional 
directional  coupler  approaches  zero,  the  effective  index  of  the  region 
between  the  two  waveguides  approaches  the  index  of  the  core  region.  After 
substituting  the  various  measured  parameters  (i.e.,  refractive  index, 
waveguide  width  and  core  thickness)  and  driving  the  waveguide  separation 
to  zero,  the  Marcatiii  model  was  fit  to  the  measured  interaction  length. 
Choosing  a  center  index  of  ng  s  1.0  and  driving  the  gap  to  zero  resulted  in 
a  predicted  interaction  length  much  longer  than  the  measured  length. 
Choosing  rig  dose  to  the  core  index  and  driving  the  separation  to  zero, 
resulted  in  a  predicted  interaction  length  significantly  smaller  than  the 
measured  length.  The  best  fit  occurred  when  the  refractive  index  of  the 
center  region  was  set  equal  to  the  index  of  the  top  and  bottom  dad.  Fig. 
6.1 .1  is  a  plot  of  the  interaction  length  versus  the  waveguide  separation  for 
Waveguide  Structure  121-300**K-1  with  05  =  02  =  04.  For  this  plot  the  ± 
values  for  the  refractive  index  were  not  induded.  The  core  index  used  was 
n^  s  3.420,  the  top,  bottom,  and  center  clad  indices  used  were  r^  =  04  =  Og 
s  3.388  and  the  side  dad  index  used  was  Og  =  1 .0.  The  width  used  for  the 


Figure  6.1.1 


Interaction  Length  versus  Waveguide  Separation 


using  the  Marcatili  Model 
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input  and  output  channels  was  b  s  2.7  pm.  the  core  thickness  used  was  t  = 

0.9  pm  and  the  wavelength  used  was  X  =  827  nm.  Fig.  6.1 .2  is  a  BPM  plot 
of  the  amplitudes  of  the  propagating  fields  of  the  ZGDC  structure  using  the 

same  parameters  along  with  the  Marcatili  modelled  interaction  length  and  a 

0.1  pm  step  size.  This  illustrates  the  predicted  behavior  of  the  device.  As 
can  be  seen,  there  is  good  agreement  between  the  two  models.  The 

interaction  length  was  also  modelled  for  a  ZQOC  using  Waveguide  Structure 
121  •300**K-1  with  b  »  2.8  pm  and  ail  other  parameters  remaining  the  same 
as  for  the  previous  case.  The  interaction  length  modelled  with  these 

parameters  was  208  pm.  Taking  the  mean  of  the  two  values  gives  an 
interaction  length  of  199  pm  ±  10  pm.  which  is  not  far  from  the  measured 
value.  Comparisons  were  made  of  all  of  the  waveguide  structures  from 

Chapter  V  for  sample  121 ,  indudkig  the  ±  values  of  the  refractive  indices. 
Table  6.1.1  Hats  the  results  of  these  comparisons.  As  can  be  seen,  the 

modeHed  value  that  best  fits  the  measured  interaction  length  is  that  using 
for  the  Ai  concentration  determined  with  the  SOO'^K  PL  measurements. 
Most  importantly,  though,  Table  6.1.1  shows  that  a  change  in  the  width  of 

foe  input  and  output  channels  by  as  little  as  0.1  pm  can  lead  to  as  much  as 
a  20  pm  change  in  foe  interaction  length.  A  change  in  the  Ai  concentration 
by  as  little  as  0.5%  can  lead  to  as  much  as  a  20  pm  change  in  the 
interaction  length.  This  could  result  in  a  total  change  in  the  interaction 

length  of  ±  40  pm.  The  voltage  independent  phase  change  could  then  vary 
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Table  6.1.1 

Interaction  Length  versus  the  Various  Waveguide 
Structures  and  Channel  Widths 

Waveguide  X  t  b  L 


Structure  n,  ng, 04,05  03  (nm)  (urn)  (^m)  Qxm) 


121-300®K-1 

1.00 

827 

0.9 

2.7 

189  ±13 

121-300*K-1 

1.00 

827 

0.9 

2.8 

208  ±13 

121-3®K-1 

1.00 

827 

0.9 

2.7 

194  ±13 

121-3«K-1 

^COI»1 

^dadl 

1.00 

827 

0.9 

2.8 

214  ±15 

121-300®K-2 

1.00 

827 

0.9 

2.7 

183  ±13 

121-300®K-2 

OOMII2 

■^ctaeB 

1.00 

827 

0.9 

2.8 

202  ±18 

121-3‘‘K-2 

'\oiiie 

1.00 

827 

0.9 

2.7 

182  ±14 

121-3®K-2 

^oUmC 

1.00 

827 

0.9 

2.8 

201  ±  20 

121-300«K-3 

•^ooii»3 

1.00 

827 

0.9 

2.7 

199  ±15 

121-300®K-3 

^eora3 

1.00 

827 

0.9 

2.8 

219  ±17 

121-3«K-3 

^ooi«3 

*\lad3 

1.00 

827 

0.9 

2.7 

202  ±15 

121-3®K-3 

•^oor»3 

^dadS 

1.00 

827 

0.9 

2.8 

225  ±17 

by  as  much  as  ±  20%.  Even  with  these  large  changes  in  the  interaction 
length  and  phase,  however,  the  output  intensity  will  only  change  by  ±  10% 
(see  Eq.  1 .5.5).  This  demonstrates  a  fabrication  tolerant  device. 


» 

Active  ZerO'Gap  Directional  Coupler 


Before  the  switching  voltage  was  analyzed,  the  measured  interaction 
length  of  the  active  At^Ga^xAs  ZGDC  was  compared  to  the  modelled  length. 
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After  substituting  the  measured  parameters  into  the  Marcatiii  model, 
considering  only  the  accuracy  of  the  known  waveguide  width,  an  interaction 

length  of117pm±7pm  was  caicuiated.  This  compares  favorably  with  the 
measured  interaction  length  within  experimental  error.  When  taking  the 

accuracy  of  the  refractive  index  into  account,  the  interaction  length  becomes 

117  pm  ±  14  pm.  The  output  of  these  longer  devices  is  much  more 
sensitive  to  the  accuracy  of  the  interaction  length  than  it  is  for  the  passive 

devices.  If  one  could  control  the  interaction  length  even  to  within  1  pm, 
which  would  hold  the  voltage  independent  phase  change  to  within  1%,  the 

output  intensity  could  vary  by  as  much  as  45%  (see  Eq.  1 .5.5).  This  is 
much  less  fabrication  tolerant  than  the  passive  devices.  Figs.  62.1  and 
622  are  BPM  plots  illustrating  the  predicted  behavior  of  active  ZQDC’s 

having  interaction  lengths  of  123  pm  and  124  pm,  respectively,  with  5000 
pm  interaction  channels  and  no  voltage  applied.  This  was  done  in  order  to 
demonstrate  the  sensitivity  of  the  output  intensity  versus  interaction  length. 

The  step  size  used  for  these  plots  was  0.1  pm.  For  the  active  device 
fabricated  for  this  experiment,  the  measured  output  intensity  with  no  voltage 

applied  was  30%  less  than  the  optimum  with  no  voltage  applied  (see  Fig. 
5.5.7).  To  achieve  maximum  coupling  for  this  device,  the  interaction  channel 

length  would  have  to  be  40  pm  longer.  With  the  degree  of  inaccuracies  for 
this  experiment,  however,  it  would  be  very  difficult  to  predict  the  exact 

interaction  length  and  extremely  difficult  to  ensure  that  it  could  be  fabricated 
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the  same  way  each  time. 

Switching  Voltage 

The  acth/e  AlxGat.xAs  ZQDC  fabricated  for  this  experiment  required  30% 
higher  vottage  than  optimum  to  achieve  switching  (see  Fig.  5.5.7).  This  was 
due  to  the  fact  that  the  interaction  channel  length  was  not  optimum.  Based 
on  the  known  accuracies  for  the  refractive  index  and  waveguide  width,  one 
could  end  up  with  a  device,  after  fabrication,  requiring  as  much  as  twice  the 
optimum  switching  voltage.  Possible  solutions  to  this  problem  will  be 
addressed  in  Chapter  Vli. 

To  evaluate  the  actual  versus  riKxJelled  voltage  required  for  complete 
switching  of  the  active  A^Gaf.xAs  ZGDC,  Eq.  4.5.6  was  used.  Since  the 
electrodes  were  placed  on  top  of  the  device,  the  vertical  distance  between 
the  top  of  the  top  dad  to  the  bottom  of  the  core  was  substituted  for  the 
electrode  separation  d.  Using  Table  5.1 .4,  and  knowing  that  0.1  pm  of 
GaAs  was  grown  between  the  core  arfo  dad  layers,  d»  1 .95  pm  for  sample 
121 .  From  electromagnetic  theory,**  the  electric  field  adually  reaching  the 
bottom  of  the  core  layer  would  be  5%  less  than  the  electric  field  applied  to 
the  electrode.  This  was  determined  from  the  lateral  distance  between  the 
center  of  the  voltage  eledrode  and  the  center  of  ground  electrode.  Table 
6.2.1  presents  the  predided  switching  voltage  of  the  adive  AlxGai.xAs  ZGDC 


for  the  various  waveguide  structures  for  sample  121 ,  as  well  as  the  linear 
EO  coefficients  from  Table  5^.1.  Comparing  the  modelled  values  for 
switching  voltage  with  the  switching  voltage  measured  in  the  last  chapter, 
one  can  see  that  the  measured  voltage  is  approximately  greater  than  the 
predicted  voltage.  This  is  due  to  the  width  of  the  voltage  electrode  and 
position  of  the  ground  electrode.  The  voltage  electrode  was  less  than  half 
the  width  of  the  interaction  channel.  The  ground  electrode  was  on  the 
Channel  2  side  of  the  device.  Because  the  voltage  electrode  was  not  over 
the  majority  of  the  Interaction  channel  and  the  ground  electrode  was  located 
only  on  one  side  of  the  channel,  it  can  be  hypothesized  that  most  of  the 
voltage  was  concentrated  in  the  Channel  2  side  of  the  interaction  channel. 
The  electric  field  distribution  between  the  voltage  and  ground  electrode  is 
further  explained  by  electromagnetic^  arxi  transmission  line  theory***  and 
has  been  analyzed  for  the  region  between  the  gate  and  drain  electrodes  in 
GaAs  metal  semiconductor  field  effect  transistofs  (MESFETs).**'  With  an 
unequal  distribution  of  electric  field  in  the  two  regions  of  the  interaction 
channel,  the  voltage  induced  change  in  the  refractive  iixiex  would,  likewise, 
be  unequal  across  the  width  of  the  channel.  The  index  change  would  be 
greater  on  the  Channel  2  side  than  on  the  Channel  1  side  and  the  device 
would  behave  more  like  a  conventional  directional  coupler  (see  Eqs.  1 .4.1 1 
and  4.5.7).  To  test  the  plausibility  of  this  theory,  the  BPM  was  programmed 


Table  62.1 

Switching  Voltage  versus  the  Various  Waveguide 
Structures  and  Electro-Optic  Coefficients 


Waveguide 
Structure 


121-300®K-1 


J21-aO(W.1  rv«., 


121-300*K-1 


121-3®K-1 


121-3®K-1 


121-3*K-1 


121-300^-2 


121-300®K-2 


121-300*K-2 


121-3n(-2 


121-3^-2 


121-3®K-2 


121-300®K-3 


121-300®K-3 


121-300*K-3 


121-3®K-3 


121-3«K-3 


(pm)  <um)  (nm)  (pmAO 


1.95 


1.95 


1.95 


1.95  5000 


1.95 


121-3*K-3  In^  1.95 


1.95 


1.95 


1.95 


1.95 


1.95  5000- 


1.95 


1.95 


1.95 


1.95 


1.95 


1.95 


1.95 


5.72  ±  0.02 


1.49  5.41 


1.57  5.14 


5.39 


1.56  5.11  ±0.01 


4.86  ±0.01 


±0.02 


1.49  5.49 


±0.01 


1.48  5.48  ±0.02 


1.56  520  ±0. 


1.64  4.94  ±0.02 


5.62  ±  0.02 


5.32  ±0. 


5.05  ±0. 


529  ±  0. 


5.02  ±0. 


1.64  4.77 


fortwo  conditions.  The  first  condition  was  to  let  the  voltage  induced  change 
in  the  refractive  index  occur  across  the  entire  width  of  the  interaction 


channel.  The  second  condition  was  to  let  the  voltage  induced  change  in  the 
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rofractive  index  occur  only  across  the  width  of  the  Channel  2  side  of  the 
interaction  channel  (i-o-  only  in  one  half  of  the  width  of  the  channeO-  The 
step  size  used  to  calculate  the  propagating  field  amplitudes  for  t>oth  plots 

was  0.1  pm.  The  BPM  model  predicted  that  for  complete  switching,  the 
change  in  the  refractive  index  for  the  case  in  which  the  index  change 

occurred  only  on  one  half  of  the  interaction  channel  would  have  to  be 
approximately  ^  higher  than  for  the  case  in  which  the  index  change 
occurred  across  the  entire  width.  Since  the  applied  voltage  is  proportionai 
to  the  change  in  the  refractive  index,  this  predicts  that  the  voltage  would  also 
have  to  be  higher  for  the  secoTKl  case.  Dividing  the  measured  switching 
voltage  by  taking  into  account  the  5%  difference  between  the  applied 

vottage  and  actual  voltage  reaching  the  core,  one  gets  5.76  VDC  ±  0.17 
VDC.  This  is  witNn  0.35  VDC  of  the  predicted  voltage  for  Sample  121- 

300*K-1  arxt  wen  within  the  experimental  accuracies.  Figs.  6J2.3  and  6.2.4 
are  BPM  plots  illustrating  the  predicted  behavior  of  the  active  ZGDC  with  no 

voltage  applied  and  with  voltage  applied,  respectively.  An  interaction 

channel  iertgth  of  5040  pm  was  chosen  for  these  plots  because  it  allows 
complete  crossing  from  one  channel  to  the  other  with  no  voltage  applied. 

Rg.  6.2.4  was  generated  for  the  case  in  which  voltage  is  applied  across  the 

entire  width  of  the  interaction  channel.  For  both  plots  a  0.1  pm  step  size 
was  used.  By  increasing  the  width  of  the  voltage  electrode  so  that  it  covers 
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a  larger  area  of  the  interaction  channei.  the  switching  voltage  should 
approach  the  predicted  value.  This  was  verified  experimentally,  where  an 
active  A^Qa^xAs  ZGDC  was  fabricated  with  a  voltage  electrode  width  of  6 

pm  (see  Rg.  6.2.5).  Using  the  analysis  technique  described  in  Chapter  V, 
a  switching  voltage  of  7.6  VDC  ±  0.3  VDC  was  measured.  This  is  the  lowest 
achieved  switching  voltage  reported  in  the  open  literature  for  an  active 

AlxGa^^xAs  ZGDC  with  a  5000  pm  interaction  channel  length  arxl  is  within 
30%  of  the  predicted  voltage.  With  better  alignment  it  is  believed  that  the 
predicted  value  could  be  reached.  Reaching  the  predicted  voltage  would 
also  provide  TTL  voltage  level  operation  without  increasing  the  length  of  the 
device.  Based  on  the  EO  properties  of  A^Gaf.xAs,  however,  5000  pm  would 
be  the  minimum  length  the  interaction  channel  could  be  to  achieve  TTL  level 
switching  for  AlxGai.xAs  ZGDC’s. 


Presented  in  Chapter  V  is  the  first  demonstrated  AND  logic  gate 
reported  in  the  open  literature  using  an  AlxGa^.^  zero^ap  directional 
coupler,  ft  can  easily  be  seen  that  by  redefining  the  functions  of  input  and 
oufout  channels  and  cascading  devices  together  that  these  simple  logic 
gates  could  be  expanded  to  more  complicated  digital  circuits  outlined  in 
Chapter  I.  Reconfigurabie  interconnects  were  also  demonstrated  in  Chapter 
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V  with  the  successful  operation  of  the  active  AlxGat..As  ZGDC.  The  input 
signal  could  be  eiectro-opticaHy  directed  to  either  output.  With  the  proof-of- 
oonoept  demonstrations  of  these  EO  circuits,  practicai  issues  then  become 
the  most  prominent  for  these  devices. 

From  a  practical  viewpoint,  one  would  have  to  consider  the  intensity 
limits  constituting  a  logic  1  and  a  logic  0.  Intensity  ranges  could  be  assigned 
to  output  logic  levels  in  much  the  same  way  that  voltage  ranges  are 
assigned  to  TTL  voltage  logic  levels.  Loss  was  high  for  the  A^Ga^^ 
ZQOCs  and  would  have  to  be  addressed  for  more  complicated  logic  circuits 
and  optical  interconnects.  Ono  would  have  to  try  to  reduce  loss  through 
improved  processing  techniques  or  Kmit  the  number  of  devices  cascaded 
together.  Integration  of  the  devices  with  sources  and  receivers  would  also 
have  to  be  considered,  in  addition,  integration  with  electronics  would  need 
to  be  investigated  since  it  is  anticipated  that  the  final  circuit  would  be  an 
electronic/electrcK)ptic  hybrid  making  use  of  the  advantages  of  both 
techrK)tQgie8.  The  voltage  induced  polarization  rotation  would  have  to  be 
adcbessed.  One  would  either  have  to  employ  a  scheme  to  detect  the 
polarization  rotation  at  the  output  or  compensate  for  it  by  integrating  active 
polarizers  at  the  outputs  of  each  device  to  detect  only  one  polarization.  The 
sensitivity  of  the  output  versus  interaction  length  brought  up  in  the  last 
section  would  have  to  be  taken  into  account.  One  would  either  have  to  find 
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ways  to  more  accurately  control  the  interaction  length  so  the  coupling  state 
could  occur  at  no  voltage  or  one  would  have  to  apply  a  voltage  constantly 
to  put  the  device  in  the  coupling  state.  All  these  considerations,  though 
primarily  concerned  with  manufacturing,  would  need  to  be  addressed  before 
practical  circuits  could  be  realized.  Possible  solutions  to  these  issues  will 
be  detailed  in  Chapter  VII. 


CHAPTER  VII 


t 


SUMMARY  AND  CONCLUSIONS 

In  summary,  the  lowest  switching  voltage  reported  in  the  open  literature 
for  an  active  AlxGai.xAs  zero-gap  directional  coupler  (ZGDC)  with  a  5000  pm 
interaction  channel  length  was  achieved.  The  first  AND  logic  gate  reported 
in  the  open  literature  using  an  AixGa^.xAs  ZGDC  was  also  demonstrated. 
The  fact  that  the  input  signal  into  the  AlxGa^.xAs  ZGDC  was  switched  from 
one  channel  to  the  other  with  applied  voltage  demonstrated  a  reconfigurable 
optical  interconnect.  These  can  be  easily  exparfoed  to  more  complicated 
digital  logic  circuits  and  reconfigurable  interconnects  by  cascading  the 
ZGDC’s  arKi  reassigning  the  functions  of  the  input/output  channels.  The 
advantages  of  these  devices,  such  as  physical  size,  packing  density  and 
reconfigurability,  may  help  enhance  electronics  as  electronic  interconnects 
reach  their  limits.  Fig.  7.1 .1  illustrates  an  envisioned  circuK  integrating  both 
electronics  and  electro-optic  Fredkin  gates. 

The  comparisons  presented  in  the  Chapter  VI  show  excellent  agreement 
between  the  models  for  refractive  index,  interaction  length,  electro-optic 
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Figure  7.1 .1 

Schematic  of  a  Reconfigurable  Integrated  Electro-Optic  Circuit 
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coefficient  and  switching  voltage,  and  the  actual  perfbnnance  of  the  devices. 
The  data  from  Chapter  V  along  with  the  analysis  from  Chapter  VI  provides 
verification  of  the  models.  These  models  can  provide  good  device  design 
as  well  as  help  predict  device  perfomnance.  By  comparing  Tables  6.1 .1  and 
6.2.1  with  the  device  performance,  it  was  determined  that  the  best  fit 

occurred  when  both  the  Al  concentration  x,  measured  at  300**K,  and  the 
average  refractive  index  n^  were  used. 

In  order  to  approach  TTL  voltage  levels,  the  interaction  channel  length 
of  the  active  AlKGai.xAs  ZGDCs  needed  to  be  5000  pm  long.  This  was 
anywhere  from  25  -  40  times  longer  than  the  optimum  interaction  length  of 
the  passive  devices,  depending  on  the  waveguide  widths.  There  was  some 
concern  over  this  until  further  analysis  was  performed.  Using  the  models, 
it  was  found  that  even  at  these  long  interaction  channel  lengths,  the  active 
AlxGa^^  ZGDC  is  5  orders  of  magnitude  shorter  than  conventional  AI^Ga^. 
xAs  directional  couplers  would  be  with  an  air  interface.  By  substituting  the 
measured  parameters  for  the  active  device  along  with  a  waveguide 

separation  of  0.5  pm,  the  minimum  resolvable  spacing  that  could  be  realized 
with  the  processing  equipment  used  for  this  experiment,  (into  the  Marcatili 

model)  yielded  a  predicted  interaction  length  of  143  m.  Even  if  one  were 
to  deposit  the  highest  index  matching  fluid  available  (n  ^  2.31 )  between  the 
waveguides,  the  interaction  length  would  be  over  7  m  or  3  orders  of 
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magnitude  longer  than  the  ZGDC.  Regrowth  of  AlxGai.xAs  between  the 
wp^veguides  was  also  considered.  If  this  was  done,  the  separation  would 

have  to  reach  1 .1  pm  before  the  interaction  length  would  be  longer  than 
5000  pm.  This  would  decrease  the  phase  shift  sensitivity  and  still  allow  us 
to  operate  at  TTL  voltage  levels.  Without  extensive  modifications  to  the 

growth  chamber,  however,  one  would  not  be  able  to  realize  adequate 
adhesion  to  the  vertical  sidewalis,  resulting  in  a  poor  interface  between  the 
core  and  side  cladding.  This  would  increase  the  propagation  loss.  In 
addition,  pattern  alignment  and  electrode  fabrication  would  be  much  more 
difficult,  as  the  number  of  processing  steps  would  significantly  Irx^rease. 
More  oontrollabie  processing  then  becomes  an  additional  advantage  of 
AI^Ga^^  ZGDC’s  over  conventional  AlxGai.xAs  directional  couplers.  Using 
the  above  analysis,  it  becomes  evident  that  an  air  interface  is  also  very 
desirable.  Waveguides  can  be  run  very  dose  together  without  the  risk  of 
coupling.  Packing  densities  for  integrated  drcuits  could  be  greatly  increased 
by  using  integrated  optics.  Before  the  AlxGai.xAs  ZGDC  can  find  a  place  in 
electronic  circuits,  questions  about  integration  and  manufacturing  must  be 


addressed. 
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SuQQeetfons  for  (mprovement  and  Further  Research 

The  major  disadvantage  of  the  5000  iim  interaction  channel  length  of  the 
active  A^Qa^^  JQOC  is  that  the  phase  retardation  and,  hence,  the  output 

intensity  are  very  sensitive  to  the  accuracy  of  the  optimum  interaction  length. 
Any  error  in  waveguide  width  and/or  Ai  concentration  is  essentially  being 
increased  by  multiples  of  the  difference  between  ^  optimum  passive 
interaction  length  and  the  actual  active  interaction  channel  length.  For  the 
active  device,  one  might  have  to  supply  up  to  twice  the  voltage  necessary 
for  switching,  depending  on  what  the  output  phase  is  with  no  voltage.  If  one 
could  operate  the  device  at  the  optimum  interaction  length  the  output 
Intensity  error  would  be  minimized.  As  mentioned  in  Chapter  VI,  with  regard 

to  the  passive  device,  that  even  with  an  error  in  the  interaction  length  of  ± 
40  pm,  the  resulting  vaHation  in  the  output  intensity  would  only  be  10%.  For 
AIxOSi^,  however,  the  voltage  necessary  to  switch  at  these  short 

distances  would  be  well  over  100  V,  causing  device  breakdown.  The  only 
solution  to  achieve  short  devices  and  TTL  voftage  level  switching  would  be 

to  use  a  material  with  a  higher  EO  coefficient.  To  realize,  say.  a  200  pm 
device  with  5  VDC  switching,  a  material  with  an  electro-optic  figure  of  merit 

of  nV^  1 .5  X 10**  m/V  would  be  required.  This  is  25  times  higher  than  that 
for  AI^Ga,^.  Potential  candidates  that  could  be  used  to  realize  shorter 
devices  include  lithium  niobate  (LiNbOg),  aluminum  gallium  nitride  (/VIGaN) 


and  nonlinear  electronic  polymers.  UNbO,,  however,  is  not  compatible 
with  standard  electronic  fabrication  processes.  It  is  expensive  to  package, 
difRcult  to  interface  with  other  materials  and  it  has  a  much  lower  data  rate 
potential  than  QaAs  due  to  its  large  dielectric  constant.  MGaN  and 
nonlinear  electro-optic  polymers  look  very  promising  but  they  have  only  just 
recently  begun  to  be  irwestigated  for  electro-optic  devices.  There  is  still 
much  development  and  characterization  that  needs  to  be  done  before  these 
materials  can  start  to  compete  with  GaAs.  These  materials  are  currently 
being  investigated  through  several  joint  Wright  Laboratory  (WL)/Advanced 
Research  Projects  Agency  (ARPA)  contracts  as  well  as  through  a  joint  effort 
between  the  Solid  State  Electronics  Directorate  (WL/EL)  and  the  Materials 
Directorate  (WUML). 

For  the  present,  however,  one  needs  to  explore  ways  to  improve  the 
accuracies  of  material  growth  and  fabrication  processing  for  Ai^Ga^^  in 
order  to  fabricate  practical  TTL  voltage  level  active  ZGDC  switches.  For 
active  Al;tGaf^  ZGDC’s,  assuming  a±  10%  variation  in  the  output  intensity 

or  a  ±  0.5%  accuracy  in  the  phase  retardation,  one  would  need  to  be  able 
fo  control  the  interaction  length  to  within  ±  0.5  pm.  In  order  to  achieve  this 
kind  of  accuracy,  the  tolerance  for  the  refractive  index  must  be  ±  0.00006 
(which  translates  into  a  ±  0.01%  accuracy  for  the  Al  concentration)  and  the 
toterance  for  the  waveguide  widths  must  ontroiied  to  within  ±  0.02  urn. 
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Aapnes.  Quinn  and  Gregory  have  recently  investigated  ways  to  improve 
on  the  accuracy  of  AI^Ga,.,^  growth/^  They  claim  that  by  using  in-situ 
spectroeNipsometnc  monitoring,  the  compositionai  precision  of  AlxGa,.xAs  can 

be  controlled  to  better  than  ±  0.1%.  WL/EL  is  tryir^g  similar  methods  to 
improve  growth  accuracies  of  their  MBE  machines.  On  a  test  run  using  a 

rTKXfified  GEN  2  MBE,  an  AI,Gai.xAs  composition  to  within  ±  0.02%  was 
achieved.^ 

For  control  of  the  waveguide  width,  there  is  a  UV  based  lithographic  tool 
being  developed  on  a  WL/ARPA  program  that  has  the  potential  to  provide 
submieron  patterning  and  nanometer  alignment  accuracies  for  electronic 
multichip  moduiee  (MCM's)  and  flat  panel  cHspiays.  This  would  also  serve 
as  an  excellent  tool  for  the  fabrication  of  Integrated  electro-optic  circuits. 

Polarisation  rotation  is  another  important  issue  that  needs  be  addressed 
before  practical  circuits  can  be  realized.  As  mentioned  in  Chapter  VI,  using 
an  active  polarizer  at  the  output  of  the  ZGDC  might  be  a  possible  solution. 

A  probable  scheme  could  be  to  use  the  EO  properties  of  the 
material  by  placing  electrodes  over  the  output  channels  of  the  device.  For 
no'ewftching  (i.e.,  no  voltage),  the  output  would  remain  at  the  same  \ 
polarization  as  the  input.  For  swKching,  the  voltage  would  be  applied  over 
the  ou^  channel  as  well  as  over  the  interaction  channel  of  the  ZGDC. 

The  signal  would  go  through  two  polarization  rotations  and  the  polarization 
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at  the  output  of  the  device  would  be  the  same  as  that  at  the  input.  This 
scheme  suggests  that  the  output  channels  would  have  to  be  the  same  length 
as  the  Interaction  channel  in  order  to  ensure  that  the  two  sections  introduced 
the  same  amourtt  of  polarization  rotation.  The  total  length  of  the  ZQDC 
switch  would  then  be  over  1  mm.  This  is  still  orders  of  magnitude  shorter 
than  the  conventional  directional  coupler.  In  addition,  the  electrodes  for  the 
output  channels  could  be  fabricated  in  the  same  steps  as  the  voltage 
electrode  for  the  interaction  channel  and  ground  electrode.  No  additional 
fabrication  steps  would  be  required. 

Loss  is  another  parameter  that  needs  to  be  investigated.  BPM  models 
were  run  using  S43erKls  as  inputs  into  the  interaction  charmeis  of  ZGDC*s. 
The  BPM  predicted  a  significantly  bwer  Ices  for  these  devices.  This  is  due 
to  the  fact  that  S4)ends  maximize  mode  confinement,  minimize  radiative  loss 
and  optimize  wavefront  tilt.^  In^siane  right  angle  comer  bends  have  also 
been  investigaled  for  both  AlxQat.xAs  and  polymer  optical  devices.^^  The 
At^iQa^xAs  comer  bends  had  air  interfaces  and  used  total  internal  reflection 
for  operation.^  These  bernts  proved  lossy.  The  polymer 

comerbends  on  the  other  hand  were  fabricated  with  metal  deposited  at  the 
comer.  The  polymer  bends  had  negligible  loss.^  Out-of>piane  right  angle 
comer  bends  have  also  been  fabricated  in  polymers.^  Out-of-plane  comer 
bends  are  beneficial  for  integrating  the  sources  and  receivers  to  the 
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waveguides  within  EO  circuits  artd  EO  MCM’s. 

Electro-optic  transceh/ers  and  modulators  are  other  areas  of 
investigation  at  WITEL  In  order  to  utilize  EO  circuits  to  their  full  potential, 
multi-GHz  transceivers  and/or  multi-GHz  modulators  will  be  required. 
Without  a  transceiver  or  rrKxiulator  capable  of  high  data  rates,  optical 
interoonnscts  may  not  progress  much  further  than  research  devices.  Many 
of  the  opticai  waveguide  materials  have  been  characterized  out  to 
frequencies  ranging  anywhere  from  40  -  100  GHz.  This  promises  good 
potential  for  achieving  high  data  rate  EO  modulators.  The  transceiver  is  the 
remaining  device  that  requires  major  development  for  EO  circuits.  Plans  to 
develop  a  100  GHz  frequency  modulated  transceiver  are  currently  underway 
atWUEL 


The  optical  interconnect  is  beginning  to  find  acceptance  within  the 
sisctronic  community.  Prototype  2  GHz  data  rate  MCM's  using  optical 
inlsrconnecls  are  currently  being  fabricated  on  WL/ARPA  programs  and 
plans  to  use  optical  interconnects  for  600  MHz  clock  distribution  in 
supercomputers  Is  planned  for  the  near  future.  These  programs  will  provide 
the  mechanism  for  insertion  of  active  ZGOCs  to  provide  electro-optic 
reconfigurability  when  they  are  fully  developed. 


APPENDIX  A 


PHOTOLUMINESCENCE  MEASUREMENTS 

The  following  pages  are  the  photoluminescence  measurements  taken 
of  A^Qai.xAs  Samples  121  (at  3^K  and  SOO'^K),  562  (at  3**K),  1 108  (at  3**K) 

and  179  (at  These  measurements  represent  the  uniformity  of  the 
aluminum  concentration  present  in  the  AlxQa,.xA$  samples  grown  by  the 

MOCVD  and  the  three  MBE  machines  used  for  this  experiment  as  well  as 

the  difference  between  PL  measurements  taken  at  the  two  temperatures. 

To  measure  the  PL  at  3°K,  a  200  mW  Spectra  Physics  375B  dye  laser 
source  pumped  with  a  15  W  Spectra  Physics  171  argon  ion  laser  was  used 

along  with  two  2  m  Bausch  and  Lomb  Spectrograph  monochrometers  placed 

end  to  end.  To  measure  the  PL  at  OOO^'K,  a  Spectra  Physics  1 65  argon  ion 
laser  source  operating  at  100  mW  was  used  along  with  a  0.5  m  Jarrell  Ash 

82*0000  spectrometer.^” 
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Figure  k2 

PhotoluminescerKse  Measurement  of  Sample  121  taken  at  300°K 


Figure  A.5 

Photoiuminescence  Measurement  of  Sample  1108  taken  at 
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Figure  A.6 

Photoluminescence  Measurement  of  Sample  179  taken  at  3°K 
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APPENDIX  B 


PROCESSING  TECHNIQUES 

Included  here  are  the  processing  techniques  used  to  f^ricate  the 
AlxQai.xAs  optical  Fredkin  gates.  The  processing  steps  start  with  grown 
material  on  top  of  a  GaAs  substrate. 

Passive  Devices 

The  passive  devices  were  fabricated  using  the  following  procedures: 

A)  Clean  the  surface  using  a  wafer  spinner: 

Acetone  (ACE)  rinse  at  1000  rpm  for  30  sec 
Methanol  (MTH)  rinse  at  1000  rpm  for  30  sec 
Isopropyl  (ISO)  rinse  at  1000  rpm  for  30  sec 
Deionized  distilled  water  (Dl)  rinse  at  1000  rpm  for  30  sec 
Blow  dry  with  dry  nitrogen  (N2) 

Bake  at  100°C  on  a  hot  plate  for  5  min 
Cool  to  room  temperature 

B)  Deposit  Shipley  1 400-27  positive  photoresist  using  a  wafer  spinner. 
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Deposit  several  drops  of  1400-27  to  the  center  of  the  top  surface  of 
wafer  using  a  syringe 

Spin  wafer  at  3000  rpm  for  30  sec  (spinning  at  3000  rpm  the  1400- 
27  photoresist  will  be  1 .7  pm  thick) 

Bake  at  100**C  on  a  hot  plate  for  5  min 

Cool  to  room  temperature 

Remove  the  edge  bead  of  resist  using  a  cotton  swab  and  ACE 

C)  Expose  the  surface  to  the  waveguide  pattern  on  the  mask  using  a 
Karl  Suss  MJB-3  mask  aligr^r  arid  505  optical  energy  controller 
with  X  =  320  nm  source: 

(Align  the  waveguide  pattern  so  it  is  perpendicular  to  the  major  flat 
(Of^  so  that  the  length  of  the  waveguides  will  be  along  the  <01 1> 
direction) 

Expose  the  photoresist  at  20  mW/cm^  for  40  sec 

D)  Develop  the  waveguide  pattern  in  the  resist  with  a  1:5  ratio  of 
Shipley  AZ351  developer  to  Dl  using  a  wafer  spinner 

AZ351:Di  (1:5)  develop  at  500  rpm  for  35  sec 

Di  rinse  at  500  rpm  for  1  min 

Blow  dry  with  dry  N2 

Bake  at  100**C  on  a  hot  plate  for  5  min 

Cool  to  room  temperature 

E)  Remove  any  photoresist  left  in  the  waveguide  channels  with  oxygen 
(O2)  plasma  using  a  Plasma-Therm  Series  70  reactive  ion  etch 
(RIE): 


Set  Oa  at  90  seem 
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Set  pressure  at  150  mT 
Set  RF  power  at  75  W 
Set  DC  bias  at  220  V 
Operate  for  1  min 

F)  Just  before  RIE  of  waveguide  pattern; 

Dip  the  wafer  in  a  1 :10  ratio  of  Shipley  buffered  oxide  etch  (BOE) 
to  01  for  15  sec 

Dip  in  running  Di  for  1  min 

Blow  dry  with  dry  Ns 

Bake  at  lOO^’C  on  a  hot  plate  for  3  min 

Cool  to  room  temperature 

G)  Etch  the  waveguide  pattern  using  a  Piasma>Therm  Series  700  RIE; 

1)  Set  helium  (He)  flow  rate  at  45  seem 

2)  Set  boron  trichloride  (BCy  flow  rate  at  60  seem 

3)  Set  pressure  at  25  mT 

4)  Set  RF  power  at  100  W 

5)  Set  DC  bias  at  230  V 

6)  Etch  for  30  min 
Pump  Down 

7)  Turn  RF  power  off 

8)  Turn  gas  flow  off 


261 


9)  Reduce  pressure  to  5  X  10*^  ml 

10)  Hold  pump  down  parameters  for  10  min 
Repeat  steps  1  - 10  for  two  (2)  additional  times 
Repeat  steps  1  -  6 

(This  gives  a  total  etch  time  of  120  min.  The  10  min  pump  down 
after  MCh  30  min  of  etch  is  to  prevent  the  resist  from  over  heating 
from  the  RF) 

H)  Remove  the  remaining  photoresist  from  the  surface  with  oxygen 
(Oj)  plasma  using  a  Plasma>Therm  Series  70  RIE: 

Set  Oj  at  90  seem 

Set  pressure  at  100  mT 

Set  RF  power  at  75  W 

Set  DC  bias  at  240  V 

Operate  for  10  min 

Active  Devices 

The  active  devices  were  fabricated  using  the  following  procedures: 

A)  dean  the  surface  using  the  cleaning  procedure  described  in  the 
previous  section 

B)  Deposit  Shipley  PMQI  positive  photoresist  using  a  wafer  spinner. 

Deposit  several  drops  of  PMGI  to  the  center  of  the  top  surface  of 
the  wafer  using  a  syringe 

Spin  at  4000  rpm  for  30  sec  (spinning  at  4000  rpm  the  PMQI 
pMoresist  will  be  0.95  pm  thick) 
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Bake  at  270**C  on  a  hot  plate  for  5  min 
Cool  to  room  temperature 

C)  Depoeit  Shipley  1 400-27  posftive  photoresist  using  a  wafer  spinner: 

Deposit  several  drops  of  1400-27  to  the  center  of  the  top  of  the 
wafer  (over  the  PMQI)  using  a  syringe 

Spin  at  4000  rpm  for  30  sec  (spinning  at  4000  rpm  the  1400-27 
photoresist  will  be  1 .35  pm  thi^) 

Bake  at  100X  on  a  hot  plate  for  5  min 

Cool  to  room  temperature 

0)  Expose  the  1400-27  photoresist  with  the  electrode  pattern  using  a 
Karl  Suss  MJB-3  mask  aligner  and  505  optical  energy  controller 
with  X  s  320  nm  source: 

(Align  the  electrode  pattern  so  that  it  is  perpendicular  to  the  msyor 
flat  (OF)  so  that  the  length  of  the  electrodes  will  be  along  the  <01 1  > 
direction) 

Expose  the  1400-27  photoresist  at  20  mW/cm^  for  30  sec 

E)  Develop  the  electrode  pattern  in  the  1400-27  resist  with  a  1:5  ratio 
of  Shipley  AZ351  developer  to  Dl  using  a  wafer  spirmen 

AZ351  :DI  (1 :5)  develop  at  500  rpm  for  25  sec 

Dl  rinse  at  500  rpm  for  1  min 

Blow  dry  with  dry  Ns 

F)  Expose  the  developed  electrode  pattern  in  the  PMQI  using  a  deep 
ultra  violet  source  with  X  s  205  nm: 

Expoee  the  PMGI  photoresist  at  900  W  for  200  sec 

Q)  Develop  the  electrode  pattern  in  the  PMGI  resist  with  Shipley  SAL 
101  developer  using  a  wafer  spinner 
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SAL  101  develop  at  500  rpm  for  50  sec 
Dl  rinse  at  500  rpm  for  1  min 
Blow  dry  with  dry  N2 

H)  Plasma  etch  the  surface  to  remove  any  remaining  developed  resist 
from  the  electrode  channels  using  a  Plasmaiine  O2  plasma  system: 

Plasma  etch  at  30  W  RF,  1 .0  Torr  of  pressure  and  a  flow  rate  of  2 
scfm  for  6  min 

Bake  at  100°C  on  a  hot  plate  for  5  min 
Cool  to  room  temperature 

I)  Just  before  metallization: 

Dip  the  wafer  in  a  1 :1  ratio  of  Shipley  BOE  to  Dl  for  15  sec 
Dip  in  running  Dl  for  1  min 
Blow  dry  with  dry  N2 

J)  Deposit  300  A  of  titanium  (Ti),  400  A  of  platinum  (R)  and  3500  A  of 
gold  (Au),  in  that  order,  to  the  wafer  using  a  Temescal  FC-1^ 
electron  beam  evaporator 

K)  Remove  all  the  unpattemed  metal  from  the  surface: 

Dip  the  wafer  in  ACE  for  53  min  to  remove  the  1400-27 
Dip  in  ISO  for  1  min 

WHh  the  surfeoe  still  wet  place  the  wafer  on  the  wafer  spinner 
ACE  rinse  at  700  rpm  for  30  sec 
MTH  rinse  at  700  rpm  for  30  sec 
ISO  rinse  at  700  rpm  for  30  sec 
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Dl  rinse  at  700  rpm  for  30  sec 
Blow  dry  with  dry  nitrogen  (isy 

(This  performs  a  metal  liftoff.  Ail  metal  on  top  of  unexposed  resist 
is  removed.  The  electrode  pattern  has  rK)w  been  defined  on  the 
surface) 

L)  Dip  the  wafer  in  90**C  Shipley  1165  remover  for  2  min  to  remove 
thePMQI 

Dip  in  running  Dl  for  2  min 
Blow  dry  with  dry  nitrogen  (Ns) 

M)  Plasma  etch  any  remaining  resist  following  step  H  above 

N)  Perform  steps  A  arfo  B  from  the  previous  section 

O)  Expose  the  surface  of  the  wafer  with  the  waveguide  pattern  using 
a  Kart  Suss  MJB>3  mask  aligner  and  505  optical  energy  controller 
with  X  *  320  nm  source: 

(Align  the  waveguide  pattern  on  the  mask  so  that  the  interaction 
len^  sections  of  the  zero-gap  directional  couplers  line  up  with  the 
electrode  patterns  on  the  wafer) 

Expose  the  1400-27  photoresist  at  20  mW/cm^  for  40  sec 

P)  CompletethedevicefabricationusingstepsD-Hfrom the  previous 
section 


APPENDIX  C 


EQUIPMENT  LIST 

Below  is  a  list  of  the  equipment  used  during  the  course  of  this 
experiment. 


Modelling 

Sun  Sparc  Station  2  computer 

IBM  personal  computer  with  Intel  486  processor 

Material  Growth 

Varian  360  molecuiar  beam  epitaxy  (MBE) 

Varian  1.5  MBE 
Varian  2  MBE 

Aixtron  AIX  200  metal-organic  chemical  vapor  deposition  (MOCVD) 

FfiferiOStion 

Omni-Oax  Design  Workshop  2000  graphics  programming  environment 
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Jeol  JBX-5011  electron  beam  fine  lithography  system 

Solitec  wafer  spinner  « 

Karl  Suss  MBJ-3  mask  aligner 

Karl  Suss  505  optical  energy  controller 

Jerry  Bachur  &  Associates  JBA  205  nm  deep  UV  iight  source  arxi  constant 
intensity  controller  with  JBA  55  shutter  controller 

Tegai  Plasmaline  O2  plasma  system 

Temescal  FC-1800  electron  beam  evaporator 

Plasma-Thenm  Series  700  reactive  ion  Etch  (RIE) 

Plasma-Therm  Series  70  RIE 

Characterization 

SPEX  Industries  1269  Spectrometer 

3^  photoluminesoence  (Spectra  Physics  375B  dye  laser,  Spectra  Physics 

171  argon  ion  laser  and  two  2  m  Bausch  and 
Lomb  Spectrograph  monochrometers  placed  end 
to  end) 

SOO'K  photoluminescence  (Spectra  Physics  165  argon  ion  laser  and  a 

0.5  m  Jarrell  Ash  82-0000  spectrometer) 

Cambridge  Stereoscan  250  Mk2  scanning  electron  microscope  (SEM) 

A 

Jeol  JSM-iC845A  SEM 

Tenoor  Alpha-Step  250  surface  profiler 

Carl  Zeiss  Axioplan  MC80  microscope 
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Cohu  6315-2001/AL16  CCD  camera 
Data  Translation  DT  3852-2  flexible  frame  processor 
Data  Translation  GLOBAL  LAB  image  analysis  software 
Hewiett-Padcard  HP  611 5A  precision  power  supply 

Miscellaneous 

ILX  Ughtwave  LDX-3207B  precision  current  source 
ILX  Lightwave  LDT-5910  temperature  controller 
ILX  Lightwave  LDM-4412  laser  diode  mount 
Alessi  MH4  probe  micropositioners 

Klinger  Scientific  MC4  programmable  stepping  motor  controller 
Klinger  Scientific  MD4  stepping  motor  driver 

Klinger  Scientific  UZ80  PP  stepping  motor  driven  vertical  translation  stage 
Klinger  Scientific  UR60  PP  stepping  motor  driven  rotation  stage 
Klinger  Scientific  UT100  PP  stepping  motor  driven  linear  translation  stage 
Klinger  Scientific  BG50  PP  stepping  motor  driven  goniometric  cradle 
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